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1.0  INTRODUCTION 

Ships  and  offshore  platforms  are  designed  and  built  to  endure 
a  rugged  environment  typified  by  extreme  events  [1-1,  1-2] 
that  are  highly  unpredictable.  The  random  seaway  is  often 
characterized  by  amassed  probabilities  [1-2,  1-3].  This 
random  seaway  acts  on  marine  structures  that  are  designed  to 
resist  the  random  loads  by  welded  structural  geometries  with 
intersecting  structural  members.  More  often  than  not  the 
loadings  are  complex;  that  is,  they  produce  two  or  three  prin¬ 
cipal  stresses  that  may  be  nonproportional  or  whose  directions 
may  change  during  a  cycle  offloading.  Structural  details  such 
as  intersections  in  longitudinal  and  transverse  framing  and 
adjacent  butt  welds  in  ships  and  nodal  joints  (K,  T)  in 
offshore  structures  are  a  few  examples.  Complex  stress 
systems  are  also  common  at  notches  or  geometric 
discontinuities.  The  variable  loading  of  these  details 
produces  a  fatigue  response  [1-4,  1-5]  under  such  stress 
systems.  The  multiaxial  fatigue  response  is  generally 
unaccounted  for  on  a  detail  level.  Computational  techniques 
for  predicting  the  state  of  stress  in  structural  elements  have 
improved  over  recent  decades  because  of  the  increased 
availability  and  capability  of  finite  element  computer 
programs.  They  also  require  substantial  effort  in  predicting 
stresses  on  a  detail  level.  The  fatigue  prediction  methods 
used  in  practice  for  designing  ship  and  offshore  structures  do 
not  generally  provide  a  comparable  level  of  detail  because 
they  are  based  on  the  nominal  stress  field  applied  to 
boundaries  of  welded  configurations. 

1.1  PROJECT  OBJECTIVE 

The  objective  of  this  project  is  to  survey  and  review 
methodologies  for  predicting  multiaxial  fatigue  performance  of 
structural  details  pertinent  to  marine  structures.  The 
research  should  ultimately  lead  to  increased  levels  of 
reliability  in  designs  and  performance  evaluations  of  existing 
structures  and  potentially  minimize  the  need  to  conduct  full- 
scale  fatigue  tests  of  ship  and  offshore  structural  elements. 

1.2  SUMMARY 

To  accomplish  the  project  objective  a  detailed  literature 
survey  was  performed  identifying  over  600  references  on 
structural  fatigue  under  complex  loading  of  various  structural 
configurations.  The  database  is  presented  in  Appendix  A.  The 
references  include  multiaxial  fatigue  approaches  used  in 
structural  systems  such  as  nuclear  reactors,  aircraft,  gas 
turbines,  automobiles  and  heavy  moving  equipment. 

The  factors  affecting  fatigue  in  marine  structures  were 
reviewed  including  the  stress  characteristics  in  marine 
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structural  details.  Also,  existing  fatigue  design  procedures 
were  reviewed  as  a  baseline  to  judge  procedures  associated 
with  multiaxial  fatigue  approaches  and  to  identify  the  extent 
multiaxial  fatigue  response  is  considered  in  the  existing 
techniques.  Essentially,  the  existing  design  approaches  are 
based  on  structural  components  tests  where  complex  stress 
distributions  are  internal  to  the  applied  nominal  stress  field 
boundary. 

Several  basic  multiaxial  fatigue  approaches  were  identified 
including  stress  and  strain-based  approaches  where  an 
equivalent  stress  or  strain  is  correlated  to  simple  uniaxial 
test  data,  critical  plane  approaches  where  crack  initiation  is 
dependent  on  a  critical  stress  or  strain  plane  and  strain 
energy  approaches  for  both  crack  initiation  and  crack 
propagation.  These  approaches  were  reviewed  and  the 
engineering  significance  discussed  in  Section  3.0.  The 
multiaxial  fatigue  approaches  were  compared  to  test  data  for 
typical  structural  details  found  in  marine  structures 
including  a  web  frame  cutout,  center  vertical  keel  (CVK) , 
hatch  corner  and  a  butt  weld  for  ship  structures  and  K  and  T 
joints  in  offshore  structures.  These  comparisons  indicate 
that  there  are  candidate  approaches  for  predicting  multiaxial 
fatigue  response  in  marine  structures.  The  amount  of 
experimental  verification  has  been  extremely  limited  and  there 
are  a  large  number  of  variables  that  have  not  been  quantified 
for  marine  environments. 

To  support  the  evaluation  of  multiaxial  fatigue  data, 
statistical  characteristics  (bias  and  scatter)  were  identified 
and  compared  to  existing  fatigue  design  approaches  by  Munse 
[1-6],  the  American  Petroleum  Institute  (API)  [1-7]  and  the 
United  Kingdom  Department  of  Energy  (UK  DOE)  [1-8].  There  are 
indications  from  the  data  that  there  are  possible  gains  in 
reliability  to  be  realized  by  applying  multiaxial  fatigue 
approaches.  However,  these  gains  must  be  evaluated  in  context 
of  the  overall  system  reliability  and  associated  level  of 
uncertainties  and  the  additional  effort  required  to  apply 
multiaxial  fatigue  procedures  as  discussed  in  Section  5.0. 

Finally,  recommendations  are  provided  to  incorporate  the 
multiaxial  fatigue  design  approaches  into  a  fatigue  design 
procedure  for  marine  structures  including  experimental  work. 
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2.0  FATIGUE  IN  MARINE  STRUCTURES 


Fatigue  in  marine  structures  is  a  function  of  the  loading 
characteristics  and  a  given  material's  capability  to  perform 
without  developing  fatigue  cracks  in  the  operational 
environment.  These  factors  will  be  examined  and  engineering 
techniques  used  to  estimate  the  load  and  response,  and  to 
provide  a  proper  perspective  for  reviewing  and  evaluating  the 
applicability  of  multiaxial  fatigue  research. 

2.1  STRESS  CHARACTERISTICS  AND  OTHER  FACTORS  AFFECTING 
FATIGUE  IN  MARINE  STRUCTURES 

2.1.1  Loading  and  Stress  Characteristics  in  Marine  Structures 

There  are  a  number  of  multiaxial  fatigue  analysis  procedures 
for  predicting  crack  initiation  and  growth.  Each  method  is 
applicable  to  a  specific  set  of  stress,  strain,  and  strain 
energy  loading  characteristics.  Each  method  must  be  evaluated 
by  comparing  the  loading  assumed  in  the  multiaxial  fatigue 
criteria  with  the  actual  loading  on  the  structure  of  interest. 
Therefore,  it  is  beneficial  to  review  the  state  of  stress  in 
typical  marine  structural  details,  many  of  which  are  where 
fatigue  cracks  are  known  to  initiate  and  propagate.  The 
actual  stress  distribution  in  the  structural  details  varies 
depending  on  the  operational  environment  encountered;  however, 
design  generalizations  are  usually  made  to  characterize  basic 
stresses  for  ships  and  offshore  structures. 

2. 1.1.1  Ship  Details 

Sources  of  fatigue  loads  on  ships  are  summarized  in  Table  2-1. 
Global  loads  are  distributed  through  the  structure  via  plates, 
girders  and  panel  stiffeners  at  welded  structural  details. 

In  the  steel  structure  of  a  ship,  the  stress  or  strain  cycles 
are  generally  caused  by  the  seaway  and  by  changing  still  water 
bending  moments.  These  loads  produce  bending  stress  and  shear 
stress  in  the  ship's  hull  girder.  These  global  stresses  are 
illustrated  in  Figure  2-1  for  a  typical  tanker  where  vertical 
bending,  lateral  bending  and  torsional  bending  stresses 
combine  in  the  primary  structural  members.  Superimposed  on 
the  hull  girder  loads  are  local  stresses  caused  by  changes  in 
hydrostatic  pressure  and  local  loading  from  ships  cargo  or 
ballast.  As  shown  in  Figure  2-1,  the  stresses  are  plane 
stresses  within  a  thin  walled  plate  structural  member.  In  a 
transverse  plane,  bending  and  shear  stresses  are  caused  by 
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TABLE  2-1 

VARIOUS  STRESSES  IN  SHIP  STRUCTURAL  DETAILS 


Residual  Stress  -  This  includes  the  locked-in  stresses  in  a 
structural  element  which  occur  during  fabrication  and  assembly  as 
well  as  the  stresses  induced  by  the  support  of  the  ship's  owr 
structure.  The  local  stress  is  then  the  state  of  stress  that 
exists  in  the  light  ship  condition. 

Initial  Mean  Stress  -  The  still  water  bending  stress  (SWBS)  may 
be  induced  by  the  addition  of  the  deadweight  which  includes 
cargo,  fuel  and  luh^  oil,  potable  water,  stores,  crew  and 
effects,  ballast  and  light  ship  bending  stress. 

Varying  Mean  Stress  -  This  refers  to  stress  changes  due  to  fuel 
burn-off,  consumption  of  consumables,  and  change  in  ballasting 
that  affect  the  total  displacement  and  attitude  of  the  ship  and, 
consequently,  the  stresses  a  structural  element  may  experience. 

Stress  Due  to  Ship's  Own  Wave  -  This  stress  is  induced  by  the 
pressure  of  the  ship's  own  wave  system.  Methods  are  available  to 
estimate  the  speed  dependent  bending  moment  contribution  and  thus 
the  stress  contribution  from  the  ship's  own  wave  system. 

Diurnal  Thermal  Stresses  -  These  stresses  arise  from  the  thermal 
expansion  of  the  topside  in  the  day  and  contraction  during  the 
night.  The  thermal  stresses  are  also  affected  by  the  amount  and 
location  of  sun  exposure  occurring  during  daylight  hours. 

~r ow  Frequency  Wave-Induced  Stresses  -  These  stresses  are  caused 
by  the  wave  forces  on  the  hull  and  the  ship  motions  due  to  these 
forces.  These  cyclic  stresses  occur  at  the  frequency  of 
encounter  of  the  ship  with  the  wave  system.  The  level  of  stress 
experienced  is  directly  related  to  (although  not  directly 
proportional  to)  the  significant  wave  height  of  the  encountered 
seaway . 

Hiah-Freouencv  Wave-Induced  Stresses  -  These  stresses  are  induced 
by  dynamic  wave  loads  which  act  on  the  ship's  structure.  The 
most  common  are  bottom  slamming,  shipping  of  water  on  deck,  and 
flare  impact.  Dynamic  loads  produce  whipping  and  springing 
elastic  motions  of  the  hull,  typically  at  higher  frequencies  than 
the  frequency  of  wave  encounter.  The  impact-induced  stresses 
will  produce  an  initial  spike  in  the  stress  records  followed  by 
high-frequency  vibrations. 
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VARYING  LATERAL 
HYDROSTATIC  LOAD 


Figure  2-1 


GLOBAL  STRESSES  DUE  TO  COMBINED  VERTICAL  AND 
LATERAL  BENDING  AND  TORSION 


differences  in  hydrostatic  pressure  and  internal  cargo  loads 
or  ballast. 

These  stress  patterns  are  transmitted  to  structural  details. 
The  stress  patterns  in  elemental  details  vary;  however,  the 
planer  character  of  stress  remains  because  the  geometry  of  a 
typical  structure.  These  stress  patterns  will  be  investigated 
further  for  specific  details  where  fatigue  cracks  are  known  to 
exist  for  a  web  frame  cutout,  a  hatch  cover  detail,  a  CVK  and 
a  bottom  plate  butt  weld. 

The  web  frame  cutout  detail  was  chosen  as  an  example  detail 
because  of  frequent  cracking  found  in  the  web  plating.  Figure 
2-2  shows  the  location  of  the  cutout  under  consideration. 

This  detail  is  commonly  found  in  tankers.  The  state  of 
stresses  was  investigated  by  ABS  [2-1]  and  is  illustrated  in 
Figure  2-3.  The  stress  patterns  in  the  cutout  were 
investigated  by  Fricke  [2-2]  for  a  similar  cutout  in  a  bulk 
carrier.  The  constituent  loadings  are  illustrated  in  Figure 
2-4  and  the  resulting  stress  patterns  are  shown  in  Figure  2-5. 
This  cutout  detail  was  also  fatigue  tested  by  Munse  [2-3]. 

The  stress  distribution  was  measured  by  Munse  with  strain 
gauges  as  shown  in  Figure  2-6.  The  load  was  applied  as  a 
concentrated  load  between  cutouts,  however,  the  stress 
patterns  are  consistent  with  those  shown  in  the  previous 
figures.  The  constituent  stress  characteristics  in  the 
vicinity  of  the  cutout  are  biaxial.  However,  the  stress  is 
uniaxial  on  the  extreme  fiber  of  the  cutout.  This  is  because 
a  free  surface  cannot  support  stress  normal  to  the  surface. 
Shear  strains  can  exist  near  the  free  surface.  Another 
important  characteristic  is  the  stress  concentration  that 
exists  around  the  cutout.  This  concentration  produces  local 
stresses  that  can  exceed  yield  while  adjacent  biaxial  field 
stresses  are  well  below  yield.  The  behavior  has  been 
confirmed  by  stress  predictions  for  this  cutout  conducted  by 
ABS  [2-1]  and  can  result  in  strain-controlled  fatigue 
cracking.  Fricke  [2-2]  also  describes  this  strain-controlled 
phenomenon  for  a  similar  structural  detail. 

Hatch  corners  often  experience  excessive  stresses  that  lead  to 
fatigue  cracking.  Compressive  and  tensile  stresses  result 
from  longitudinal  and  lateral  bending  due  to  the  hogging  or 
sagging  condition.  Also,  torsionally  induced  hull  stresses 
are  high  and  common  in  large  hatch  ships.  Figure  2-7 
illustrates  a  hatch  opening  through  a  deck  with  associated 
stresses  present.  These  stresses  concentrate  around  the 
corners  of  the  hatch  opening.  Figure  2-8  illustrates  stress 
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1  Initiation  of  Crack  in  Flat  Bar  Stiffanar 
}  Crack  at  Traa  Edge  of  Cut-out 

3  Crack  in  Side  Shall  Mating 

4  Crack  at  Radius  of  Cut-out 


Figure  2-2.  Cutout  in  a  Tanker  Web  Frame 
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Figure  2 


-3.  Predicted  Stress  Distribution  in  a  Tanker  Web  Frame 
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Figure  2-4.  Constituent  Loadings 
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Figure  2-6.  Stress  Distribution  in  Clearance 
Cutouts  for  Fatigue  Tests 
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Stress  Distribution  in  a  Containership  Hatch  Corner 
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patterns  that  were  determined  for  a  class  of  containerships 
from  finite  element  analysis. 

The  stresses  that  act  on  an  element  of  the  deck  plate 
surrounding  a  hatch  corner  are  tensile,  compression,  and 
shear.  The  proportionality  of  principal  stress  is  variable, 
depending  on  the  status  of  mean  stresses  and  encountered 
waves.  Stresses  in  hatch  corners  were  measured  as  part  of  the 
well  known  and  well  documented  SL-7  instrumentation  program 
[2-4].  Stresses  in  excess  of  yield  were  measured  in  the 
vicinity  of  the  hatch  corner.  Fatigue  cracking  occurred  early 
in  the  ships  service  life  as  documented  by  Stambaugh  [2-5]  and 
others  [2-6]  and  shown  in  Figure  2-9.  The  fatigue  crack 
growth  is  along  a  path  perpendicular  to  the  maximum  principal 
stress  field  indicative  of  Mode  I  crack  extension.  Again,  a 
local  strain  controlled  condition  exists  where  local  stresses 
are  in  excess  of  yield  and  nominal  stresses  are  below  yield  in 
the  vicinity  of  a  stress  concentration. 

Fatigue  cracking  has  been  reported  in  longitudinal  CVK  details 
where  they  bracket  to  a  transverse  bulkhead.  Figure  2-10 
illustrates  the  global  stresses  acting  on  the  CVK  due  to 
bending  and  lateral  hydrostatic  pressure.  As  shown,  axial 
bending  stresses  and  vertical  shear  are  present  in  the  CVK. 
Local  shearing  is  present  due  to  the  longitudinal  deflection 
of  the  transverse  bulkhead  which  acts  on  the  CVK  through  the 
bracket.  Figure  2-10  shows  stresses  on  the  detail  level,  and 
stresses  present  on  an  elemental  level.  Bending  and  shear  act 
along  the  two  principal  axes  forming  a  biaxial  stress  field  in 
the  CVK  web. 

Stambaugh  [2-5]  reported  on  fatigue  cracking  in  transverse 
butt  welds  located  in  the  bottom  shell  of  a  ship.  The  fatigue 
stresses  in  the  butt  weld  result  from  primary,  secondary  and 
tertiary  loading  all  with  mean  stress  and  residual  stresses 
present  as  shown  in  Figure  2-11.  The  net  response  is  a 
biaxial  field  with  "trying  longitudinal  principal  stresses. 

In  summary,  stresses  in  ship  structure  are  uniaxial  on  the 
edge  of  plates  and  cutouts  and  biaxial  in  plate  details. 
Principal  stresses  generally  align  with  the  major  longitudinal 
axis  of  the  hull  or  in  the  direction  of  major  transverse 
frames.  Secondary  and  tertiary  stresses  are  produced  by 
lateral  hydrostatic  loadings  and  induce  biaxial  components. 
However,  the  direction  of  principal  stresses  is  highly 
variable.  Phasing  of  the  stresses  is  also  variable.  All 
variations  are  random  at  wave  encounter  frequencies.  Long¬ 
term  and  short-term  load  distributions  have  been  investigated 
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Figure  2-9.  Fracture  of  the  Forward  Hatch  Cutout  on  the  High 
Speed  Containership  During  the  Second  Winter 
Season 


FLUCTUATING  LATERAL 
HYDROSTATIC  LOAD 


Figure  2-10 

Bl— AXIAL  STRESS  DISTRIBUTION  IN  A  TANKER  CVK 


2-13 


extensively  with  little  agreement  between  researchers  as  to 
the  exact  type  of  distribution  that  will  account  for  random 
load  effects.  Mean  stresses  should  also  be  taken  into 
account.  Also,  large  residual  stresses  are  present  from  the 
weld  process,  in  many  instances  on  the  order  of  magnitude  of 
yield,  and  effect  the  mean  stress  level. 

2. 1.1. 2  Stresses  in  Offshore  Structures 

Stress  analyses  of  tubular  offshore  structures  (see  Figure 
2-12)  begin  with  a  global  analysis  of  the  jacket  as  shown  in 
Figure  2-13  and  are  generally  made  through  space  frame 
analyses.  The  nature  of  member  end  loads  changes  with 
direction  of  wave  forces,  as  shown.  Other  loads  on  offshore 
structures  are  summarized  in  Table  2-2.  The  common  joint 
stress  analysis  procedure  utilizes  parametric  equations  [2-7, 
2-8,  2-9]  to  estimate  the  hot  spot  stresses  at  selected 
locations.  The  parametric  equations  are  based  on  either 
finite  element  analysis  or  on  laboratory  testing.  Currently, 
parametric  formulations  are  available  for  X,  T,  K,  and  TK 
joints  subjected  to  simple  loadings  of  axial,  in-plane  bending 
and  out  of  plane  bending  applied  separately  to  the  joints  as 
shown  in  Figures  2-14  and  2-15.  The  principal  stress 
distribution  of  a  simple  T  joint  is  illustrated  in  Figure  2- 
16.  Stress  distributions  around  the  weld  toe  of  a  Y  joint  are 
shown  in  Figure  2-17.  The  stress  state  is  biaxial  on  the 
surface  and  triaxial  through  the  thickness  with  significant 
shear  stresses  developing  across  the  thickness  of  the  chord 
and  the  brace  at  the  intersection.  Stress  "hot  spots"  occur 
at  the  saddle  point  (±90*  to  chord  axis)  of  the  intersection 
where  stress  concentrations  of  greater  than  six  times  the 
nominal  axial  stress  in  the  brace  may  occur  depending  on  the 
geometry  of  the  joint.  Peak  stresses  which  are  higher  than 
those  at  the  intersection  occur  at  the  toe  of  the  weld. 
Principal  stresses  at  the  saddle  points  lie  within  16*  to  30* 
of  the  normal  to  the  weld  line.  At  crown  points  (0*,  180*  to 
chord  axis)  the  angle  to  the  normal  is  about  8*.  Stress 
decays  rapidly  away  from  the  weld  toe. 

The  state  of  stress  in  tubular  joints  is  biaxial  on  the 
surface.  Stresses  will  generally  be  compressive  on  one  side 
of  the  chord  brace  intersection  and  tensile  on  the  other 
depending  on  direction  of  applied  in-plane  bending  moment. 
Stress  hot  spots  will  be  located  near  or  at  crown  points  of 
the  intersection.  Stress  concentrations  vary  from  four  to 
seven  times  the  nominal  bending  stress  in  the  brace.  Maximum 
principal  stresses  at  the  hot  spots  generally  lie  within  8* 
from  normal  to  the  weld  line;  however,  the  magnitude  and 
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Figure  2-12.  Illustration  of  an  Offshore  Platform 
with  Multibrace  Connections 
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•I.  HEAO  ON  WAVE 


Figure  2-13.  Typical  Joint  Geometry /Loading  Configuration 


TABLE  2-2 


INFLUENCE  OF  VARIOUS  8TRES8E8  ON  FATIGUE 
IN  OFFSHORE  STRUCTURAL  DETAIL8 


TYPE  OF  LOADING 
Wave 


Current  and  Wind 
(Static  and 
Quasi-Static) 

Current  and  Wind 
(Dynamic) 

Deadweight 


Fabrication 


Launching/ Installation 
Live  Loads 

Transportation 

Buoyant  and  Hydrostatic 

Foundation  Movement/ 
Earthquake 

Floating  Ice  Impacts 


CONTRIBUTION 


Fatigue  Crack  Initiation  and 
Propagation 

Mean  Stress  (R  ratio)  in  Crack 
Initiation  and  Propagation 


Fatigue  Crack  Initiation  and 
Propagation 

Mean  Stress  (R  ratio)  in  Crack 
Initiation  and  Propagation 

Mean  Stress  (R  ratio)  in  Crack 
Initiation  and  Propagation 
Fatigue  Crack  Initiation 
Initial  Flaw  Size  in  Crack 
Propagation 

Crack  Initiation 

Fatigue  Crack  Initiation  and 
Propagation 

Low-Cycle  Fatigue  Crack  Initiation 
and  Propagation 

Mean  Stress  (R  ratio)  in  Crack 
Initiation  and  Propagation 
Crack  Initiation 

Low-Cycle  Crack  Initiation 


Crack  Initiation 
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Figure  2-17.  Stress  Distribution  at  a  Welded  Y-Joint  Connection  with  an  Axial  Load 


direction  vary  during  the  complete  cycle.  Out-of-plane 
bending  loads  tend  to  produce  higher  stress  than  in-plane 
loads;  however  the  location  and  magnitude  of  stress  hotspot 
and  the  direction  of  the  principal  stress  at  that  location 
vary  with  loading  and  joint  geometry  parameters. 

The  interaction  between  members  of  in-plane  joints  has  been 
investigated  by  Gulati  [2-7]  where  stress  concentration 
factors  are  provided  for  stresses  normal  to  the  weld  toe  in 
complex  KT  joints  typical  of  jacket  structures;  however, 
little  is  presented  in  the  literature  on  the  multiaxial  stress 
distributions  between  multi  planar  joints.  Rather,  loading 
interaction  effects  are  accounted  for  by  superimposing 
moments.  Stress  concentration  factors  are  then  determined. 
Typical  stress  concentration  factors  are  presented  by  Gulati 
[2-7],  Marshall  [2-8,  2-9]  and  Rodabaugh  [2-10].  Gulati  also 
discusses  the  fact  that  principal  stress  in  joints  subjected 
to  combined  loading  is  not  orthogonal  to  the  weld  lineal 
direction.  This  condition  is  shown  in  Figure  2-18  for  a 
simple  T  joint.  Furthermore,  the  direction  of  the  principal 
stress  changes  with  the  changing  magnitudes  of  constituent 
loadings.  The  effect  of  the  non-perpendicularity  to  the  weld 
lineal  direction  of  the  principal  stress  (at  the  hotspot)  on 
the  crack  trajectory  and  fatigue  life  of  an  as  welded  joint  is 
not  known  at  this  time. 

In  summary,  stress  characteristics  in  offshore  structures  are 
triaxial  within  the  tube  wall  thickness.  Biaxial  stress 
fields  exist  at  the  tube  free  surface  at  the  weld  toe  in  chord 
and  brace  intersections  with  large  shear  stresses  in  the 
trough  thickness  of  the  brace  wall.  Principal  stresses  vary 
in  proportion  with  the  random  nature  of  wave  encounter.  The 
location  of  maximum  stresses  changes  with  the  encounter  and 
passing  of  waves.  Mean  stresses  are  present  in  terms  of  dead 
loads,  current  loads  and  residual  stresses. 

2.1.2  Factors  Influencing  Fatigue  Response  (General) 

In  addition  to  the  state  of  stress,  there  are  other  factors 
that  influence  fatigue  response  that  must  be  considered  in  the 
evaluation  of  multiaxial  fatigue  analysis  and  design 
procedures.  But  first  a  brief  description  of  the  fatigue 
process  will  be  extremely  helpful  in  understanding  and 
evaluating  multiaxial  fatigue  mechanisms. 

Fatigue  cracks  initiate  in  local  slip  planes  or  in  the  plane 
of  a  discontinuity  in  welded  details  depending  on  the  size  and 
shape  of  the  defect  and  local  stress  conditions.  In  the 
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Strain  Gauge 
Locations  for 
Tested  T  Joint 
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*  Pomn*  ftutm  on  chord  tOondUmm  from  maid  tot 

*  Romm  §oupm  on  bneo  10.  tS.  SOond  100mm  from  mid  tot. 


Figure  2-18.  Principal  Stress  Distribution  in  a  T-Joint.  Note  Stress  Direction 
Relative  to  Weld  Toe 


absence  of  defects,  fatigue  cracks  tend  to  grow  in  a  plane  of 
maximum  shear  stress  range.  This  growth  is  quite  small, 
usually  on  the  order  of  several  grains.  As  cycling  continues, 
fatigue  cracks  tend  to  coalesce  and  grow  along  planes  of 
maximum  tensile  stress  range.  When  defects  are  present  the 
mode  of  crack  growth  is  more  complex.  The  two  stages  of 
fatigue  crack  growth  are  called  Stage  I  and  Stage  II  as  shown 
in  Figure  2-19.  The  stages  are  important  because  multiaxial 
fatigue  prediction  techniques  are  generally  applicable  to  one 
stage  or  the  other  and  some  attempt  to  account  for  both. 

Fatigue  life  is  influenced  by  numerous  factors.  Bea  [2-11] 
has  presented  Table  2-3  listing  the  factors  influencing 
fatigue  in  offshore  structures  indicating  the  level  of 
complexity  and  amount  of  effort  required  to  predict  fatigue 
response  in  offshore  structures.  Of  a  total  of  almost  40 
general  topics,  multiaxial  fatigue  is  included  in  1  or  2 
areas.  Other  factors  affecting  fatigue  response  in  welded 
marine  structures  are  considered  as  subcategories  to  Table 
2-3.  They  are  as  follows: 

1.  Material  properties 

a.  Base  metal 

b.  Heat  affected  zone 

c.  Weld  metal 

2 .  Stress  characteristics 

a.  Stress  gradients  in  weld  geometry 

b.  Stress  proportionality 

c.  Elastic  and  plastic  strain  relationships 

d.  Mean  and  residual  stresses 

e.  Stress  phasing 

f.  Random  loading 

g.  Stress  relieving 

3 .  Corrosion 

4.  Thickness  effects 

5 .  Flaws 

6.  Fabrication  procedures 

7.  Surface  finish 
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Figure  2-19.  Schematic  of  Stages  I  and  II 
Transcrystalline  Microscopic 
Fatigue  Crack  Growth 
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TABLE  2-3 


A  SUMMARY  OF  FACTORS  AMD  CONSIDERATIONS 
RELATED  TO  FATIGUE  IN  WELDED  JOINTS  OF  OFFSHORE  PLATFORMS 


I.  FATIGUE  BEHAVIOR  OF  WELDED  JOINTS 

1.  Definition  of  fatigue  failure  in  S-N  data 

2.  Size  effect  in  S-N  data 

3.  Effect  of  weld  profile 

4.  Effect  of  corrosion  and  cathodic  protection 

5.  Assumption  of  a  linear  model  and  lognormal  distribution 
for  N 

6.  Classification  of  joint  on  the  basis  of  geometry  rather 
than  load  pattern 

7.  Relationship  between  stress  at  joint  and  stress  used  to 
obtain  S-N  curve 

8.  Ignoring  possible  stress  endurance  in  S-N  curve 

9.  Compatibility  of  determination  of  hot  spot  stress  with 
S-N  curve 

II.  MANUFACTURING  CONSIDERATIONS 

1.  Fabrication  uncertainties 

2.  Requirements  on  weld  contours  not  met 

III.  DEFINITION  OF  THE  ENVIRONMENT 

1.  Use  of  full  scatter  diagram  of  wave  height  and  wave 
period 

2.  Variations  in  wave  period 

3.  %  occurrence  estimates 

4.  Wave  directionality 

5.  Interaction  of  waves  and  currents 

6.  Theoretical  model  used  for  ocean  waves 

IV.  HYDRODYNAMIC  LOADS  ON  STRUCTURE 

1.  Inertia  and  drag  coefficient 

2.  Directional  wave  spectra  which  accounts  for  wave 
spreading 

3 .  Marine  growth 

4.  Sheltering  effects 

V.  STRUCTURAL  RESPONSE  TO  HYDRODYNAMIC  LOADS 
1.  Assumptions  made  in  spectral  analysis 

a.  linear  response  during  transfer  function 
development 

b.  linearization  of  drag  term 

c.  at  joints 

i.  no  flexibility 

ii.  effect  of  can 

iii.  center  to  center  coordinates 
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TABLE  2-3  (continued) 


A  SUMMARY  OF  FACTORS  AMD  CONSIDERATIONS 
RELATED  TO  FATIGUE  IN  WELDED  JOINTS  OF  OFF8HORN  PLATFORMS 


d.  soil  stiffness  in  dynamic  model 

e.  damping  effects  in  structural  response 

f.  dynamic  response  not  accounted  for  in  analysis 

VI.  FATIGUE  STRESSES  AT  JOINT 

1.  Method  of  analysis  to  evaluate  stress  concentration 
factors  (SCF) 

2.  Parametric  equations  used  for  SCF 

3.  Point  at  intersection  where  failure  occurs 

VII.  FATIGUE  DAMAGE  EQUATIONS 

1.  Assumption  of  Miner's  Rule 

2.  Assumption  of  narrow  band  damage  equation  in  spectral 
approach 

3.  Assumption  of  Weibull  distribution  for  stress  ranges  in 
stress  distribution  approach 

VIII.  OTHER  CONSIDERATIONS 

1.  Errors  by  designers 

2.  Bad  judgment  during  towing  and  installation 

3.  In  service  loads 
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These  factors  should  be  kept  in  mind  in  reading  the  review  of 
multiaxial  fatigue  prediction  approaches.  Few  have  been 
addressed  adequately  if  the  intent  is  to  incorporate  the 
multiaxial  fatigue  approaches  in  design  of  welded  marine 
structures . 

As  we  will  present  later,  there  are  gains  to  be  made  in 
overall  structural  reliability  by  considering  multiaxial 
effects  if  the  designer  has  the  proper  design  tools, 
information  and  financial  resources  to  do  so. 

2.2  CUMULATIVE  DAMAGE  AND  CRACK  GROWTH  APPROACHES  FOR 
PREDICTING  FATIGUE  RESPONSE 

As  with  stress  distributions,  it  is  important  to  understand 
existing  fatigue  design  procedures  before  reviewing  multiaxial 
fatigue  research.  Most  multiaxial  fatigue  approaches  are 
extensions  of  fatigue  life  and  fracture  mechanics  approaches 
developed  for  uniaxial  loading. 

The  fatigue  life  of  a  structural  detail  is  determined  by  the 
sum  of  the  elapsed  cycles  required  to  initiate  a  fatigue  crack 
and  propagate  the  crack  from  subcritical  dimensions  to  a 
critical  size.  (Note  that  the  critical  crack  size  and 
criteria  for  failure  often  differs  from  one  set  of  data  to 
another.)  The  two  basic  fatigue  prediction  approaches  that 
are  most  widely  used  in  structural  design  include  the  Miner's 
Linear  Rule  for  fatigue  life  (crack  initiation  and  growth  of 
short  cracks)  and  the  fracture  mechanics  theory  for  crack 
growth. 

Miner's  [2-12]  approach  is  based  on  knowledge  of  the 
structural  loading  and  the  resistance  of  the  structure  in 
terms  of  stress  range  and  number  of  cycles  to  failure.  This 
method  is  developed  from  test  data  (S-N  curves)  together  with 
the  hypothesis,  that  fatigue  damage  accumulates  linearly. 
According  to  this  hypothesis,  the  total  fatigue  life  under  a 
variety  of  stress  ranges  is  the  weighted  sum  of  the  individual 
lives  at  the  various  stress  ranges,  S,  as  given  by  the  S-N 
curves,  with  each  being  weighted  according  to  the  fractional 
exposure  to  that  level  of  stress  range.  To  apply  this 
hypothesis,  the  long-term  distribution  of  stress  range  is 
replaced  by  a  stress  histogram,  consisting  of  a  convenient 
number  of  constant  amplitude  stress  range  blocks,  S,,  and  a 
corresponding  number  of  stress  cycles,  r7,.  The  constraint 
against  fatigue  fracture  is  then  expressed  in  terms  of  a 
nondimensional  damage  ratio,  rj: 
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B  rii 

z  —  <  nL 

i=l  N, 


where  6  =  number  of  stress  blocks 

rjf  =  number  of  stress  cycles  in  stress  block  i 

N,  =  number  of  cycles  to  failure  at  a  constant 

stress  range,  s, 

rjL  =  limit  damage  ratio 

The  limit  damage  ratio  t?l  depends  on  the  number  of  cycles  at 
each  stress  level. 

Theories  of  crack  propagation  associated  with  most  crack 
growth  methods  are  based  on  linear  elastic  or  elastic-plastic 
fracture  mechanics  methods. 

In  linear  elastic  fracture  mechanics,  the  stress  field  at  the 
tip  of  a  crack  is  described  in  terms  of  stress  intensity  with 
a  function  of  the  following  form: 

K  =  F»a  JK  (Ksiyin) 

where  a  -  gross  stress  (ksi) 
a  =  crack  length  (in) 

F  =  geometry  correcting  factor,  dependent  upon 
crack  and  part  geometry,  and  stress 
gradient 

K  =  stress  intensity  factor 

In  the  literature  there  exists  stress-intensity  solutions  for 
a  wide  variety  of  crack  shapes  and  loading  cases. 

The  stress  intensity  is  used  in  all  known  schemes  for 
calculating  crack-growth  rates.  For  constant-amplitude  crack 
growth,  it  has  been  shown  that  crack  growth  is  primarily  a 
function  of  AK  where  AK  =  K*.,  -  K*,,,.  That  is,  da/dn  =  f(AK). 

When  constant-amplitude  crack-growth  rate  (da/dn)  is  plotted 
against  A K  on  log/log  paper,  an  s-shaped  curve  results  as 
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shown  in  Figure  2-20.  Traditionally,  this  curve  is  broken 
into  three  regions  called  Stage  I,  II  and  III  crack  growth. 
Stage  I  crack  growth,  sometimes  called  threshold  crack  growth, 
applies  when  crack-growth  rates  are  very  low.  Small  changes 
in  AK  substantially  alter  the  crack-growth  rates.  Stage  II 
crack  growth,  sometimes  referred  to  as  steady-state  crack 
growth,  is  characterized  by  a  nearly  linear  relation  between 
log  da/dn  and  log  AK,  while  Stage  III  crack  growth  is 
characterized  by  a  rapidly  increasing  crack-growth  rate  and 
results  from  crack  instability. 

Several  models  have  been  presented  which  attempt  to  describe 
the  relation  between  da/dn  and  A  K.  The  simplest,  and  limited 
strictly  to  Stage  II  crack  growth,  is  the  Paris  equation: 

da 

—  =  C  (A  K) " 

dn 

where  C  and  m  are  material-dependent  constants  and  are 
determined  by  experiment.  Other  relations  account  for  stress 
ratio  effects  or  include  Stage  I  or  Stage  III  crack  growth  but 
are  modifications  of  this  approach. 

Both  the  fatigue  life  and  fracture  mechanics  approaches  have 
either  been  applied  or  proposed  for  use  in  marine  structures 
as  described  next. 

2.2.1  Stress  Concentration  Factor  (BCD  Approach  in  Offshore 
Structures 

The  American  Petroleum  Institute  (API)  [2-13]  and  the  American 
Welding  Society  (AWS)  [2-14]  provide  a  method  of  evaluating 
the  fatigue  life  of  offshore  platform  tubular  joints  using 
either  an  experimental  or  theoretically  determined  "hot-spot 
stress"  range  or  "hot-spot  strain"  range.  The  hot-spot 
stress/strain  is  the  maximum  stress/strain  in  a  given  tubular 
joint  due  to  a  specified  load  range,  but  does  not  include  the 
"peak"  stress/strain  which  arises  due  to  the  geometrical 
discontinuity  at  the  toe  of  the  fillet  welds.  The  hot-spot 
stress/strain  with  this  definition  can  be  measured  by  strain 
gauges  (the  peak  stress/strain  cannot  be  so  measured)  or  it 
can  be  calculated  by  finite-element  analysis;  the  model  used 
in  the  finite-element  analysis  does  not  include  the  fillet 
weld  contour  or  the  extremely  fine  element  mesh  needed  to 
determine  peak  stresses. 
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The  AWS  Code  states  that  the  X-curve  is  to  be  used  with 
"greatest  total  range  of  worst  hot-spot  stress  or  strain  on 
the  outside  surface  on  intersecting  members  at  the  toe  of  the 
weld  joining  them  —  measured  after  shakedown  in  model  or 
prototype  connection  or  calculated  with  best  available 
theory" . 

The  British  have  also  developed  an  approach  for  predicting 
fatigue  life  of  welded  details  and  offshore  structural  joints 
[2-15].  The  non-nodal  joints  are  classified  according  to  weld 
configuration  and  geometry.  The  nodal  joints  are  all  repre¬ 
sented  by  one  S-N  curve.  Unlike  the  API  code,  this  approach 
recognizes  combined  stresses  and  presents  data  in  terms  of 
principal  stresses  at  the  weld  toe.  This  approach  is  based  on 
nominal  field  stresses  and  tests  of  actual  joint  con¬ 
figurations. 

During  the  past  twenty  years,  a  substantial  amount  of  work  has 
been  conducted  and  results  published  on  hot-spot  stresses  in 
simple  tubular  joints.  This  work  includes  stresses  as 
experimentally  determined  by  strain  gauges  on  tubular  joint 
models;  stresses  as  experimentally  determined  by  photo-elastic 
tubular  joint  models,  and  by  finite-element  analysis  of 
tubular  joint  models.  These  results  have  been  used  to  develop 
correlation  methods  for  estimating  elastic  hot-spot  stresses 
in  a  wide  variety  of  types  and  joints  and  parameters. 

Comparisons  between  measured  hot-spot  stresses  and  finite- 
element  analysis  involve  considerations  of  exactly  where  the 
hot-spot  stresses  are  located  and  which  type  of  data  is  more 
appropriate  for  use  with  the  AWS  X-curve.  Potvin  [2-16] 
discussed  this  aspect  using  the  illustration  in  Figure  2-21. 

Potvin  et  al  [2-16]  noted  that  there  is  a  ±20%  difference 
between  their  analytical  hot-spot  stresses  and  measured  hot¬ 
spot  stresses  and  they  attribute  this  to  the  difference  in  the 
analytical  and  experimental  hot-spot  locations.  This 
variation  in  hot-spot  location  is  a  result  of  two  factors:  the 
absence  of  the  weld  fillet  in  their  finite  element  model  and 
the  difficulties  involved  with  obtaining  consistent 
experimental  data.  Their  finite-element  model  predicts  that 
the  hot-spot  occurs  at  the  intersection  of  the  midsurface  of 
the  brace  and  chord.  The  inconsistency  in  the  experimental 
data  involves  the  actual  location,  with  respect  to  the  weld 
toe,  at  which  the  experimental  database  is  obtained.  The 
experimental  uncertainty  is  compounded  by  the  size  and  number 
of  strain  gauges  used  (i.e.,  is  the  gauge  small  enough  to  be 
placed  close  to  the  toe  of  the  weld  and/or  are  there  a 
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sufficient  number  of  gauges  so  that  an  accurate  extrapolation 
of  stresses  to  the  toe  of  the  weld  is  possible  as  indicated  in 
Figure  2-21) . 

The  test-data  available  and  correlation  methods,  in  addition 
to  being  restricted  to  the  elastic  regime,  are  also  restricted 
in  the  sense  that  they  are  based  on  "isolated"  joints  (i.e., 
joints  in  a  chord  where  any  other  joint  is  sufficiently  far 
away  so  that  the  two  joints  are  not  influenced  by  each  other) . 
In  offshore  platforms,  quite  often  two  or  more  joints  are 
placed  at  the  same  axial  location  on  a  chord.  There  are  no 
available  data  on  such  joints.  There  are  "rule  of  thumb" 
approaches  for  joint  spacing  such  that  data  on  isolated  joints 
would  be  reasonably  applicable. 

Correlation  methods  for  joint  hot-spot  stresses  based  on  joint 
configuration  parameters  have  been  developed  by  a  large  number 
of  researchers.  These  correlation  methods  include: 

1.  Marshall's  [2-17]  presentation  of  the  Kellogg  [2-18] 
equation  for  hot-spot  stresses  in  the  chord  and  his 
equations  for  hot-spot  stresses  in  the  brace; 

2.  Bijlaard's  [2-19]  method  for  hot-spot  stresses  in 
the  chord; 

3.  Kuang  et  al  [2-20]  method  for  hot-spot  stresses; 

4.  Wordsworth  [2-21]  for  chords; 

5.  Gibstein  [2-22]  for  chords. 

Other  authors  [2-23  through  2-26]  have  developed  correlation 
methods.  These  correlation  methods  are  developed  for  in-plane 
joints  as  shown  by  Kuang  in  Figures  2-14  and  2-15.  Gulati 
[2-7]  reviewed  several  methods  mentioned  above  and  compared 
them  to  finite  element  analyses,  the  results  of  which  are 
shown  in  the  same  figures. 

The  interaction  of  stress  fields  of  neighboring  connections  is 
accounted  for  by  direct  super  position  of  stresses.  However, 
results  have  been  proven  suspect  in  magnitude  and  location  of 
hot-spot  stress.  Additionally,  local  interaction  formulations 
are  based  on  parametric  formulations  and  rules  of  thumb  and 
these  limitations  will  be  discussed  in  Section  2.3. 
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Figure  2-21.  Illustration  of  Stress  Concentration  Factors  and 
Hot-Spot  Stress  for  Offshore  Structures 


2.2.2 


Proposed  Approach  for  Faticue  Analysis  of  Ship 
Structural  Details 


While  there  is  no  approach  specifically  written  into  codes  or 
rules  for  ship  structures  (other  than  empirical  stress 
allowables) ,  there  have  been  approaches  proposed  to  account 
for  fatigue  life  in  ship  structural  details  [2-3,  2-27  through 
2-29].  The  most  notable  approach  is  that  developed  by  Munse  et 
al  [2-3]  and  presented  in  SSC  318.  This  approach  is  based  on 
calculating  a  "design  allowable"  stress  range,  Srd  for  fatigue. 
This  stress  range  is  the  maximum  peak-to-trough  stress  range 
expected  at  the  point  in  guestion  once  under  the  most  severe 
sea  state  and  during  the  entire  life  of  the  structure. 
Comparing  that  stress  range  to  the  allowable  stress  for  other 
failure  modes  indicates  the  controlling  mode  of  failure.  In 
any  case,  the  maximum  stress  computed  from  the  fatigue  design 
stress  range,  Srd,  must,  generally  speaking,  be  less  than  the 
nominal  permissible  stress  permitted  once  by  the  basic  design 
rules.  According  to  the  Munse  approach,  the  design  stress 
range,  Srd,  is  found  using  the  following  equation: 

Srd  -  Sn  •  Rf  •  e 

where  S„  =  mean  value  of  the  constant  amplitude 

stress  range  at  the  design  life,  Nd 

Rf  =  reliability  factor 

e  =  random  load  factor. 

The  mean  value  of  the  constant  amplitude  stress  range,  Sn,  is 
found  by  entering  the  S-N  curve  of  the  structural  detail  of 
interest  at  the  number  of  cycles  expected  in  the  design  life, 
Nd .  The  probabilistic  nature  of  the  design  method  is 
introduced  by  the  other  two  factors  in  the  equation. 

The  reliability  factor,  R-,  is  meant  to  account  for 
uncertainties  in  the  fatigue  data,  workmanship,  fabrication, 
use  of  the  equivalent  stress  range  concept,  errors  in  the 
prediction  of  load  history  and  errors  in  the  associated  stress 
analysis.  The  factor  comes  from  the  assumption  that  fatigue 
life  is  a  random  variable  with  a  Weibull  distribution,  and  the 
use  of  a  relationship  for  the  probability  of  survival  through 
N  loading  cycles.  The  effect  of  the  reliability  factor  is  to 
reduce  a  mean  constant  stress  range  to  an  equivalent  stress 
range  which  corresponds  to  a  designated  probability  of 
survival  greater  than  the  50%  level  of  the  mean  stress  range. 
The  random  load  factor,  c,  is  introduced  in  the  design 


procedure  to  make  possible  the  use  of  existing  constant-cycle 
fatigue  data  in  designing  for  variable  loading  service  condi¬ 
tions. 

The  work  of  Munse  et  al  represents  a  significant  step  forward 
in  the  design  of  ship  structures.  It  presents  the  first  truly 
probabilistic  approach  to  fatigue  design  for  ship  structural 
details.  This  approach  integrates  well  with  the  ship  design 
process  where  a  nominal  level  of  stress  is  developed  and  no 
detail  stress  calculations  are  required.  Multiaxial  stress 
distributions  are  not  accounted  for  within  the  detail, 
however,  reliability  factors  are  presented  which  are  to 
account  for  these  effects  such  as  scatter  about  the  mean  line 
of  the  S-N  curve.  When  detailed  stress  analysis  is  warranted, 
other  approaches  are  required  for  correlation  to  the  basic  S-N 
curve  for  the  materials. 

2.2.3  Proposed  Fracture  Mechanics  Approaches  for  Marine 
Structures 


Thayamballi  et  al  [2-30]  proposed  a  fracture  mechanics 
approach  for  ship  structures  based  on  the  Paris  equation 
presented  above.  This  approach  takes  into  account  various 
factors  influencing  fatigue  response  and  reliability 
considerations.  There  are  approaches  presented  to  account  for 
multiaxial  loading  sources  and  complex  stress  fields. 

However,  the  authors  emphasize  Mode  I  stress  intensity  factors 
being  typical  of  ship  structural  loading  and  fatigue  crack 
growth . 

Chen  used  this  approach  for  the  fatigue  crack  growth  analysis 
of  the  SL-7  hatch  corner  cracking  [2-6].  An  integration  of 
the  Paris  equation  and  use  of  the  equivalent  stress  concept 
yield  the  cycles  to  failure.  As  described  earlier, 

Thayamballi  found  that  initial  crack  length  is  proven  a  major 
consideration  and  source  of  uncertainty. 

Fracture  mechanics  analysis  for  offshore  structures  has  been 
used  as  a  fitness  for  purpose  tool  for  some  time  and  has  been 
proposed  for  fatigue  life  analysis  by  Rhee  [2-31]  and  Haung 
[2-32]  among  others.  Again,  the  Paris  equation  for  linear 
elastic  fracture  mechanics  is  utilized.  The  important  steps 
associated  with  the  procedure  are  summarized  as  follows: 

1.  calculate  storm  member  forces  through  frame  analysis 
for  a  jacket  under  a  given  environment; 
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2.  determine  the  crack  location  for  fatigue  crack 
growth  simulation  analyses  based  on  the  stress 
analysis  results  and  defect  distribution  status; 

3.  calculate  the  stress  intensity  factors  for  the 
fatigue  crack  growth  simulation  analyses,  and  other 
fracture  parameters  required  for  crack  instability 
analyses.  Calibrate  the  material  fatigue  and 
fracture  properties  required  consistently  with  the 
fracture  parameters; 

4.  perform  fatigue  crack  growth  simulation  analyses; 

5.  perform  crack  instability  analyses  to  determine  the 
critical  crack  size; 

6.  correlate  the  results  of  the  fatigue  crack  growth 
and  crack  instability  analyses  to  determine  the 
component  fatigue  life. 

2.3  IMPETUS  FOR  A  MULTIAXIAL  FATIGUE  APPROACH 

Over  the  j.  nst  decade,  the  application  of  finite  element 
analysis  to  structural  details  has  proliferated,  especially 
with  the  advent  of  personal  computers  with  memory  capacity  for 
finite  element  computer  programs.  Increasing  use  of  this  tool 
has  made  it  possible  to  analyze  the  state  of  stress  in  a 
structural  detail  to  a  level  of  detail  beyond  that  required  by 
existing  nominal  stress  fatigue  design  procedures.  This  tool 
also  makes  it  quite  possible  to  predict  combined  stresses  from 
multiple  loads  and  from  geometric  concentrations.  The 
designer  now  has  a  choice:  1)  to  use  techniques  that  he  has  to 
predict  nominal  stresses  in  joint  boundaries  and  accept 
limitations  on  conservatism  associated  with  the  correlation 
equations  or  tesc  data  from  similar  welded  components;  or  2) 
cross  the  nominal  stress  boundary  and  predict  the  state  of 
stress  within  the  detail.  More  structural  designers  for  both 
ships  (ABS  [2-1]  and  Columbia  Research  Corporation  (CRC)  [2- 
34])  and  offshore  platforms  (Gulati  [2-7]  and  Rhee  [2-31])  are 
choosing  the  latter  option.  However,  the  new  knowledge  has 
also  raised  additional  questions. 

Existing  S-N  curves  and  crack  growth  data  is  typically 
presented  from  uniaxial  tests.  Finite  element  techniques 
estimate  principal  stresses  in  multiple  directions  resulting 
from  various  load  and  geometric  effects.  The  designer  now 
requires  additional  knowledge  on  how  to  characterize  the  state 
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of  stress  to  correlate  with  uniaxial  fatigue  life  predictions 
and  improve  design  reliability. 

This  situation  is  evident  as  a  result  of  a  finite  element 
study  conducted  by  Gulati  [2-7]  for  offshore  structures.  As 
shown  in  Figure  2-18,  the  principal  stress  in  the  neighborhood 
of  the  concentration  point  at  the  chord-brace  junction  of  a  T- 
joint  subjected  to  combined  loading  is  not  orthogonal  to  the 
weld  lineal  direction  assumed  by  the  current  fatigue  life 
prediction  methods.  Furthermore,  the  direction  of  the 
principal  stress  changes  with  the  changing  magnitudes  of  the 
constituent  loadings.  The  effect  of  the  non-perpendicularity 
to  the  weld  lineal  direction  of  the  principal  stress  (at  the 
"hot  spot")  on  the  crack  trajectory  and  fatigue  life  of  an  as 
welded  joint  is  not  known. 

Because  designers  are  asking  increasingly  difficult  guestions 
on  multiaxial  fatigue  response  in  real  structures,  it  became 
necessary  to  review  the  general  multiaxial  fatigue  research 
and  those  approaches  that  have  been  applied  to  welded 
structural  details  and  evaluate  to  their  applicability  to 
marine  structures. 
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3.0  REVIEW  OF  LITERATURE  ON  MULTIAXIAL  FATIGUE  RESEARCH 


Interest  in  problems  of  metal  fatigue  and  fracture  in 
engineering  has  a  long  history  dating  back  well  into  the 
last  century.  There  was  certainly  an  early  awareness  of  the 
problem  of  dealing  with  the  effect  of  combined  loads,  but 
the  major  thrust  of  early  fatigue  development  was  in  the 
field  of  uniaxial  fatigue  testing  and  the  development  of 
uniaxial  data  for  material  fatigue  characteristics  in  the 
form  of  the  well-known  S-N  diagrams.  Interest  was  centered 
in  the  area  of  long  material  lives,  that  is,  High  Cycle 
Fatigue  (HCF) .  In  this  regime  the  nominal  stresses  in  the 
material  are  in  the  elastic  range,  only  a  fraction  of  the 
yield  stress  for  the  material,  and  plasticity  is  confined  to 
the  region  in  the  immediate  vicinity  of  the  fatigue  crack 
tip.  The  discipline  of  fracture  mechanics  was  not  born 
until  after  World  War  I,  but  it  remained  until  World  War  II 
for  this  field  to  experience  really  active  development  [3-1 
through  3-8].  As  understanding  of  the  mechanisms  of  crack 
initiation  and  growth  improved,  interest  increased  in  the 
area  of  higher  loadings  and  shorter  lives  —  that  is,  the 
Low  Cycle  Fatigue  (LCF)  regime.  As  the  demands  of 
technology  made  clear  the  need  for  a  more  sophisticated 
understanding  of  the  entire  fatigue  problem,  testing 
machines  were  developed  which  had  the  capability  of  testing 
in  two  modes,  most  often  in  torsion  and  axial  push-pull, 
sometimes  in  torsion  and  bending,  and  less  frequently  in 
in-plane  linear  orthogonal  loads.  Multiaxial  fatigue 
research  has  been  active,  lively  and  growing  for  the  past 
twenty-five  years,  but  it  is  not  by  any  means  a  mature 
discipline.  There  are  a  multitude  of  proposed  multiaxial 
fatigue  criteria,  the  most  prominent  of  which  will  be 
discussed  in  subsequent  sections  of  this  chapter,  and  there 
are  codes  and  formulations  (e.g.,  American  Society  of 
Mechanical  Engineers  (ASME)  Pressure  Vessel  Code)  which  are 
in  place  and  may  be  used  on  a  provisional  basis,  but  there 
does  not  exist  a  single  unified  and  validated  criterion, 
accepted  by  the  research  community,  which  can  be  used  for 
engineering  design  applications. 

In  the  following  sections,  the  general  multiaxial  fatigue 
research  and  candidate  approaches  to  marine  structure  design 
and  evaluation  will  be  reviewed. 

3.1  GENERAL  MULTIAXIAL  FATIGUE  RESEARCH 

3.1.1  Fatigue  Life  Estimates  for  Crack  Initiation 

The  first  general  category  of  multiaxial  fatigue  research  is 
for  those  approaches  based  on  the  cumulative  damage  approach 
for  estimating  fatigue  life.  An  equivalent  stress  is  used 
to  correlate  multiaxial  stresses  or  strains  to  uniaxial  test 
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data.  Each  approach  is  applicable  to  a  specific  type  of 
fatigue  response.  Generally,  shear  stress  based  criteria 
are  applicable  to  fatigue  crack  initiation  in  the  high  cycle 
regime.  Principal  stress  criteria  are  applicable  to  fatigue 
life  approaches  that  include  crack  growth.  Combinations  of 
approaches  have  been  developed.  Approaches  have  also  been 
developed  for  estimating  fatigue  life  under  low  cycle 
conditions  where  plastic  work  dominates.  Traditional 
interest  had  been  in  long-lived  behavior  in  which  nominal 
stresses  in  the  material  remained  in  the  elastic  range  well 
below  the  tensile  yield  stress.  Figure  3-1  taken  from 
Reference  [3-8]  illustrates  schematically  the  partition 
between  low  cycle  fatigue  and  high  cycle  fatigue  regimes. 

The  differentiation  between  high  cycle  and  low  cycle  is  im¬ 
portant  in  further  categorizing  multiaxial  fatigue 
approaches . 

In  the  following  sections  the  general  multiaxial  fatigue 
research  developed  for  applications  in  other  industries  such 
as  pressure  vessels,  aircraft,  moving  vehicles  and  heavy 
equipment  is  reviewed.  Much  of  this  research  is  based  on 
extension  of  existing  uniaxial  research  with  newer 
approaches  based  on  plastic  work  estimates. 

3. 1.1.1  Stress  Based  Criteria 

Early  research  in  multiaxial  fatigue  was,  quite  naturally, 
an  extension  of  uniaxial  work  and  dealt  with  combined 
stresses  in  the  HCF  regime.  The  work  of  Mason  [3-9]  and 
Mason  and  Delaney  [3-10]  is  representative. 

Two  well-known  and  widely  used  yield  criteria,  those 
attributed  to  Tresca  and  to  Von  Mises  were  adapted  early  to 
cyclic  stresses. 

Tresca:  Tm„  =  (aj-a3)/2  =  constant 

Von  Mises:  roetlh.dr.,  =  1/3  j\ai-a#)*+  (cr2-a3)2 

1  1/2 

+  (a3-ffi)z  J  =  constant 

where  a3  >  az  >  a3  are  the  principal  stress  and  the  value  of 
the  constant  is  specified  for  a  given  fatigue  life. 

A  significant  contribution  was  made  by  Gough  and  colleagues 
in  the  mid-30's  through  the  early-50's  [3-11  through  3-14]. 
Gough  found  that  the  arc  of  an  ellipse  provided  a  good  fit 
to  his  experimental  data  from  in-phase  bending  and  torsion 
tests  of  brittle  materials  and  ductile  materials.  Guest 
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[3-15,  3-16]  noted  that  Gough's  empirical  formulations  could 
be  reduced  to  a  single  expression: 

0i-02  ffi+Oa 

-  +  k  -  =  constant  (3-1) 

2  2 

where  ax  >  a2  >  a3  are  principal  stresses  and  k  is  a  measure 
of  the  relative  ductility  or  brittleness  of  the  material. 
Guest's  relationship  is  a  statement  that  constant  fatigue 
lives  are  a  function  of  the  maximum  shear  stress,  (a,-a2)/2, 
modified  by  a  fraction,  k,  of  the  normal  stress  on  the  plane 
of  maximum  shear  stress,  (a!+a3)/2.  Variations  of  this 
theme  in  terms  of  both  stress  and  strain  will  be  seen  to 
occur  repeatedly  in  the  following  years.  Tipton  and  Nelson 
[3-17]  have  expressed  the  Guest  criteria  for  combined  in- 
phase  bending  and  torsion  as: 


where  a  =  amplitude  of  cyclic  bending  stress 

r  =  amplitude  of  cyclic  torsion  shear  stress 

bn  =  fatigue  limit  for  notched  specimens  in 

bending  (determined  experimentally) 

tn  =  fatigue  limit  for  notched  specimens  in 
torsion  (determined  experimentally) 

Figure  3-2,  taken  from  Reference  [3-17],  shows  an  excellent 
correlation  of  data  taken  from  multiaxial  notched  specimen 
tests  in  torsion  and  bending  with  the  torsion-bending 
fatigue  limit  predicted  by  (3-2)  above.  Also  shown  in 
Figure  3-2  are  fatigue  life  predictions  using  the  well-known 
Von  Mises  distortion  energy  yield  criteria  applied  to 
constant-amplitude  bending  and  torsional  stresses.  Stresses 
in  the  notch  of  the  specimen  have  been  calculated  in  two 
ways:  (1)  using  elastic  stress  concentration  factors,  K*; 

and  (2)  experimentally  determined  stress  concentration 
factors,  Kf.  In  this  case,  using  the  elastic  stress  concen¬ 
tration  factors  gives  a  very  poor  correlation.  Other  cases 
displayed  in  [3-17]  using  K*  are  in  better  agreement.  K, 
predictions  are  generally  better,  but  in  most  cases 
unconservative.  Further  confirmation  of  the  form  of  the 
Guest  Law  was  provided  by  Findley  [3-18].  McDiarmid  [3-20] 
also  proposed  a  criterion  for  fatigue  failure  based  on  the 
critical  range  of  shear  stress  modified  by  the  normal  stress 
acting  on  the  plane  of  maximum  shear  stress.  Sines  [3-19] 
proposed  a  tensorial  method  of  combining  biaxial  alternating 
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DIFFERENT  R  VALUES 
USEO  IN  TESTS 


Figure  3-2.  Comparison  of  Predicted  Combined  Bending-Torsion  Fatigue 
Limit  with  Test  Data  for  a  Circumferentially  Notched 
Round-Bar  Specimen  [3-17] 
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(in-phase)  stresses  with  a  static  mean  stress.  The  method 
is  general  in  nature  and  could  be  extended  to  triaxial 
stress  systems  with  a  superimposed  static  stress. 

Langer  [3-21]  has  discussed  the  application  of  Tresca's 
criterion  to  the  design  of  pressure  vessels  involving 
multiaxial  fatigue.  Tresca's  criterion  was  favored  over  Von 
Mises'  "criterion"  because  of  its  simplicity,  conservatism 
and  adaptability  to  more  situations.  This  particular 
criterion  was  adopted  as  the  basis  for  evaluating  multiaxial 
fatigue  in  the  ASME  Boiler  and  Pressure  Vessel  Code 
procedures  [3-22].  The  approach  is  as  follows: 

S,1t  =  the  value  of  (1/2)  (S1D-S3D)  maximized 
with  respect  to  the  least  and  greatest 
principal  stress  values  of  the  stress 
cycle  where  SID  and  S3D  are  the  maximum 
and  minimum  stresses  respectively. 

The  parameter  S,1t  is  called  "stress-intensity"  which  is  used 
as  the  parameter  correlating  fatigue  life  (to  crack 
initiation) .  It  is  the  maximum  range  (over  the  complete 
"cycle"  of  loading)  of  shear  stress  acting  on  a  particular 
pair  of  planes  and  in  a  particular  direction;  its  absolute 
value  is  important  and  not  the  planes  and  directions  along 
which  it  acts.  This  is  a  stress  based  criterion.  This 
approach  has  been  refuted  under  out-of-phase  multiaxial 
loading  because  of  the  path  independence  assumption. 

Reviews  of  the  ASME  Boiler  and  Pressure  Vessel  Code  exist 
[3-23,  3-24]  indicating  mixed  results  for  out-of-phase 
loading.  In  the  case  of  biaxial  in-phase  constant  amplitude 
loadings,  the  principal  stresses  will  act  in  a  constant 
direction,  but  with  oscillating  amplitudes.  The  out-of- 
phase  monofrequency  constant  amplitude  problem  is  much  more 
realistic  for  many  rotating  machinery  and  pressure  vessel 
applications.  Note  that  most  of  the  following  equivalent 
stress  approaches  are  presented  in  a  format  in  which 
sinusoidal  loading  functions  are  known,  for  example,  in 
applications  such  as  rotating  machinery. 

More  recently,  ASME  adopted  the  following  formulation  in  the 
1974  ASME  Boiler  and  Pressure  Vessel  Code  [3-25]  for  an 
equivalent  stress,  SALT  (S.1t)  which  for  fully-reversed 
bending  and  torsion  can  be  expressed  as: 
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SALT 


B 

J~~Z 


1  +  K*+  /  1  +  2K*COs20  +  K4 


1/2 


(3-3) 


where  B 
T 
K 
<P 


amplitude  of  cyclic  bending  stress 
amplitude  of  cyclic  torsion  stress 
2T/B 

phase  angle  between  bending  and  torsion. 


In  a  modification  of  the  above  method,  an  equivalent  stress 
parameter,  SEQA  (S.qulw.) ,  was  presented  as  an  ASME  Code  case 
[3-26].  For  fully  reversed  out-of-phase  bending  and  torsion 
the  parameter  is  given  by: 


B 

SEQA  =  - 

/"2 


3  /  3  9 

l  +  -Kz  +  /l  +  -  KZCOS20  +  —  K4 

4  2  16 


1/2 
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It  can  be  shown  that  the  SEQA  formulation  reduces  to  the  Von 
Mises  criteria  for  in-phase  loadings.  Lee  [3-27]  developed 
a  criterion  in  the  spirit  of  Gough's  ellipse  quadrant 
applicable  to  constant  amplitude,  monofrequency,  fully 
reversed,  out-of-phase  torsion  and  bending  loads. 


where 


a 

fi 


1/a 


SLEE  =  b  I  (B/b)  +  (T/t) 

b 


(3-5) 


bending  fatigue  strength  for  a  given 
life  N 

torsional  fatigue  strength  for  the  same 
life  N 

variable  power  dependent  on  phase 
difference,  <p  and  material 

2(1  +  Bsin0) 

material  constant 


Lee,  in  Reference  3-27,  reports  a  more  favorable  correlation 
with  test  data  using  Equation  3-5  than  with  Equations  3-3  or 
3-4. 


The  ASME  Codes  are  widely  considered  to  be  quite 
conservative.  Tipton  and  Nelson  [3-17]  confirm  this  but 
note  that  out-of-phase  loadings  are  more  damaging  than 
in-phase  loadings,  a  feature  not  well  predicted  by  the  SEQA 
criterion.  Lee  [3-27]  also  found  out-of-phase  loadings  to 
be  more  damaging.  For  90°  out-of-phase  loadings  and  for 
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certain  ranges  of  the  parameter  K,  Lee  found  that  both  SALT 
and  SEQA  underestimated  the  danger  of  fatigue  failure. 

The  central  theme  in  the  development  of  stress-based 
criteria  has  been  the  reduction  of  a  complex  stress  state  to 
a  single  scalar  parameter  which  can  be  related  to  the  cyclic 
material  properties  and  lives  determined  from  uniaxial 
fatigue  tests  and  used  to  predict  life  under  the  more 
complex  loading  conditions.  While  cyclic  applications  of 
Tresca  and  Von  Mises  yield  criteria  are  the  most  prevalent, 
there  are  a  number  of  others,  each  with  their  proponents. 
Stress-based  criteria  find  their  greatest  success  in  the  HCF 
(long-lives)  regime  where  nominal  stresses  remain  in  the 
elastic  range  and  plastic  straining  effects  remain  very 
localized. 

Brown  and  Miller  (3-28]  note  that  a  common  feature  of  all 
biaxial  HCF  criteria  is  that  they  contain  two  constants. 

For  example,  in  the  case  of  torsional  and  bending  loading 
systems  the  constants,  t  and  b,  determined  experimentally 
from  separate  tests,  are  present.  But  for  bending  tests  the 
value  of  b  depends  on  the  physical  size  of  the  specimen  and 
will  vary  from  one  test  situation  to  another.  Standard 
push-pull  fatigue  tests  do  not  exhibit  this  size  dependent 
behavior.  Brown  and  Miller's  comment  does  not  bode  well  for 
the  eventual  development  of  a  completely  general  criteria. 

3. 1.1.2  Strain  Based  Criteria 

Strain-based  life  prediction  methods  gained  momentum  in  the 
1960's  with  increased  interest  in  the  LCF  regime.  Again, 
the  use  of  an  equivalent  amplitude,  in  this  case  a  strain 
amplitude,  is  determined  from  measured  or  calculated  strain 
histories  as  an  argument  to  enter  an  e-N  curve  to  determine 
fatigue  crack  initiation  life.  On  one  hand,  strain-based 
methods  are  inherently  better  suited  for  applications  where 
loading  is  such  that  relatively  large  regions  of  plasticity 
are  involved.  On  the  other  hand,  this  requires  a  knowledge 
of  the  cyclic  plastic  behavior  of  the  material  and  component 
including  variation  in  Poisson's  Ratio,  u  ,  with  strain 
level.  Thus  while  strain-based  methods  may  be  more 
realistic,  strain-based  life  predictions  are  not  as 
tractable  for  engineering  purposes  as  stress-based  methods. 

Application  of  Tresca  and  Von  Mises  criteria,  written  in 
terms  of  strains  rather  than  stresses,  is  an  obvious 
starting  point.  The  Tresca  method  assumes  that  the  maximum 
shearing  strain  amplitude  under  multiaxial  cr  3itions  will 
correlate  fatigue  life  with  that  determined  lion  the  maximum 
shearing  strain  in  uniaxial  tests.  For  bending  and  torsion 
this  is  given  by  Tipton  and  Nelson  [3-17]  as: 


3-8 


(1+  )2  + 


=  greater  of 


(1+u) 


1/2 

(in-phase  loading) 


(90°  out-of -phase 
loadings) 


where 


u  = 


torsional  shear  strain  amplitude 
bending  strain  amplitude 
Poisson's  Ratio 


The  equivalent  strain  amplitude  is  given  by: 


e.  = 


.y(l+u) 


In  the  case  of  a  standard  SAE  uniaxial  bending-torsion 
fatigue  specimen,  the  fatigue  strain  life  for  a  normalized 
1045  steel  is  given  by  Tipton  and  Nelson  as: 


0.00481(2N,) 


+  0.182(2^) 


where  e,  =  strain  amplitude 

2N/  =  cyclic  reversals  to  failure  (two 

reversals  per  cycle) 

Results  are  shown  in  Figures  3-3  through  3-5. 

Note  that  the  coefficients  and  exponents  in  the  above 
equation,  known  as  the  "strain-life"  equation,  are  material- 
specific  and  at  this  time  are  available  for  only  a  limited 
number  of  materials. 

Tipton  and  Nelson  have  also  evaluated  a  Von  Mises  strain- 
based  approach  under  assumptions  of  eyy  =  -ue„,  Figure  3-4, 
and  u  =  u ,  Figure  3-5,  where  u  is  a  variable  Poisson's  Ratio 
calculated  by  Gonyea's  method  [3-29].  The  effective  strain 
used  for  entering  the  uniaxial  e-N  curve  is  given  by: 


=  [  <«,- 

2(l+u)  L 


ez)z  +  (ea-c,)8  +  (e 


J 


where  €,,  e2,  e3  are  principal  strain  amplitudes  and  u  is 
Poisson's  Ratio  =  .29  or  .50. 


Tipton  and  Nelson  comment  on  several  problems  that  are 
characteristic  of  fatigue  testing.  First,  there  is  not  a 
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TRESCA  EQUIVALENT  STRAIN 


Figure  3-3.  Comparison  of  Life  Predictions  Based  on  the  Tresca  Equivalent 
Strain  Criterion  with  SAE  Test  Data  [3-17] 


IN-PHASE  VON  MISES  STRAIN;  EVV  =  — I'KXX 


Figure  3-4.  Comparison  of  Life  Predictions  Based  on  the  Von  Mises  Equivalent 

Strain  Criterion  with  In-Phase  SAE  Test  Data,  Taking  Circumferential 
Notch  Strain  as  Poisson's  Ratio  times  Notch  Bending  Strain  [3-17] 
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IN-PHASE  VON  M1SES  STRAIN;  u  =  T? 


Figure  3-5.  Comparison  of  Life  Predictions  Based  on  the  von  Mises  Equivalent 
Strain  Criterion  with  In-Phase  SAE  Test  Data,  Using  a  Variable 
Poisson's  Ratio  [3-17] 


Figure  3-6.  Actual  versus  Predicted  Lives: 

Maximum  Principal  Strain  Theory  [3-8] 
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commonly  accepted  standard  for  the  size  of  a  crack  that 
constitutes  initiation.  A  crack  detectable  by  visual 
inspection  is  one  definition.  In  other  cases,  optical 
microscopes  are  used  or  ultrasonic  acoustic  transducers  may 
be  employed.  0.1mm  and  1mm  crack  lengths  are  other 
definitions  which  are  used.  "Life"  may  be  life  to  crack 
initiation,  or  life  to  failure  which  may  be  defined  as  a 
percentage  drop  in  load-carrying  capability.  Ten  percent, 
twenty  percent  and  one  hundred  percent  drop-offs  are  often 
used.  Tipton  and  Nelson  use  a  fifty  percent  drop-off  in 
Reference  [3-17]. 

A  further  difficulty  in  using  uniaxial  data  is  that, 
depending  on  the  loading  modes,  their  relative  amplitudes 
and  the  multiaxial  strain  state,  the  type  of  cracking 
experienced  by  a  multiaxial  specimen  may  be  quite  different 
than  that  experienced  by  the  uniaxial  specimen  under 
"equivalent"  conditions. 

It  may  well  be  that  application  of  equivalent  uniaxial  data 
to  predict  multiaxial  fatigue  phenomena  will  be,  at  best,  an 
interim  procedure  which  will  be  useful  until  more  general 
multiaxial  criteria  are  developed  together  with  specifically 
applicable  supporting  material  data. 

The  use  of  a  constant  value  for  Poisson's  Ratio  is  also  a 
problem  which  must  be  dealt  with  when  large  plasticity  or 
high  temperatures  are  present.  Gonyea  [3-29]  has  proposed 
an  iterative  method  for  determining  a  variable  value  for 
Poisson's  Ratio. 

Brown  and  Miller  [3-30  through  3-32]  have  presented  an 
important  strain-based  multiaxial  fatigue  life  criterion 
which  is  a  natural  extension  of  Guest's  earlier  work.  The 
approach  is  one  of  several  known  as  "critical  plane" 
methods.  It  assumes  that  initiation  is  dominated  by  the 
amplitude  of  the  maximum  shear  strain,  Y  and  that 
propagation  is  strongly  influenced  by  the  strain,  e„,  normal 
to  the  plane  on  which  "Ym.,  is  acting,  a  plane  which  they 
called  the  T-plane.  The  general  form  of  the  criterion  they 
recommended  formed  a  family  of  ellipses  with  life  as  a 
parameter. 


The  specific  application  is  to  strain-controlled  in-phase 
combined  tension  and  torsion  tests,  with  experimentally 
determined  values  of  g,  h  and  j  which  were  found  to  vary 
with  life.  Figure  3-6  shows  these  results. 
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Fash  et  al  [3-8]  conducted  experiments  on  a  thin-walled  tube 
specimen  in  tension-torsion  and  on  a  solid  notched  shaft  in 
bending-torsion  and  compared  results  with  predictions  based 
on  five  different  theories.  They  concluded  that: 


1.  maximum  principal  multiaxial  strain  correlates 
life  with  maximum  principal  uniaxial  strain 
(Figure  3-6)  ; 

2.  the  effective  strain  (Von  Mises)  may  be  used  to 
obtain  life  from  uniaxial  data  (Figure  3-7) ; 


3.  maximum  multiaxial  shear  strain  may  be  correlated 
by  converting  uniaxial  data  into  shear  strains 
(Tresca  method)  (Figure  3-8) ; 

4.  a  theory  proposed  by  Lohr  and  Ellison  [3-33]  is 
based  on  the  maximum  shear  strain, Y  ,  which  occurs 
on  a  "critical  plane"  which  in  the  case  of 
tension-torsion  loading  of  thin-walled  tubes 
occurs  on  planes  intersection  the  surface  at  45°. 
The  maximum  shear  strain  range,  Y  ,  is  modified 
by  a  fraction  of  the  strain  normal  to  the  critical 
planes,  en  (Figure  3-9) ; 

5.  a  theory  proposed  by  Kandil,  Brown  and  Miller  [3- 
34],  essentially  an  extension  of  the  earlier  Brown 
and  Miller  work  [3-31,  3-32],  is  similar  in 
appearance  to  that  of  Lohr-Ellison  but  the  methods 
differ  in  the  definition  of  the  shear  strain  and 
normal  strain  components  (Figure  3-10) .  Lohr  and 
Ellison  and  Kandil,  Brown  and  Miller  approaches 
yield  essentially  equivalent  correlations  of  thin- 
walled  tube  specimen  fatigue  life. 

» 

The  relationship  between  fatigue  life  and  applied  strain 
amplitude  for  uniaxial  tests  is  given  by  the  strain-life 
equation: 


AC  Of  ' 

—  =  -  (2N/)b  +  6/  (2N,)C 

2  E 


where  A  c 

°f' 

«/' 

b 


axial  strain  range 
fatigue  strength  coefficient 
fatigue  ductility  coefficient 
fatigue  strength  exponent 
fatigue  ductility  exponent 


3-13 


100,05 


Figure  3-7.  Actual  versus  Predicted  Lives: 

Effective  Strain  Theory  r 3—8 ] 


Figure  3-8.  Actual  versus  Predicted  Lives: 

Maximum  Shear  Strain  Theory  [3-8] 


Figure  3-9.  Actual  versus  Predicted  Lives: 

Lohr  and  Ellison  Theory  [3-8] 


Figure  3-10.  Actual  versus  Predicted  Lives: 

Kandil,  Brown  and  Miller  Theory  [3-8] 
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E 

Young's  Modulus 

N /  = 

cycles  to  failure 

2N f  = 

reversals  to  failure 

Multiaxial  strains  as  determined  in  the  above  theories  are 
equated  to  the  strain  life  equation  above  with  coefficients 
and  exponents  appropriate  to  the  material.  For  the  work  of 
Fash  et  al,  failure  was  defined  as  a  10%  load  drop-off  for 
the  thin-walled  tube  tests  and  the  occurrence  of  a  1  mm 
crack  for  the  notched  shaft  specimen  tests. 

In  the  uniaxial  test  the  maximum  principal  strain  and 
amplitude  is  equal  to  the  applied  axial  strain  amplitude. 
Life  estimates  based  on  the  range  of  maximum  principal 
strain  are  based  on  a  similar  form  of  the  strain-life 
equation: 


Ac  Of ' 


(2N,)b  +  6/  (2N,)C 


2  E 

where  Ac  '  is  the  range  of  maximum  strain. 
Effective  (Von  Mises)  strains  are  given  by: 


2 

e  =  — 
3 


(er-e2)2  +  (e2-e3)2  +  (€,-€3)  = 


then 


A  €  Of 

—  =  —  (2Nf) b  +  ef  (2N,)C 
2  E 


For  uniaxial  data: 


Y 


1/2 


(u  =.5) 


»«  AC 

—  =  (1+u)  — 
2  2 


For  multiaxial  loadings,  the  maximum  shear  strain  from 
Mohr's  circle  is: 


Ym 


=  A  (€,-€3) 


(€,  >  C2  >  €3) 
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Substitution  of u=.3  and  u=.5  for  elastic  and  plastic  values 
of  Poisson's  Ratio  gives  the  maximum  shearing  strain 
(Tresca^  criterion; 


AY  ~x  Of  • 

-  =  1.30  -  (2N/)b  +  1.50  €/'(2N/)e 

2  E 

For  the  Lohr  and  Ellison  criterion: 


—  =  A  (Ci-e2) 

2 

A  (ei+e2) 

and  en  =  - 

2 

where  the  asterisk  is  used  to  designate  the  Lohr  and  Ellison 
parameters,  substituting  appropriate  values  of  Poisson's 
Ratio  gives  the  Lohr  and  Ellison  form  of  the  strain-life 
equation. 


Y* 

—  +  .4  €„  =  1.44  (2N/)b  +  1.60  ef  (2N f)c 
2 

The  Kandil,  Brown  and  Miller  criterion  differs  in  the 
definition  of  en  and  adds  a  term  modifying  the  maximum  shear 
strain  range, 

\  max. 


(€,-e3) 
6n  =  A - 

2 


Substituting  appropriate  values  of  Poisson's  Ratio  yields 
the  Kandil,  Brown  and  Miller  variant  of  the  strain-life 
equation. 


y  max  Of  f 

-  +  k€n  =  1.65  -  (2N,)b  +  1.75  € f  (2N;)C 

2  E 

A  comparison  of  actual  lives  versus  the  predicted  lives  by 
the  five  theories  is  shown  for  the  smooth  thin-walled  tube 
tests  in  Figures  3-6  through  3-10.  Correlation  for  all  five 
methods  falls  within  a  scatter  band  of  three  of  the  perfect 
correlation  line,  except  for  the  torsion  data  which  is 
attributed  by  the  authors  to  anistropy  measured  in  the  1045 
steel  specimens.  The  Lohr  and  Ellison  theory  appears  to 
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give  the  slightly  better  results  but  not  enough  to  justify 
any  strong  preferences. 

The  results  for  the  notched  specimen  tests  are  shown  in  the 
Fash  et  al  paper  but  are  not  included  here.  Correlations 
based  on  all  five  theories  were  rather  poor  and 
unconservative  showing  that  there  is  still  much  to 
understand  about  the  effect  of  stress  concentrations  in 
notches . 

At  this  point  it  appears  that  the  Kandil,  Brown  and  Miller 
approach  and  the  Lohr  and  Ellison  approach  are  emerging  at 
the  head  of  the  pack  in  multiaxial  fatigue  theories,  thus 
the  emphasis  placed  in  the  preceding  paragraphs.  Leese  [3- 
35]  in  a  recent  state-of-the-art  review,  shows  the 
difference  in  the  critical  planes  assumed  in  the  two 
theories  on  a  Mohr's  circle  representation  of  the 
three-dimensional  strain  state,  Figure  3-11.  Qualitatively, 
the  difference  between  the  two  theories  is  that  Lohr  and 
Ellison  reflect  tendencies  for  the  crack  to  be  driven 
through  the  specimen  whereas  the  Kandil,  Brown  and  Miller 
approach  reflects  the  strains  which  drive  surface  cracking. 

Socie,  Waill  and  Dittmer  [3-36]  conducted  standard  uniaxial 
tests  using  Inconel  718  specimens  and  biaxial 
strain-controlled  tests  under:  (a)  axial  loading  only;  (b) 
torsional  loading  only;  and  (c)  combined  axial  and  torsional 
loading.  Results  of  the  uniaxial  tests  are  given  in  Figure 
3-12  for  three  different  criteria  for  failure:  occurrence 
of  a  crack  0.1mm  long,  occurrence  of  a  crack  1.0mm  long,  and 
separation  of  the  specimen.  For  this  material,  strain  at 
the  .1mm  crack  level  was  15%  of  that  at  separation,  while 
strain  at  the  1.0mm  crack  was  90%  of  that  at  separation. 
Socie  et  al  also  examined  modifications  to  the  Lohr  and 
Ellison  parameter  and  to  the  Kandil,  Brown  and  Miller 
parameter  for  the  effect  of  mean  stress,  an  important 
c'  -jsideration  in  ship  fatigue.  Both  modified  theories  gave 
reasonably  good  correlations,  but  Lohr  and  Ellison  was 
slightly  the  better.  The  comparison  is  shown  in  Figure 
3-13.  The  constants  in  the  Lohr  and  Ellison  equation  may  be 
obtained  from  uniaxial  tests. 

Brown  and  Miller  [3-37]  have  identified  two  distinct 
cracking  systems  which  they  define  as  Case  A  and  Case  B, 
each  of  which  gives  rise  to  different  fatigue  lives.  The 
two  cases  are  depicted  schematically  in  Figure  3-14.  In 
Case  A,  strains  €3  and  €3  are  parallel  to  the  surface  (as  in 
the  case  of  combined  tension  and  torsion)  and  tend  to 
generate  shallow  surface  cracks.  Under  other  loadings  e3  is 
normal  to  the  surface  tending  to  drive  deep  cracks.  This  is 
indicated  in  Case  B.  Different  functions  were  proposed  for 
each  case. 
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Figure  3-11.  Difference  in  Critical  Planes  of  Brown  and  Miller 

and  of  Lohr  and  Ellison  for  Tension/Torsion  Loading 
on  Mohr's  Circles  of  Strain  (from  Ref.  [3-35]) 


3-19 


2/*V‘*Pn*Hdu,V  UI°JIS 


3-21 


Figure  3-13.  Correlation  of  Test  Data  Including  Mean 

Stress  Effects  (a)  Lohr  and  Ellison  Para¬ 
meter  and  (b)  Kandil,  Brown  and  Miller 


MMdPAt.  STRAINS  PLANES  Of  MAXIMUM  SMEAR  PLANE  ft  DIRECTION 

C  »C.»C.  STRAW  AND  MRCCTON  Of  Of  STAGE  S 

'  STAGE  1  CRACK  GROWTH  CRACK  GROWTH 


SURfACE  PLAKE 


Figure  3-14.  Stage  I  and  Stage  II  Crack  Growth  Systems 
Under  General  Multiaxial  Cyclic  Strains 
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Other  authors  also  have  made  contributions  in  the  area  of 
multiaxial  strain-based  methods.  Konter,  Janssen  and 
Husslage  [3-38]  have  presented  a  quadratic  method  for 
equivalent  strain  range  using  plastic  strain  components: 

£.q2  =  (16-96)  cn2  +  6  (Y  mx/2  ) 2 

Different  values  of  the  constant,  6,  apply  to  Cases  A  and  B. 
Zamrik  and  Frishmuth  [3-39]  proposed  use  of  maximum  total 
strain  to  correlate  out-of-phase  fully  reversed  bending- 
torsion  tests: 

«T  =  J  €i2  +  «22  +  e32 

Grubsic  and  Simburger  [3-24]  proposed  examining  all  planes 
in  a  body  and  selecting  the  most  unfavorable  combination  of 
mean  and  alternating  shear  stress  on  each  plane  as  the  basis 
for  a  life  criterion. 

3. 1.1. 3  Plastic  Work  Methods 

The  difficulties  associated  with  extending  stress-based  or 
strain-based  criteria  to  the  more  general  cases  of 
cumulative  damage  under  random  loadings  of  mixed  frequency 
and  phase  has  led  a  number  of  investigators  to  consider 
methods  in  which  the  material's  resistance  to  fatigue  under 
multiaxial  loading  conditions  can  be  related  to  the  work 
expended  in  each  stress-strain  cycle.  Early  foundations  for 
energy  approaches  were  laid  by  Feltner  and  Morrow  [3-40], 
Morrow  [3-41]  and  Halford  [3-42]  during  the  1960's.  In 
1977,  Leis  [3-43]  proposed  a  damage  parameter  utilizing 
strain  energy  per  cycle  which  was  intended  to  account  for 
multiaxiality  of  stress  and  stain,  mean  stress,  and  viscous 
deformation  response.  The  method  assumed  isothermal  stress- 
strain  cycles  of  constant  strain  amplitude  obtained  from 
uniaxial  tests.  Using  existing  published  data,  Leis 
obtained  an  excellent  data  collapse  for  uniaxial  data  using 
a  log  energy  versus  a  log  life  plot,  but  the  collapse  of 
biaxial  data  was  much  less  impressive.  Leis'  definition  of 
failure  was  the  termination  of  Stage  I  crack  growth.  In 
1981,  Garud  [3-44]  proposed  an  energy  method,  again  based  on 
evaluation  of  the  plastic  work  per  strain  cycle,  but  now 
introducing  a  modified  strain  hardening  method  into  the 
stress-strain  constitutive  relations  for  the  material. 

Garud  used  the  axial-torsional  data  of  Kanazawa,  Miller  and 
Brown  [3-45]  to  correlate  his  results.  Figure  3-15  shows 
the  uniaxial  stress-strain  curve  used  to  determine  the 
plastic  work  per  cycle.  Figure  3-16  shows  log  plastic 
energy  versus  log  life  for  the  various  tests  in  the  Kanazawa 
data.  Note  that  the  axial  test  data  falls  on  the  lower 
bound  and  torsional  data  falls  on  the  upper  bound  of  a 
factor-of-three  scatter  band  with  the  phase  angle  biaxial 
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Axial  Strain  Amplitude  (percent) 


Figure  3-15.  Uniaxial  Stress-Strain  Approximation  for  1  Percent 
CR-Mo-V  Steel  (segments  4,  5  and  6  extrapolated; 
segments  5  and  6  not  shown) 
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Plastic  WorV Cycle,  Wc  (MJ/m3) 


Figure  3-16.  Plastic  Work  Per  Cycle  (Calculated)  Versus 
Life  to  Failure  (Observed) 
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tests  generally  falling  in  between.  When  a  weighing  factor 
of  h  was  applied  to  the  plastic  work  done  by  the  shear 
straining,  the  correlation  improved  as  shown  in  Figure  3-17. 
The  effect  of  phase  angle  is  shown  in  Figure  3-18.  At  the 
higher  axial  strains  (above  the  point  C) ,  the  90°  out-of¬ 
phase  loading  absorbs  the  most  plastic  work,  confirming  the 
findings  of  other  investigators  that  this  condition  is  the 
most  damaging  for  this  type  of  loading.  Garud  states  that 
hydrostatic  stress  and  mean  stress  effects  may  be  included 
in  the  method;  however,  strain  rate  and  creep  effects  are 
excluded. 

Ellyin  and  his  colleagues • [3-46  through  3-49]  have  also 
worked  on  the  development  of  energy  based  multiaxial  fatigue 
criteria  throughout  the  1980's.  The  criterion  proposed  in 
Ellyin' s  recent  paper,  Reference  3-49,  is  based  on  the 
theory  that  the  damage  caused  as  a  result  of  cyclic  loading 
is  a  function  of  the  mechanical  work  input  to  the  material. 
The  total  work  per  cycle,  Aw‘,  is  partitioned  into  elastic, 
W*,  and  plastic,  AWP,  components  with  separate  formulations 
for  each.  The  elastic  and  plastic  components  are  depicted 
graphically  for  a  uniaxial  cyclic  loading  case  in  Figure 
3-19.  The  Ellyin  criterion  includes  the  effect  of 
hydrostatic  pressure  and  mean  stress  but  in  its  present  form 
is  limited  to  proportional  or  nearly  proportional  loadings 
in  the  biaxial  mode.  This  implies  that  the  principal 
stresses  are  nearly  constant  in  direction  through  the  cycle. 
Experimental  results  for  two  strain  ratios  are  reguired  to 
determine  the  constants  in  Ellyin 's  equation;  however,  these 
coul^  be  determined  from  relatively  simple  tests  such  as  a 
standard  uniaxial  push-pull  fatigue  test  and  a  simple 
torsion  fatigue  test.  Results  of  such  tests  are  generally 
available  for  a  wide  range  of  materials.  Ellyin  and  Golos 
have  applied  the  criterion  to  correlate  test  data  obtained 
from  thin-walled  tube  tests  with  alternating  internal 
pressure.  Results  for  four  different  strain  ratios  are 
shown  in  Figure  3-20. 

Plastic  strain  energy  methods  have  great  appeal  because  of 
their  intrinsic  flexibility.  However,  there  are  practical 
considerations.  They  are  basically  better  suited  for  the 
LCF  regime  where  high  straining  is  present.  For  the  HCF 
regime  where  straining  levels  are  much  lower,  the  amount  of 
plastic  work  is  low  and  much  more  difficult  to  quantify.  An 
assumption  present  in  plastic  work  methods  is  that,  after 
the  first  few  cycles  of  straining,  the  constitutive 
relationships  remain  constant.  For  materials  or  under 
loading  conditions  for  which  this  assumption  is  less 
realistic,  the  method  would  be  much  more  complex.  A  summary 
of  the  multiaxial  fatigue  life  approaches  is  presented  in 
Table  3-1. 
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Plastk:  WocVCycte,  We  (MJ/m3) 


Life  to  Failure.  N  (cycles) 


Figure  3-17.  Plastic  Work  Per  Cycle  (Calculated  with  a 
Weighing  Factor  on  Shear  Work)  versus  Life 
to  Failure  (Observed) 
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Figure  3-18.  In-Phase  vs.  Out-of-Phase  Straining  (Comparison) 
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Figure  3-19.  Elastic  and  Plastic  Strain  Energy  Densities  for 
a  Uniaxial  Cyclic  Loading  Case 
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Strain  Energy  Density.  AW1  (MJ/m3) 


Cycles  to  Failure.  N; 


Figure  3-20.  The  Predicted  Strain  Energy  Density  Wc  versus  the 
Number  of  Cycles  to  Failure,  N/,  for  Various  Strain 
Ratios 
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3.1.2  Fatigue  Crack  Propagation  Research 


General  research  on  fatigue  crack  propagation  using  fracture 
mechanics  has  been  reviewed  and  is  presented  here.  While 
few  approaches  have  been  proven  for  incorporation  into 
design  procedures,  they  do  provide  insight  for  future 
directions  in  crack  growth  studies. 

3. 1.2.1  Observations  of  Mixed  Mode  Crack  Growth 

Smith  and  Pascoe  [3-50]  show  that  in  HY100  steel,  cracks  can 
either  initiate  and  grow  in  shear  (Stage  I-Mode  II)  or 
change  to  grow  in  the  opening  mode  (Stage  I I -Mode  I)  under 
different  stress  status.  They  tested  thick  plate  specimens 
under  the  applied  strain  conditions  in  Table  3-2. 

Under  Mode  II  dominated  conditions,  crack  bifurcation  was 
observed  which  had  a  significant  effect  on  crack  growth 
rates  as  shown  in  Figure  3-21.  Under  these  conditions  the 
Paris  equation  is  invalid.  They  also  noted  that  crack  tip 
plasticity  had  a  large  effect  on  crack  growth  rates  for 
HY100 .  The  most  significant  finding  of  their  study  is  the 
stable  crack  growth  under  Mode  II  loading.  This  finding  is 
attributed  to  crack  tip  plasticity  of  this  ductile  material. 
Most  literature  indicates  a  lack  of  full  understanding  of 
crack  propagation  in  mixed  mode  stress/strain  fields. 

Brown  and  Miller  [3-50]  show  the  importance  of  crack  tip 
plasticity  in  controlling  crack  speed.  Crack  growth  rates 
are  always  faster  than  linear  elastic  fracture  mechanics 
(LEFM)  predictions  suggest.  Their  rudimentary  approach  to 
predicting  propagation  rates  requires  further  verification. 

The  work  of  Kitagawa  et  al  [3-51]  confirms  that  LEFM  is  only 
applicable  below  one-third  of  yield  stress  and  that  at 
higher  stress  levels  there  is  a  need  to  account  for  biaxial 
effects.  They  also  confirm  the  observation  that  long  cracks 
under  mixed  mode  loading  branch  to  follow  the  direction 
where  K„  is  zero  or,  in  other  words,  in  a  direction  normal 
to  the  principal  stress. 

In  contrast,  Gao  [3-52]  observed  stable  mixed  mode  growth  at 
stress  levels  close  to  threshold  in  four  materials.  Smith 
and  Pascoe  [3-50]  show  that  stable  Mode  II  growth  is 
possible  for  long  cracks  and  high  strain  amplitudes.  This 
shows  that  at  least  one  question  still  to  be  resolved  is  why 
shear  mode  cracks  should  change  to  Mode  I.  Again,  large- 
scale  yielding  (or  lack  of)  in  shear  planes  is  the  suspected 
cause . 

Mode  III  crack  growth  studies  were  conducted  by  Ritchie  et 
al  [3-53],  Hourlier  et  al  [3-54]  and  Pook  [3-55].  The 
researchers  found  it  difficult  to  maintain  Mode  III  growth 
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TABLE  3-2 


ANGLE  OF  CRACK  ON  SURFACE 

APPLIED  STRAIN  STATE 

£  =  A£2 
A£i 

-1  (shear) 

-0.5 

0  (plane  strain) 

+0.5  (uniaxial  strain) 

+1  (equibiaxial) 


AS  FUNCTION  OF  STRAIN  STATE 

CRACK  ANGLES  FROM 
PRINCIPAL  AXIS 

45°  and  135° 

55 *  and  125“ 

90* 

Between  45°  and  135° 
Any  angle 
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Figure  3-21.  Crack  Propagation  Plotted  as  Crack  Length  versus  Crack 
Growth  Rate  to  Compare  Mode  I,  Mixed-Mode,  and  Mode  II 
Growth 
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for  long  cracks  except  at  the  highest  stress  levels  where 
presumably  large  crack-tip  plastic  zones  permitted  Mode  III 
extension.  The  fracture  surface  morphology,  frequently 
observed  in  torsional  loading  situations,  indicates  a 
preference  for  Mode  I  cracking.  However,  the  crack  path 
appears  to  be  mean  stress  and  material  dependent.  Hourlier 
et  al  tested  structural  steel  and  other  materials  and 
suggested  that  cracks  will  always  seek  the  mode  that  can 
generate  the  greatest  propagation  rate,  thus  adding 
additional  observational  information  on  the  direction  of 
crack  propagation  under  mixed  mode  loading. 

3. 1.2.2  Prediction  of  Mixed  Mode  Crack  Growth  Rates 

The  inclined  crack,  shown  in  Figure  3-22,  has  served  as  a 
useful  configuration  for  biaxial  fatigue  testing.  Research 
on  inclined  and  branched  cracks  is  extensive.  Kfouri  [3-56] 
presents  a  review  of  this  research.  Nearly  all  the  studies 
have  been  in  the  context  of  LEFM.  Briefly,  the  main 
criteria  are:  (a)  the  maximum  tangential  stress  criterion; 
(b)  the  criterion  of  local  symmetry;  (c)  the  minimum  strain 
energy  criterion;  and  (d)  the  maximum  energy  release  rate 
criterion.  A  criterion  based  on  the  maximum  normal  strain 
has  been  proposed  and  another  uses  the  crack-tip  opening 
displacement  (CTOD)  to  predict  the  direction  of  the  onset  of 
crack  extension. 

The  work  of  Erdogan  and  Sih  [3-57]  can  be  considered  among 
the  initial  efforts  in  the  area  of  mixed  mode  cracking. 

They  treated  the  problem  of  a  plate  under  uniform  tension  a 
with  an  angled  central  crack  of  length  2a  to  determine  the 
crack  direction.  They  assumed  that  the  crack  grows  in  a 
direction  eo  for  which  the  hoop  stress,  a9,  at  the  crack  tip 
is  maximum.  Williams  and  Ewing  [3-58]  modified  this  theory 
by  including  the  non-singular  terms  in  the  series  expansion 
for  better  correlation.  Later,  Finnie  and  Saith  [3-59] 
pointed  out  that  the  proposed  modification  neglected  the 
contribution  of  the  normal  stress  to  the  crack. 

Among  the  models  that  were  proposed  to  deal  with  mixed  mode 
cracking,  Swedlow's  Model  [3-60]  was  considered  as  a  basic 
model  for  the  biaxial  stress  investigation.  Swedlow 
proposed  a  compression  model  to  account  for  the  crack 
closing  and  frictional  effect  in  cracks  under  compression. 

He  used  Williams  analysis  [3-61]  to  derive  general 
expressions  for  the  stress  distribution  near  the  crack  tip. 
These  expressions  are: 
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Fracture  Angle  -Qb  degree 
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Figure  3-22.  Initial  Fracture  Angle  versus  the  Ratio  of  Stress  Intensity 
Ranges  Kjj  and  Kj 
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k„  19  3  e 

+  ~~ZZT  ~  cos  —  +  —  cos3  —  -a  sin8cos9 

J  2irr  4  2  4  2 

_  _ 

where : 

K,  =  aj  ?ra  (sin26) 

Kn  =  oj  jra(sinBc'  sB) 

a  =  a(cos2S  -  sinz6) 

where  a  is  the  semi-elliptical  crack  length  and  a  is  the  applied 
stress. 
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To  take  into  account  the  crack  closing  and  the  frictional 
effect,  Swedlow  suggested  that  a  should  not  be  replaced  by 
-a  to  express  a  compressive  model.  He  postulated  that  when 
Kt  has  negative  values,  it  means  that  the  crack  is  closing 
on  itself  and  by  such  action,  it  causes  friction  on  the 
surface  of  the  crack.  Thus,  in  compression,  Kx  can  be 
neglected  and  the  equations  can  be  expressed  as: 


’5  ©  3  9 

-  sin  -  +  -  sin3  - 

4  2  4  2 


K„ 

Or  =  ——— 

J  (2nr) 


+  <y(sinz0  -  jifsin20)  sinzB  +  atc os20 


3K1X 

a6  =  — — — - 
4y  (2?rr) 


0  0 

sin  —  +  sin3  — 
2  2 


+  a(cosz0  -  |i^sin20)  sinzB  +  atsinz0 


K 


1: 


T  rG  = 


J  (27rr) 


1  0  3  0 

—  cos  —  +  —  cos  — 
4  2  4  2 


where  a 


and 


+  o(cos0sin  0  -/LifCos20)  sinzB  +  Mtsin©cos0 
-1,  nt  is  the  coefficient  of  friction,  and 
Kn  =  a/“lra(sin  B cos  B  -  /ifsinzB) 

at  =  a(coszB  -  sinzB) 


These  equations  are  applicable  along  the  crack  faces.  The 
coefficient  of  friction  may  take  a  static  or  dynamic  value 
depending  on  the  motion  of  the  crack  faces.  Using  the 
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maximum  hoop  stress  criterion  of  Erdogan  and  Sih  and  taking 
into  account  the  crack  closing  and  the  friction  effect, 
Swedlow  was  able  through  his  model  to  predict  the  initial 
fracture  angle,  0O,  for  an  inclined  crack  under  uniaxial 
compression. 

Another  model  related  to  this  problem  is  one  proposed  by  Woo 
and  Ling  [3-62].  They  reviewed  a  number  of  fracture 
criteria  dealing  with  angled  crack  initiation  e0  and  its 
relation  to  crack  inclination  angle  B.  They  examined  the 
biaxial  loading  effect  on  the  parameter  used  in  each 
criterion  for  the  prediction  of  the  fracture  behavior  under 
mixed  mode.  In  their  experimental  work,  they  used  cruciform 
type  specimens  of  1/8  in.  thickness  with  an  initial  crack 
length  of  1-1/4  in.  The  material  used  in  their 
investigation  was  polymethyl-methacrylate  (PMMA) .  Their 
analysis  was  an  extension  of  Swedlow' s  approach  to  the  case 
of  biaxial  tension-compression  stressed  models.  This  was 
accomplished  by  modifying  Swedlow* s  near-tip  stress 
equations  to  the  following  forms: 


+  a„sin2e  +  atsin28 


3K„  e  e 

-  -  sin  —  +  sin3  —  +  atsin20 

4/  (2wr)  2  2 


+  ancos20  -  atsin20 


Ki _ 

4 J  ( 2nr ) 


e  e 

sin  -  +  sin3  - 
2  2 


+ 


K 


XX 


J  (2nr) 


19  3  e 
—  cos  —  +  —  cos3  — 
4  2  4  2 


-  otsinecos  e 


+  a„sin  0  cos  0  +  afcos20 


where 

a„  =  normal  compression  stress  on  the  crack  faces 

af  =  frictional  resistance  due  to  the  compressive 

stress  (=/ifan) 

Hi  =  friction  coefficient 

at  =  a(l-a)cos2fl 

They  stated  that  the  terms  in  the  above  equations  will  not 
appear  at  the  same  time  and  the  selection  of  the  proper 
terms  for  a  given  problem  was  listed  in  a  tabulated  form  as 
shown  in  Table  3-3.  On  the  basis  of  their  experimental 
work,  they  concluded  that  the  maximum  hoop  stress  criterion 
was  the  most  appropriate  one  to  use  for  both  the  open  and 
closed  crack  conditions. 

Zamrick  [3-63]  developed  a  biaxial  stress  model  to  predict 
crack  initiation  angle,  0O,  with  respect  to  an  inclined 
crack  angle  6: 


(o  +  tanzB) 


(1-a) 


1-3COS0 


sinO 


tanS 
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Selection  of 


3-42 


sin 


16  X 


e 

2 


3  L  tan©  J 


(1  +  atan2B)  =  0 


This  equation  is  the  proposed  model  for  relating  8  to  8  and 
a  in  a  biaxial  stress  state.  The  roots  of  the  equation  are 
denoted  by  0O. 

A  comparison  of  Zamrick's  model  to  the  Woo  and  Ling  approach 
is  shown  in  Figure  3-23.  The  two  approaches  are  in  close 
agreement.  However,  comparisons  to  experimental  data  shown 
in  Figure  3-23  are  not  as  convincing. 

Griffith's  energy  release  rate,  G(a),  which  is  equal  to  the 
crack  seprration  energy  rate,  GA(a) ,  when  the  material  is 
elastic  depends  on  the  direction  of  the  crack  extension 
given  by  the  angle  0.  When  the  extension  is  coplanar  (0  = 

0)  and  the  material  is  elastic  we  have: 

J  =  G (o)  =  GA  (O) 

where 

J  =  (K,z  +  K„2)/E  * 

and 

E  '  .  E/(l- u 2) 

The  energy  rates,  G(a)  and  G  (a) ,  can  be  expressed  as  the 
sum  of  the  Mode  I  and  Mode  II  components: 


G(a)  =  Gj(a)  +  Gn(a);  and 

GA(a)  =  G,A(o)  +  GjjA  (a) 


Estimates  for  G  (a)  for  non-planar  crack  extension  may  be 
obtained  from  finite  element  analyses  on  cracks  of  different 
orientations  undergoing  coplanar  extension.  Kfouri  and 
Miller  [3-64]  report  the  maximum  G  fracture  criterion  to  be 
applicable  to  linear  elastic  and  elastic-plastic  materials. 

Zamrick  and  Kitagawa  compared  these  criteria  to  experimental 
results  for  aluminum  and  steel  respectively  as  shown  in 
Figures  3-24  and  3-25. 
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Figure  3-23 .  Comparison  Between  the  Proposed  Method 
and  Ling  and  Woo  Method  for  a  «  -1.0 
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Figure  3-24.  -  Relationship  for  Biaxility  Ratio,  a  * 


Figure  3-25 


Crack  Angle,  6  (degrees) 
Relationship  for  Biaxiality  Ratio,  a 


The  agreement  of  any  one  approach  is  encouraging;  however, 
no  one  approach  is  significantly  highlighted  by  these 
correlations.  Also,  few  data  points  are  presented  for 
correlation.  Therefore,  further  development  of  these 
approaches  will  be  required  before  design  applications  can 
be  considered. 

3.2  FATIGUE  RESEARCH  ON  WELDED  DETAILS 

Few  multiaxial  fatigue  approaches  have  been  developed  for 
welded  details  due  to  the  number  of  variables  associated 
with  the  weld  process  and  geometry.  Nonetheless,  a  few 
researchers  have  attempted  correlations  for  relatively 
simple  weld  configurations  and  loadings  and  have  had  some 
success.  Rev  lew  of  this  important  research  is  presented 
next. 


3.2.1  Fatigue  Life  Estimates 

Munse  and  Stallmeyer  [3-65]  investigated  the  fatigue  of 
welded  plate  girders  with  panel  stiffeners  shown  in  Figure 
3-26.  Various  stress  formulations  were  tried  to  correlate 
S-N  data.  The  most  important  finding  is  that  most  of  the 
fractures  occurred  at  stiffeners  that  were  not  in  the  region 
of  pure  moment.  In  fact,  most  fractures  occurred  at 
stiffener's  ends  where  the  flexural  stresses  were 
considerably  lower  than  the  maximum  flexural  stresses. 

An  S-N  curve  is  presented  in  which  the  maximum  bending 
stress  on  the  extreme  fiber  at  the  location  of  failure  is 
used  as  the  ordinate.  These  data  are  shown  in  Figure  3-27. 
It  is  apparent  from  this  plot  that  the  maximum  bending 
stress  at  the  fracture  section  does  not  provide  a  consistent 
relationship  in  these  tests  either.  A  much  better 
correlation  is  obtained  when  the  data  are  analyzed  on  the 
basis  of  the  maximum  principal  tensile  stress  (including  the 
effect  of  shear)  at  the  point  of  failure.  Data  plotted  on 
the  latter  Vasis  are  presented  in  Figure  3-28.  Although 
there  is  still  scatter  in  the  test  results,  there  is  a  well- 
defined  scatter  band  for  the  entire  range  of  fatigue  lives. 
It  appeared  to  the  investigators  that  the  shear  and  flexural 
stress  in  the  web  of  the  welded  members  and  the  resulting 
principal  tensile  stresses,  along  with  the  stress 
concentrations  at  the  stiffener  welds,  produced  the  most 
critical  fatigue  conditions  in  the  weldments  with 
stiffeners.  Maximum  principal  tensile  stress  predicted 
fatigue  life  including  initiation  and  propagation. 

Later,  Moyar  and  Garg  [3-66]  re-analyzed  some  of  the  Munse 
data  using  other  fatigue  criteria.  They  calculated  the 
stresses  in  the  vicinity  of  fatigue  cracking  and  correlated 
the  data  based  on  maximum  principal  stress,  modified  range 
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Fig,  3—26  —Details  of  various  types  of  stiffener 
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of  shear  stress,  and  maximum  octahedral  shear  stress 
approaches.  The  results  of  these  tests  are  summarized  in 
Table  3-4.  Both  the  modified  range  of  shear  stress  and 
maximum  octahedral  shear  stress  criterion  provided  better 
correlation  with  observed  data  for  stiffened  welded  plate 
beams.  However,  this  analysis  is  based  on  a  few  tests  and 
correlations  are  not  very  convincing  in  the  low  cycle 
region. 

Lawrence  [3-67]  recently  proposed  an  approach  for  predicting 
crack  initiation  and  growth  in  simple  weldments  under 
combined  loading  (tensile  and  flexural  moment) .  This 
approach  is  based  on  the  ability  to  determine  a  stress 
fatigue  notch  factor  for  axial  and  bending  load  conditions 
in  the  vicinity  of  the  weld  toe  and  predicting  the  state  of 
stress  at  the  notch  by  superimposing  tensile  stress.  This 
approach  also  includes  the  effects  of  mean  stress  and 
residual  stress. 

Lawrence  proposed  that  for  long  fatigue  lives  and  constant 
amplitude  loading  conditions,  the  notch  root  stresses  are 
mostly  elastic  and  the  residual  stresses  can  be  considered 
not  to  relax.  Under  these  conditions,  the  Basquin  equation 
can  be  used  to  estimate  fatigue  life: 

S.Kf  =  (cr,  -  KfSm  -  or)  (2N,) b 


where  a  f  is  the  fatigue  strength  coefficient  ( a\  =  Su  +  50 
(ksi  units)),  Sm  is  the  remotely  applied  mean  stress,  ar  is 
the  notch-root  residual  stress  and  Kf  is  the  appropriate 
fatigue  notch  factor.  Expanding  the  mean  stress  (SJ  to 
include  both  (applied  or  induced)  axial  and  bending  mean 
stresses  (SnA  and  SmB) ,  and  considering  both  applied  and 
induced  cyclic  axial  and  bending  stresses  through  the 
following  equation: 

eff  A  B 

K  =  (l-x)K  +  xK 
fmax  fmax  fmax 


B  T 

x  =  S  /S 
a  a 


A  B 

where  K  and  K  are  the  worst-case-notch  fatigue  notch 
fmax  fmax 

factor  for  axial  and  bending  load  conditions,  respectively; 
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Table  3-4 


Re-Analysis  of  Plate  Girder  Data  by  Moyer 


SPEC. 

# 

AA- 

+SC 

OBS. 

CYCLES 

TO 

FAIL 

MAX.  PRIN. 
STRESS  THEORY 

MAX. SHEAR 
STRESS  THEORY 

MAX. OCT. 
STRESS  THEORY 

PREDICT  DIFF. 

PREDICT  DIFF. 

PREDICT  DIFF. 

29(0)B  0.6278 

0.1426  -0.4852 

0.3868  -0.2410 

0.3456  -0.2822 

30(0)B 

0.11467 

0.0731  -0.3736 

0.2836  -0.1631 

0.2355  -0.2112 

31(0)B  0.5236 

0.0744  -0.4492 

0.2884  -0.2352 

0.2399  -0.2837 

4l(0)B 

1.1103 

0.4438  -0.6665 

1.0261  -0.0842 

0.9188  -0.1915 

43(0)B 

1.2719 

0.8111  -0.4608 

1.2522  -0.0197 

1.1994  -0.0725 

MEAN  -0.4871 

-0.1486 

-0.2082 

STD 

.DEV.  0.1086 

0.0962 

0.0864 

COMPARISON  OF  PREDICTED  CYCLES  TO  FAILURE  (IN  MILLIONS) 

Of  Type  b  beams  using  properties  from  type 

C  BEAM  DATA 

AND  BIAXIAL  FATIGUE  THEORIES 
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B  T 

S  and  S  are  the  bending  stress  and  the  total  stress 
a  a 

amplitude,  respectively.  From  the  above,  one  can  derive  an 
expression  for  the  fatigue  strength  of  a  weldment  subjected 
to  axial  and  bending  mean  and  constant-amplitude  cyclic 
stresses: 


(1/Kf„x)  (Of  -(jr)  -S„-XSm* 

S.T  =  -  (2N1)b=  SF  (2Nj) b 

l-x(l-x) 

B  A 

where  X  is  the  ratio  of  K  to  K 

fmax  fmax 

If  the  assumptions  of  the  approach  are  valid,  then  this 
expression  should  predict  the  constant  amplitude  fatigue 
strength  at  long  lives  (NT  >  2x10®  cycles)  for  which 
initiation  is  thought  to  dominate  the  total  fatigue  life. 

The  variables  in  the  equation  can  be  divided  into  either 

A  B 

constants  of  known  quantities  (S  ,  S  ,  2N,)  and  random 

m  m 

A  A 

variables  (af ,  ar,  K  ,  x)  .  Of  the  random  variables  K 

fmax  fmax 

and  x  are  considered  as  determining  the  variation  in  fatigue 
strength.  Neither  material  property,  the  fatigue  strength 
coefficient  (a,')  which  is  proportional  to  the  UTS  or  the 
residual  stress  (ar)  which  is  equal  to  the  base  metal  yield 
strength  (Sy) ,  vary  greatly  for  one  material  and  welding 
process. 

To  show  the  improvement  associated  with  the  notch  stress 
approach,  Lawrence  calculated  mean  and  standard  deviations 
for  experimental  data  and  compared  it  to  Munse's  approach 
considering  nominal  field  stresses.  This  comparison  is 
discussed  in  Section  5.0. 

Yung  and  Lawrence  [3-68]  conducted  a  study  of  fatigue 
analysis  of  plate  and  tube  weldments  under  torsion  and 
bending  as  shown  in  Figure  3-29.  The  combined  bending  and 
torsion  loadings  studied  gave  lives  about  a  factor  of  three 
shorter  than  pure  bending  stresses,  showing  that  multiaxial 
stress  states  do  influence  the  fatigue  resistance  of 
weldments  despite  the  two-dimensional  nature  of  the  weld  toe 
notch  and  its  relative  insensitivity  to  pure  torsional 
loading.  Residual  stresses  were  found  to  play  an  important 
role  in  determining  the  portion  of  the  fatigue  life  period 
devoted  to  crack  initiation  and  early  growth.  Lawrence 
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found  that  the  multiaxial  fatigue  theories  of  Lohr  and 
Ellison  and  Kandil,  Brown  and  Miller  gave  best  estimates  of 
fatigue  crack  initiation  life,  as  shown  in  Table  3-5,  when 
the  notch  root  strains  could  be  determined  by  finite  element 
analysis  as  shown  in  Figure  3-30.  Estimates  of  initiation 
life  based  on  Basquin's  equation  modified  for  multiaxial 
stresses  (principal  stress  and  an  equivalent  stress)  gave 
predictions  for  long-life  high-cycle  fatigue  which  were 
within  a  factor  of  three  of  the  observed  total  lives,  as 
shown  in  Table  3-5. 

Siljander  et  al  [3-69]  recently  investigated  the  propor¬ 
tional  and  non-proportional  multiaxial  fatigue  of  tube-to- 
plate  weldments.  This  work  was  conducted  to  sup_  art  the 
development  of  a  multiaxial  fatigue  design  procedure  for 
built-up  plate  bridge  girders.  Bridges  experience  non¬ 
proportional  out-of-phase  stresses  as  vehicles  move  across 
the  bridge.  One  of  the  most  interesting  aspects  of  the 
study  is  the  analysis  of  the  principal  stresses  in  girders 
with  moving  loads  to  determine  representative  out-of-phase 
loadings.  They  correlated  their  data  using  various  damage 
criteria.  As  in  many  structural  codes,  principal  stress 
range  used  von  Vises  effective  shear  stress  range  and  a 
critical  shear  plane  approach  originally  developed  by 
Findley.  Of  these  approaches  the  critical  shear  plane 
approach  correlated  the  proportional  and  non-proportional 
data  quite  well  as  seen  in  Figures  3-31  through  3-33. 

Several  of  their  findings  in  achieving  this  degree  of 
correlation  are  noteworthy.  First,  the  test  specimens  were 
stress  relieved  to  eliminate  residual  weld  stresses.  This 
produced  a  significant  reduction  in  data  scatter.  Second, 
the  investigators  used  stress  concentration  factors  at  the 
weld  toe  that  are  developed  in  part  by  a  complex  3-D  finite 
element  analysis  with  the  added  advantage  of  having 
knowledge  of  weld  profiles,  information  that  is  not 
available  with  certainty  to  structural  designers. 

3.2.2  Fatigue  Propagation  in  Welded  Structures 

Multiaxial  fatigue  procedure  for  crack  propagation  have  been 
examined  because  weld  flaws  are  often  present  at  crack 
initiation  sites  as  discussed  by  Thayamballi  [3-70]  and  many 
others . 

Procedures  have  been  identified  for  estimating  the  stress 
intensity  factor  for  Mode  I,  KIt  propagation  in  complex 
stress  fields.  The  first  approach  identified  is  the 
influence  function  technique  presented  by  Bueckner  [3-71]: 


K  = 


h(x)a(x)dx 
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Table  3-5 


Research  on  Butt  Welds  &  Cruciform 
Fillet  Welds  by  Lawrence 


Statistical  S usury  of  tho  Departures  of 
Predicted  Lives  froa  Fatigue  Test  Data 


Munse's 

Miner's 

Gurney's 

RMS  Method 

I-P  Model 

No.  of 

Cases 

29 

29 

29 

13 

29 

fp 

1.061 

1.015 

0.894 

0.906 

1.016 

% 

0.124 

0.093 

0.081 

0.052 

0.067 

Fp  :  Mean  value  of  Fp;  Fp  -  - logl0P(N^Ct)°n  '  *  unlty  of  Fp  v*lu® 

Test 

represents  the  perfect  agreeaent  between  the  prediction  and  fatigue 
data. 

Qp^  :  Coefficient  of  Variation  of  Fp. 
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Figure  3-30 •  Finite  Element  Mesh  of  Tube-to-Plate  Welds 
(Ratio  of  Notch-Root  Element  Size  to  Notch 
Radius  =  1:4) 


Worst  case"  maximum  principal  stress  range  (MPa) 


Figure  3-31.  The  Maximum  ’'Worst  Case"  Local  Principal  Stress 
Range  versus  Cycles  to  Failure 


Bending  only  (A  =  <») 

Torsion  only  ( A  =  0) 

In-phase  combined  bending-torsion  (A  *  7.34) 
In-phase  combined  bending-torsion  (A  *  2.39) 
In-phase  combined  bending-torsion  ( A  »  1) 
Out-of-phase  combined  bending-torsion  (  A  *  2.39) 


Cycles  to  failure 


Figure  3-32. 


The  "Worst  Case"  Local  von  Mises'  Effective 
Stress  Range  versus  Cycles  to  Failure 


Nominal  Findley  shear  stress  range  (MPa) 


Bending  only  (A  =  °°) 

Torsion  only  (A  =  0) 

In-phase  combined  bending-torsion  ( A=  7.34) 
In-phase  combined  bending-torsion  (A  =  2.39) 
In-phase  combined  bending-torsion  ( A  =  1) 
Out-of-phase  combined  bending-torsion  ( A  =  2.39) 
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Figure  3-33.  Findley's  Equivalent  Nominal  Shear  Stress  Range  versus 
Cycles  to  Failure 


where 

h(x)  =  intensity  arising  from  a  unit  load  at 
location  x  on  the  crack  face 

a(x)  =  uncracked  stress  distribution 

This  approach  has  the  advantage  of  being  relatively  simple 
to  develop.  It's  disadvantages  include  the  assumption  that 
the  load  is  symmetric  about  the  crack  face. 

The  second  approach  for  estimating  stress  intensities  in 
complex  stress  patterns  from  geometric  configurations  is 
given  by  Gurney  [3-72]: 

MkMfMb 

K  =  -  a  J  *r  a 

0 

where 

K  -  a  yira  is  baseline  stress  intensity  for  the 
through  thickness  crack 

Mk  =  stress  concentration  factor  at  welded  toe 

Mf  =  front  face  correction  for  secondary  bending 
at  crack  tip 

M„  =  back  face  correction  for  finite  thickness 
plate 

0  =  geometric  factor  for  semi-elliptical  surface 

flaw 

The  advantages  of  this  approach  are,  again,  its  relative 
ease  of  application.  However,  factors  must  be  developed  for 
specific  crack  geometries  and  stress  conditions. 

The  third  approach  for  estimating  stress  intensity  factors 
is  by  finite  element  techniques  with  a  crack  element  in  the 
model.  This  approach  is  illustrated  in  Figure  3-34.  The 
advantage  of  this  approach  is  that  it  accounts  for  complex 
stress  patterns  and  crack  geometries  and  elastic  and  plastic 
conditions  around  the  crack  tip.  Disadvantages  of  this 
technique  are  that  it  is  relatively  time  consuming  and 
expensive  at  this  time.  This  situation  is  expected  to 
change  in  the  future  as  more  computational  techniques  are 
made  available  and  computer  capabilities  increase. 

Of  the  three  approaches,  the  later  has  been  proposed  most 
recently  for  applications.  Most  importantly,  crack  growth 
under  multiaxial  loading  has  been  shown  to  occur  along  a 


Finite  Element  Model  of  Structure  Adjacent  to  the  Crack  Tip 


Crack  Tip  Element 


Figure  3-34.  Illustration  of  Estimating  Stress  Intensity  by 
Finite  Element  Analysis 
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path  that  maximizes  Mode  One  component  K,  (Hellen  [3-73]) 
and  strain  release  rate  6  (McDonald  [3-74],  Haung  [3-75]). 
This  is  a  controversial  finding  with  conflicting 
experimental  results.  However,  if  validated  for  use  in 
marine  structures,  it  would  facilitate  development  of  an 
approach  for  predicting  multiaxial  fatigue  response  in 
welded  structures.  This  factor  alone  makes  prediction  of 
fatigue  crack  growth  in  welded  structures  worth  further 
consideration . 

There  is  evidence  that  multiaxial  stress  affects  the  fatigue 
life  of  welded  joints  as  indicated  by  Sih  [3-76].  The  crack 
trajectory  emanating  from  a  flaw  at  an  angle  to  the 
principal  stress  is  neither  exclusively  a  function  of 
principal  stress  nor  the  stress  component  perpendicular  to 
the  flaw.  This  is  also  true  for  welded  tubular  joints  by 
the  experimental  observations  of  McDonald  et  al  [3-74]. 
Furthermore,  there  is  evidence  [3-76]  to  indicate  that 
cracks  grow  along  a  curved  path  through  thickness  under 
chord-brace  intersections.  The  direction  of  a  growing 
fatigue  crack  in  a  complex  stress  field  has  been  the  subject 
of  several  hypotheses.  It  is  argued  by  Hellen  [3-73]  that 
the  crack  develops  in  such  a  way  as  to  maximize  the  Mode  I 
component  K,,  or  that  the  path  maximizes  the  strain  energy 
release  rate  G  [3-77,  3-78].  Insight  into  the  direction  of 
crack  growth  can  be  obtained  by  considering  the  effect  of  a 
small  angle  kink  on  the  crack  tip.  The  local  Mode  I  and 
Mode  II  stress  intensity  factors  on  the  kink  are  denoted  K, 
and  Kn  can  be  expressed  in  terms  of  the  Mode  I  and  Mode  II 
components  Kj  and  Kn  on  the  main  crack  following  the  work  of 
Bilby  [3-79];  Hussain  et  al  [3-80];  and  Masahiro  et  al  [3- 
81].  Although  the  analyses  differ  on  detail  they  give 
broadly  similar  results.  The  angle  of  the  crack  for  a 
tubular  joint  weld  is  shown  in  Figure  3-35.  The  angle  of 
the  kink  which  maximizes  Ki  and  makes  K„  zero  is  close  to  15 
degrees  for  crack  length  to  thickness  ratio  a/T=.2,  20 
degrees  for  a/T=.6,  and  approximately  55  degrees  for  a/T=.9. 
The  local  strain  energy  release  rate  G  for  the  kink  can  be 
determined  by  combining  the  local  Kt  and  K„,  or  be  obtained 
by  aft  axis  virtual  crack  extension.  For  the  case  shown  in 
Figure  3-35  crack  paths  which  maximize  Kx  and  G  are  similar. 
These  are  plotted  as  three  piece  wise  linear  segments  and 
compared  with  the  path  observed  experimentally  by  Rice  as 
presented  by  Haung  et  al  [3-82]. 

Rhee  [3-76]  has  presented  two  techniques  for  computing 
stress  intensity  factors  for  welded  structural  details.  The 
first  method  is  based  on  development  of  an  effective  stress 
intensity  factor  (K.)  for  mixed  mode  crack  propagation  and 
the  second  is  by  using  a  special  finite  element  program  for 
estimating  stress  intensity  factors. 
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EXPERIMENTAL  DATA 


Predicted  Crack  Path  at  a  Welded  Tubular  Intersection 


The  equivalent  stress  intensity  factor,  K.,  is  converted 
from  the  appropriate  energy  release  rate,  G: 


K. 


=  G  = 


E 


1— u 

E 


K,2  +  K„2  + 


Kni 


1-u 


This  equivalent  stress  intensity  factor  represents  the  crack 
driving  force  which  incorporates  all  the  KI#  Kn  and  Km 
components  of  a  mixed  mode  problem.  For  a  mixed  mode 
problem,  any  single  component  of  the  stress  intensity 
factors  cannot  represent  the  crack  driving  force  properly 
since  crack  propagation  is  a  combined  result  of  the 
contributions  of  all  three  stress  intensity  factors  of  the 
problem.  A  comparison  between  K,  and  any  individual 
component  of  K, ,  K„  and  Kin  can  indicate  the  contribution  of 
such  an  individual  component  to  crack  propagation.  For  the 
brace  tension  cases,  it  is  apparent  from  Figures  3-36 
through  3-39  that,  throughout  the  crack  front,  the  fatigue 
crack  propagation  will  dominate  the  Mode  I  stress  intensity 
factor  K,.  However,  for  the  in-plane  bending  cases,  the  Kj 
contribution  to  K.  is  negligible  except  for  a  small  length 
near  the  surface.  Therefore,  if  the  considered  flaw 
geometry  is  practical  for  a  fatigue  load  system  which 
consists  of  a  significant  contribution  of  in-plane  bending, 
a  fatigue  crack  growth  analysis  should  consider  all  three  — 
K, ,  K„  and  KIM. 

A  crack  driving  force  parameter  such  as  K.  has  yet  to  be 
established  as  a  parameter  through  which  the  fatigue  crack 
growth  behavior  of  a  mixed  mode  flaw  can  be  directly 
calculated  in  a  manner  similar  to  the  Kj  of  a  Mode  I 
problem.  Expressions  similar  to  the  present  K.  have  been 
studied  for  limited  mixed  mode  problems  [3-83,  3-84]  such  as 
problems  with  K,  and  Kn  mixed  and  with  K,  and  Kn,  mixed. 
However,  few  fully  mixed  mode  problems,  such  as  those 
present,  have  been  studied. 

Rhee  has  also  presented  a  method  to  calculate  the  stress 
intensity  of  an  elliptical  flaw  using  a  specialized  finite 
element  computer  approach.  This  method,  called  the  Finite 
Element  Alternating  method,  can  solve  the  problem  of  an 
elliptical  flaw  under  arbitrary  normal  and  shear  stresses 
which  can  be  expressed  in  higher-order  polynomials  [3-85 
through  3-87].  The  Finite  Element  Alternating  method 
iterates  through  the  following  steps. 

1.  Using  the  solution  of  an  elliptical  crack  in  an 
infinite  solid  and  subjected  to  a  constant 
pressure,  the  surface  traction  ayy,  ryi  and  ryz  are 
computed  on  y=-h  plane. 
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Figure  3-36. 


Stress  Intensity 
Factors  of  Weld 
Toe  Surface  Flaw 
(In-Plane  Bending) 
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Figure  3-37. 
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Figure  3-38.  Stress  Intensity 
Factors  of  Weld 
Toe  Surface  Flaw 
(Axial  Tension) 
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Figure  3-39.  Stress  Intensity 
Factors  of  Weld 
Toe  Surface  Flaw 
(Out-of-Plane 
Bending) 
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2.  The  residual  tractions  acting  on  y=-h  plane  are 
then  freed  by  applying  opposing  surface  tractions 
on  the  plane  of  an  uncracked,  semi-infinite  solid. 
The  resulting  residual  normal  stress,  azz,  on  the 
crack  surface  location  is  calculated. 

3.  The  residual  normal  stress,  azz,  on  the  crack 
surface  is  removed  by  an  opposing  stress.  This 
opposing  stress  results  in  the  stress  intensity 
factors  along  the  crack  front  and  the  surface 
tractions  on  y=-h  plane,  which  can  be  calculated 
using  Segedin's  potential  function  [3-88]. 

4.  Steps  2  and  3  are  repeated  until  the  resulting 
residual  stresses  on  the  crack  tip  surface  become 
negligible  when  compared  to  the  applied  stress  of 
r.  The  result  of  all  the  iterations  are 
superimposed  to  obtain  the  stress  intensity 
factors  for  the  elliptical  crack  of  the  original 
problem. 

This  approach  has  been  successfully  applied  to  obtain 
accurate  stress  intensity  factors  in  finite  thickness  plate 
[3-87],  in  pressure  vessels  [3-89]  and  for  quarter- 
elliptical  corner  crack  emanating  from  a  pinhole  in  plates 
and  in  aircraft  attachment  lugs  [3-90]. 

In  the  finite  element  alternating  method,  the  solutions  of 
an  elliptical  flaw  in  an  infinite  body  are  obtained  through 
the  Trefftz's  potential  function  method  [3-91].  In  this 
procedure,  extensive  numerical  calculations  are  involved  to 
evaluate  the  various  forms  of  elliptical  integrals  required 
for  the  solutions.  The  finite  element  solution  method  is 
used  to  obtain  the  solutions  of  the  uncracked  body  under 
arbitrary  boundary  tractions.  The  boundary  tractions,  which 
result  from  the  above  elliptical  flaw  problem,  are  first 
converted  into  the  boundary  nodal  forces  of  the  finite 
element  model  without  a  crack  and  then  applied  as  the 
external  loads  to  calculate  the  stresses  on  the  crack 
location.  These  solutions  are  used  following  the  sequences 
discussed  earlier  to  obtain  the  stress  intensity  factors. 

Another  finite  element  based  approach  for  calculating  crack 
propagation  in  welded  details  was  presented  by  Haung  [3-92]. 
His  approach  is  based  on  the  work  of  Rice  and  Levy  [3-93] 
where  the  crack  is  represented  by  a  series  of  generalized 
line  springs  which  act  across  a  discontinuity  in  a  thin 
shell.  The  approach  was  originally  developed  for  Mode  I 
loading  but  has  been  further  generalized  by  Parks  [3-94]  and 
Desvaux  [3-95]  to  incorporate  both  Mode  II  and  Mode  III 
loadings  and  is  implemented  in  the  finite  element  code 
ABAQUS  [3-96]. 
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Haung  [3-92]  presented  results  comparing  the  line  spring 
model  to  three-dimensional  finite  element  analysis  and  to 
experimental  crack  growth  rates  for  a  chord  and  brace 
intersection. 

A  fracture  mechanics  fatigue  analysis  method  explicitly 
considers  the  flaw  geometry,  which  is  one  of  many  critical 
characteristic  of  the  problem.  The  necessity  to  consider 
the  flaw  geometry  in  a  fatigue  analysis,  in  turn,  can  make  a 
rigorous  crack  growth  analysis  prohibitively  expensive.  The 
effort  required  for  the  evaluation  of  the  stress  intensity 
factors  can  be  significant,  even  for  a  simple  flaw  geometry. 
Especially  for  a  surface  flaw,  which  is  a  common  form  of 
weld  toe  defects  in  a  marine  structure  and  whose  stress 
intensity  factors  can  only  be  practically  evaluated  through 
a  numerical  method  such  as  the  finite  element  method,  the 
investment  required  for  the  evaluation  of  the  stress 
intensity  factor  solution  is  enormous. 

To  perform  a  rigorous  fatigue  crack  growth  analysis  of  a 
surface  flaw,  which  changes  its  size  and  shape  continuously 
under  fatigue  loading,  the  stress  intensity  factors  of  a 
group  of  flaws  with  different  dimensions  are  required  to 
cover  a  certain  range  of  flaw  shapes  and  sizes.  For  the 
stress  intensity  factor  calculation  using  the  finite  element 
method,  which  is  the  most  popular  at  present,  the  majority 
of  the  analysis  effort  is  devoted  to  modeling  the  finite 
element  meshes. 
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4.0  APPLICATION  OF  MULTIAXIAL  FATIGUE  RESEARCH  TO  MARINE 
STRUCTURES 


Marine  structures  are  subjected  to  various  loads  and  local 
effects  that  are  combined  in  structural  details.  These 
loadings  contribute  to  fatigue  in  the  details.  In  this 
section,  the  stress  characteristics  described  in  Section  2.1.1 
and  multiaxial  fatigue  approaches  presented  in  Section  3.0  are 
compared  to  identify  candidate  multiaxial  fatigue  approaches 
for  marine  structures  and  to  help  identify  areas  where 
additional  research  is  required.  Table  4-1  has  been  prepared 
to  summarize  the  proposed  approaches  identified  during  the 
course  of  this  project  and  to  highlight,  for  quick  reference, 
areas  in  which  gaps  in  technology  exist.  Further  discussion 
of  the  applicability  of  multiaxial  fatigue  research  to  marine 
structures  follows. 

4.1  SHIP  STRUCTURE  APPLICATIONS 

Returning  to  our  examples  of  representative  structural  details 
in  ship  structures,  we  will  present  relevant  methods  that  have 
been  proposed  for  predicting  multiaxial  fatigue  response. 

The  first  detail  considered  is  the  cutout  in  a  transverse  web 
frame.  The  cutout  is  to  allow  longitudinal  panel  stiffeners 
to  pass  through  continuously  in  order  to  carry  longitudinal 
primary  hull  bending  stresses.  There  have  been  studies 
conducted  on  this  detail  [4-1]  to  determine  the  state  of 
stress.  As  discussed  in  Section  2.0,  there  are  biaxial 
nominal  stresses  adjacent  to  the  cutout  caused  by  distributed 
lateral  load  on  the  girder  resulting  in  flexural  bending  and 
shear  stresses.  However,  at  the  edge  of  the  cutout,  where 
cracks  are  known  to  propagate,  the  stress  is  axial  tension  or 
compression.  The  principal  stress  field  has  mean  stress 
components  from  hydrostatic  loads  and  varying  principal 
stresses  from  variations  in  hydrostatic  loads  as  waves  are 
encountered.  Therefore,  on  a  detailed  level,  the  stresses  are 
not  at  all  multiaxial,  but  on  a  nominal  stress  level,  the 
stresses  are  biaxial  as  indicated  by  Munse  and  others.  Munse 
suggested  shear  stress  as  the  correlating  nominal  stress  for 
this  flexural  member  with  cutouts. 

The  second  example  detail  is  a  hatch  corner  cutout.  The  hatch 
corner  is  a  source  of  stress  concentration.  The  state  of 
stress  in  the  vicinity  of  this  cutout  is  biaxial  resulting 
from  combined  axial  torsional  hull  loading.  However,  as  with 
the  previous  cutout,  the  stress  state  on  the  face  of  the 
cutout  is  axial.  Research  by  Chen  [4-2]  and  Munse  [4-3] 
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TABU  4-1  (continued) 


indicates  that  the  axial  stress  is  sufficient  for 
characterizing  the  fatigue  response  of  the  hatch  corner 
cutout. 

The  third  example  detail  is  a  Center  Vertical  Keel  (CVK) .  The 
CVK  is  subjected  to  varying  axial  load  along  its  axis,  lateral 
hydrostatic  load  with  a  mean  component  and  often  internal 
static  or  dynamic  (sloshing  liquids)  loads  opposing  the 
hydrostatic  load.  As  indicated,  mean  stresses  are  present  and 
often  are  random  over  the  voyage  from  one  voyage  to  another. 
The  phase  relationship  of  the  axial  and  lateral  load  varies  in 
a  random  nature. 

Munse  [4-4]  suggested  that  the  principal  stress  correlates  the 
fatigue  life  (initiation  and  propagation).  Moyar  and  Garg  [4- 
5]  re-analyzed  a  few  tests  data  and  found  the  octahedral  shear 
stress  and  maximum  shear  stress  correlated  data  best.  More 
interestingly,  Moyar  included  a  mean  stress  correction  based 
on  the  Goodman  Diagram. 

Recently,  Silijander  [4-6]  introduced  a  new  approach  to 
predict  the  multiaxial  fatigue  response  in  plate  bridge 
girders  starting  with  an  analysis  of  the  proportionality  of 
principal  stresses  in  the  girder  web  and  flange  intersection. 
They  found  the  critical  plane  approach  of  Findley  to  best 
correlate  the  data.  However,  as  discussed  in  Section  3.0,  the 
application  of  this  technique  depends  on  the  availability  of 
information  on  weld  stress  concentration  factors  and  cyclic 
material  properties.  This  type  of  information  may  be  known  by 
laboratory  researchers  but  is  not  often  available  to 
structural  designers. 

The  transverse  butt  weld  is  a  fourth  example  and  a  very  common 
weld  detail,  there  is  little  correlatable  data  for 
applications  similar  to  those  in  ship  structures.  The 
transverse  butt  weld  is  generally  located  in  a  biaxial  stress 
field  with  varying  mean  stress  and  fluctuating  principal 
stress  proportionality.  Tests  were  conducted  on  butt  welded 
HY-80  plates  subjected  to  lateral  load?  however,  the  biaxial 
nature  of  stresses  was  not  considered  —  only  stresses  normal 
to  the  weld  toe.  Similarly,  Silijander  [4-6]  tested  butt 
weld's  subjected  to  axial  and  bending  loads  and  again  only 
considered  the  maximum  stress  normal  to  the  weld  toe. 

Because  the  transverse  butt  weld  detail  is  quite  common,  the 
stress  field  easy  to  estimate  and  a  large  quantity  of  fatigue 
data  available,  the  investigators  applied  biaxial  correlation 
techniques  to  HY-80  tests  data.  Details  of  these  calculations 
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are  presented  in  Appendix  B.  Equivalent  shear  stress 
techniques  modified  for  mean  stress  effects  provided  very  good 
predictions  of  fatigue  life  of  a  laterally  loaded  plate  with  a 
butt  weld  based  on  uniaxial  data.  It  was  found,  however,  that 
multiaxial  approaches  based  on  cyclic  strain  fatigue 
properties  were  not  easy  to  apply  because  at  this  time  these 
properties  are  generally  not  known  for  ships  and  offshore 
structures.  Additionally,  stress  concentration  factors  for 
weld  geometry  have  a  large  influence  on  fatigue  response 
depending  on  the  weld  profile  and  quality.  The  HY-80  plate 
butt  had  a  quality  weld  and  low  weld  profile  resulting  in  a 
stress  concentration  factor  of  about  1.15,  not  a  typical 
stress  concentration  for  welds. 

The  foregoing  discussion  illustrates  the  variability  of  load 
characteristics  and  geometry  of  ship  structures  and  the 
different  multiaxial  fatigue  approaches  that  are  applicable  to 
each  type  of  detail. 

Additionally,  only  a  few  data  points  were  generalized  to 
support  the  research  and  are  insufficient  for  use  in 
developing  confidence  in  statistical  parameters  (bias  and 
scatter)  needed  for  reliability  analysis. 

4.2  OFFSHORE  STRUCTURE  APPLICATIONS 

Currently  there  are  two  distinct  approaches  for  predicting 
fatigue  response  in  offshore  structures.  The  first  is  a 
fatigue  life  approach  based  on  S-N  curves  [4-7]  developed  for 
various  joint  configurations.  This  approach  is  the  most 
widely  used  in  the  industry  and  uses  the  maximum  stress  normal 
to  the  weld  toe  and  neglects  interaction  effects  for  out-of- 
plane  joint  configurations.  The  British  DOE  rules  [4-8] 
address  complex  stress  distributions  by  correlating  fatigue 
data  using  maximum  principal  stresses  at  detail  "hotspots"  as 
discussed  in  Section  2.2.1  of  this  report.  Yung  and  Lawrence 
[4-9]  have  used  maximum  shear  stress  and  octahedral  shear 
stress  to  correlate  tube  and  plate  high  cycle  fatigue  response 
under  multiaxial  loading. 

The  second  approach  for  fatigue  response  prediction  that  is 
gaining  wider  acceptance  is  the  fracture  mechanics  approach. 
Multiaxial  stress  distributions  are  being  considered. 

However,  validation  efforts  are  very  limited  in  number,  detail 
stress  distributions  must  be  known,  and  the  effort  is  more 
time-consuming  than  fatigue  life  approaches. 
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Rhee  presented  estimates  of  stress  intensity  factors  for 
multiaxial  loading  based  on  finite  element  analyses  and  using 
an  equivalent  stress  intensity  factor  as  discussed  in  Section 
3.2.2.  Figures  shown  in  Section  3.2.2  present  the  resulting 
stress  intensity  factor  solutions  of  two  flaw  geometries. 

These  plots  show  that  the  in-plane  bending  load  cases  are 
predominately  mixed  mode.  All  the  solutions  along  the  crack 
front  location  where  the  negative  Mode  I  stress  intensity 
factor  is  developed  (for  the  in-plane  bending  cases)  are  not 
accurate  since  a  negative  K,  indicates  the  crack  surface 
contact  and  penetration  which  developed  in  the  analysis  were 
physically  impossible.  These  negative  values  could  be 
eliminated  with  a  proper  multicomponent  contact  algorithm.  A 
comparison  between  K.  and  any  individual  component  of  Kj,  K„ 
and  Kni  can  indicate  the  contribution  of  such  an  individual 
component  to  crack  propagation. 

While  Rhee's  K,  is  attractively  simple,  the  approach  is 
obviously  limited  to  linear  elastic  fracture  mechanics 
assumptions  in  the  high  cycle  range  of  crack  growth.  It 
obviously  does  not  account  for  crack  tip  plasticity  observed 
in  mixed  mode  tests  as  discussed  earlier. 

Huang  [4-10]  presented  a  more  simplified  approach  for 
calculating  mixed  mode  crack  growth  in  welded  tube  joints 
using  finite  element  analysis,  line  spring  elements,  and  the 
energy  release  rate  G.  Haung  compared  his  approach  to 
experimental  results  and  found  good  agreement  for  cases  where 
Mode  I  dominates  and  cracks  are  not  long,  0.8  >  (a/T)  >  0.2. 
This  technique  is  again  limited  to  basic  limitations  of  linear 
elastic  fracture  mechanics  and  does  not  account  for  mean  and 
residual  stresses,  random  effects  or  interactions  with 
corrosive  environments.  Additionally,  there  are  conflicting 
reports  on  the  accuracy  of  assuming  Mode  I  dominant  crack 
growth . 

Another  type  of  offshore  structure  in  which  multiaxial  fatigue 
considerations  are  important  is  riser  connectors.  These 
connector  systems,  illustrated  in  Figure  4-1,  are  often 
subject  to  high  pressure  loads  and  tension  and  torsion  loads 
that  fluctuate  about  a  mean  stress.  There  are  several 
multiaxial  fatigue  approaches  that  are  good  candidates  for 
this  type  of  loading,  especially  those  in  which  tension  and 
torsion  type  loads  dominate.  The  automotive  and  heavy 
equipment  industries  have  been  very  active  in  research  on 
tension  and  torsion  multiaxial  fatigue  with  several 
specialized  approaches  investigated  by  Tipton  and  Nelson  [4- 
11]  which  have  been  described  previously.  However,  these 
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Figure  4-1.  Tension  Member  and  Connector  with  High  Pre-Load 
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approaches  should  be  subject  to  verification  using  offshore 
materials  and  the  exact  geometries  used  in  connectors. 

The  equivalent  stress  techniques  have  been  used  most  often 
because  they  are  easiest  to  apply  and  required  material  data 
is  available.  Of  the  equivalent  stress  techniques  reviewed, 
maximum  shear  stress  techniques  generally  predict  crack 
initiation  and  equivalent  principal  stress  is  most  useful  for 
predicting  fatigue  life  where  crack  propagation  life 
dominates . 

However,  from  the  foregoing  discussion  and  the  summary  in 
Table  4-1,  there  is  no  universal  technique  for  predicting 
multiaxial  fatigue  response  for  each  varied  application  in 
marine  structures.  This  is  because  marine  structural  details 
vary  in  geometry,  applied  load,  stress  state  and  fabrication 
procedures.  Unlike  pressure  vessels  or  rotating  machinery 
(where  these  approaches  have  been  applied)  there  are  a  great 
number  of  variables  for  each  structural  detail. 

4.3  FACTORS  INFLUENCING  MULTIAXIAL  FATIGUE  RESPONSE  IN  MARINE 
DETAILS 

The  many  factors  which  influence  fatigue  response  in  marine 
structures  were  discussed  in  Section  2.1.  The  factors  that 
have  influence  in  fatigue  response  of  welded  details  are 
discussed  here  along  with  the  research  conducted  to  determine 
their  influence  on  multiaxial  fatigue  response  and  prediction. 
These  factors  include  stress  gradients,  stress  proportionality 
and  phasing,  elastic  and  plastic  strain  relationships,  mean 
and  residual  stress,  random  loading,  corrosion,  thickness 
effects,  materials  and  fabrication  procedures.  A  data  set  has 
been  identified  where  several  of  these  effects  were 
investigated.  The  data  set  is  from  welded  HY-80  and  HY  series 
steel  plates  subjected  to  lateral  pressure  load.  Effects  from 
various  sources  are  easier  to  separate  from  one  another  by 
comparing  different  data  sets.  The  intent  of  this  review  is 
to  illustrate  the  influence  these  factors  have  on  multiaxial 
fatigue  response. 

Tests  of  a  variety  of  welded  plate  laboratory  specimens  of  HY 
steels  are  also  summarized  in  this  study.  The  general 
configurations  of  these  various  laboratory  specimens  are 
presented  in  Table  4-2.  A  summary  of  the  fatigue  resistance 
of  these  various  members  is  presented  in  Table  4-3  and  Figures 
4-2  through  4-19.  Table  4-4  gives  a  brief  description  of  the 
test  members  included  in  the  study  of  HY  steels. 
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Table  4-2.  General  Configuration  of  Various  Specimen  Types 


Steel 

Refs/Notes 

HY-80 

HY-100 

HY-130/ 

150 

4-12,  4-13, 
4-14,  4-15, 
4-16,  4-17, 
4-19 

HY-80 

HY-100 

4-12,  4-13 
4—14,  4-15 
4-16,  4-17, 
4-18,  4-19, 
4-20,  4-21 

HY-80 

HY-100 

4-12,  4-13, 
4-17 

HY-80 

HY-130/ 

150 

4-54 

IV(A)-Slmply  sup¬ 
ported  on  4  edges 

IV(B)- Simply  sup¬ 
ported  on  2  edges, 
free  on  2  edges 

t  -  1  1/2,  1  5/8 
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Table  4-2.  General  Configuration  of  Various  Specimen  Types  (cont'd.) 


Type 

Specimen 

29  1/2"  — 

■ 

32" 

V 

■ 

j 

J 

l _ _ _ i 

\f 

I11H11H1 

It 

t  =  1,  1  11/16 

29  1/2*  — ► 

f 

j 

i 

! 

i 

36" 

f 

or 

cc* 

VI 

| 

I 

i 

OO 

j 

Steel 

Refs /Notes 

HY-80 

4-22,  4-23, 

HY-130/ 

150 

4-24,  4-25, 
4-26,  4-27 
4-28,  4-29, 
4-30,  4-53, 
4-32 

V(A)— Simply  sup¬ 
ported  on  4  edges 

V(B)-Slmply  sup¬ 
ported  on  2  edges, 
free  on  2  edges 

HY-80 

4=33,  4-22, 

HY— 100 

4-34,  4-35, 
4-36,  4-37, 
4-38,  4-31, 
4-41,  4-42, 
4-45,  4-30, 
4-55 

V!(A)-S!mply  sup¬ 
ported  on  4  edges 
(29  1/2  x  57  1/2 
x  1  1/2) 

VI(B)— Simply  sup¬ 
ported  on  2  edges, 
free  on  2  edges 

VI(C)-Cast  tees 
(No  welding) 

HY-80 

4-33,  4-47 
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Table  4-2.  General  Configuration  of  Various  Specimen  Types  (cont'd.) 


Steel 

Refs/Notes 

MS 

(71.2ksi 
tensile) 
HY-80 
HY-100 
HY— 140 

4-49,  4-50 
4-51,  4-52 

For  2"  Plate  box 
dimensions  are 

12  1/4"x12  1/4" 
x  28  3/4" 

HY-80 

4-12,  4-13 

HY-80 

4-12 

HY-80 

HY-100 

4-12,  4-13 
4-17,  4-18 

'2 
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Table  4-2.  General  Configuration  of  Various  Specimen  Types  (cont'd.) 
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FATIGUE  STRENGTHS  OF  VARIOUS  SPECIMEN  TYPES 
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TABLE  4-4 


HY  STEEL  TEST  MEMBERS 

USED  TO  INVESTIGATE  MULTIAXIAL  FATIGUE  RESPONSE 


'ype  I .  Axially  loaded  flat  plate  specimens  in  thicknesses  of 
/2,  3/4  and  1-1/2-inch.  Materials  included  are  HY-80,  HY-100 
ind  HY-130/150.  The  members  have  been  studied  under  the  stress 
atios  of  R=0  and  R=-l.  The  results  of  the  tests  are  found  in 
'igures  4-2  and  4-3. 

’ype  II.  These  members  are  again  flat  plate  members  with 
:ransverse  butt  welds  and  studied  under  axial  loading  in 
thickness  of  3/4  and  1-1/2-inches.  Members  included  in  the 
itudies  are  of  HY-80  and  HY-100  steels  and  were  subjected  to 
itress  ratios  of  R=0  and  R=-l,  and  R=+l/2.  The  results  of  the 
:ests  are  found  in  Figure  4-4. 

ype  III.  These  are  also  flat  plate  specimens  subjected  to  axial 
.oadings  and  fabricated  with  transverse  stiffeners  welded  to  both 
iurfaces  of  the  plate  with  full  penetration  welds.  Material  was 
.-1/2  inches  thick  and  test  members  fabricated  of  HY-80  and  HY- 
.00  steels. 

ype  IV.  The  type  IV  specimens  consisted  of  flat  rectangular 
)lates  that  were  subjected  to  transverse  pressure  thereby 
jroducing  a  flexural  loading  of  the  plates.  Type  IVA  specimens 
te re  simply  supported  on  four  edges.  In  the  case  of  type  IVB 
members,  the  plates  were  simply  supported  on  two  opposite  edges 
ind  free  on  the  remaining  two  edges.  In  this  instance,  most  of 
:he  tests  have  been  conducted  on  HY-80  or  HY-130/150  steels.  The 
'esults  of  the  test  can  be  fund  in  Figures  4-5  and  4-6. 

ype  V.  These  are  flat  plate  members  with  transverse  butt  welds 
:hat  were  subjected  to  transverse  bending  loads  and  supported 
iither  on  four  or  two  edges  of  the  plate.  In  these  studies, 
jrimarily  conducted  on  HY-80  and  HY-130/150  steels,  defects  such 
is  undercut  and  lack  of  fusion  were  introduced  in  the  welds.  The 
•esults  of  these  various  studies  can  be  found  in  Figures  4-7,  4-8 
ind  4-9. 

ype  VI .  These  are  flat  type  specimens  which  were  simply 
;upported  on  either  two  or  four  edges  of  the  plates  and  subjected 
:o  a  uniform  transverse  pressure.  However,  in  this  instance  a 
itiffener  was  welded  to  the  surface  of  the  plate  with  full 
>enetration  welds.  In  this  instance  most  of  the  tests  were 
:onducted  on  HY-80  and  HY-100  steels.  However,  a  variety  of 
rariables  have  been  included  in  the  studies  and  cast  steel 
ipecimens  have  also  been  studied.  The  results  of  these  various 
itudies  can  be  found  in  Figures  4-10  and  4-11. 
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TABLE  4-4  (continued) 


HY  8TEEL  TEST  MEMBERS 

USED  TO  INVESTIGATE  ULTIAXIAL  FATIGUE  RESPONSE 


Type  VII.  The  Type  VII  members  were  beam-type  plates  on  which  a 
transverse  stiffener  had  been  welded  with  full  penetration  welds. 
Again,  these  members  are  fabricated  of  HY-80  stee.l  and  subjected 
to  a  flexural  loading  which  provided  a  stress  ratio  of  R=l.  In  a 
limited  number  of  tests  a  stress  ratio  R=0  was  also  included. 

The  toes  of  the  fillet  welds  of  these  members'  were  treated  in 
various  ways  including  mechanically  peening,  grinding,  and  shot- 
peening.  The  results  of  these  various  tests  will  be  found  in 
Figures  4-12  and  4-13. 

Type  VIII.  These  specimens  consisted  of  box  sections  fabricated 
of  plates.  The  boxes  were  fabricated  of  mild  steel,  HY-100,  and 
HY-140  steels.  The  various  box  members  were  subjected  to  a 
stress  ratio  R=0.  Only  a  limited  number  of  tests  of  this  type 
have  been  conducted.  Nevertheless,  results  of  these  various 
tests  are  summarized  in  Figure  4-14. 

Type  IX.  These  members  consisted  of  flat  plates  specimens  with  a 
longitudinal  full  penetration  butt  weld  in  the  member.  Results 
of  tests  on  HY-80  steel  are  presented  in  Figure  4-15. 

Type  X.  These  members  consisted  of  flat  plates  with  transverse 
attachments  on  the  surfaces  of  the  plate  and  attached  with  fillet 
welds.  The  tests  on  HY-80  steel  are  summarized  also  in  Figure  4- 
15.  However,  only  a  very  limited  number  of  tests  were  conducted 
on  such  members. 

Type  XI.  These  specimens  consisted  of  flat  plates  subjected  to 
axial  loadings  and  fabricated  with  a  transverse  attachment  on  one 
surface  attached  with  full  penetration  welds.  Material  was  HY-80 
or  HY-100  steel  and  the  members  were  subjected  to  stress  ratios 
of  R=-l,  R=0  and  R=l/2.  Results  of  these  various  tests  are 
summarized  in  Figure  4-16. 

Type  XII.  This  specimen  consisted  of  a  flat  plate  with 
longitudinal  attachments  fillet  welded  to  the  surface  of  the 
plate.  Only  one  tests  of  this  type  was  conducted  on  HY-80  steel 
and  is  summarized  in  Figure  4-17. 

Type  XIII.  These  specimens  were  fabricated  as  a  cruciform  T- 
joint.  The  flat  plates  were  welded  to  a  transverse  bar  with  full 
penetration  welds  and  the  transverse  bar  subjected  to  a  loading 
normal  to  the  surface  of  the  bar.  Both  HY-80  and  HY-100  steels 
were  used  in  these  T-joints.  The  results  of  the  few  tests  of 
this  type  are  presented  also  in  Figure  4-17. 
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table  4-4  (continued) 


HY  STEEL  TEST  MEMBERS 

USED  TO  INVESTIGATE  MULT I AXIAL  FATIGUE  RESPONSE 


Type  XIV.  These  specimens  were  flat  plate  members  that  were 
subjected  to  a  transverse  pressure  and  simply  supported  on  four 
edges.  These  were  smaller  than  the  specimens  discussed  above 
under  Types  IV,  V  and  VI.  They  were  also  thinner  plates,  only 
3/4-inch  thick.  Nevertheless,  they  were  loaded  in  a  manner 
similar  to  some  of  the  Type  IV,  V  and  VI  flat  plate  specimens. 

In  this  instance,  however,  a  Type  XIVA  member  consists  of  a  plain 
plate  and  a  Type  XIVB  specimen  is  plate  notched  on  the  surface 
with  a  (charpy-V)  notch.  These  results  of  these  various  tests 
will  be  found  in  Figures  4-6  and  4-18. 
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4.3.1  Multiaxial  Stress  Fields 


Tests  of  plain  plates,  butt-welded  plates,  flat  plates  with 
attachments,  and  welded  boxes  will  be  examined  in  detail  to 
obtain  a  general  indication  of  the  effects  of  multiaxial 
stresses  on  the  fatigue  behavior  of  these  types  of  structural 
members . 

Flat  plate  members  will  be  discussed.  Tests  of  Type  I,  Type 
IVB  and  Type  XIVA  members  can  be  examined  to  provide  some 
indication  of  the  effect  of  biaxial  stresses  on  flat  plates. 
The  results  of  fatigue  studies  on  these  three  types  of  members 
are  presented  in  the  following  tabulation.  The  detailed  data 
(range  of  fatigue  stress)  are  presented  in  Table  4-3  and 
Figures  4-2  through  4-6. 


Fatigue  Life  of  HY  Steels 


Fatigue  Strength  at  Life  of 
Type  104  cycles  105  cycles 


I  (HY-80) 

120  ksi 

72  ksi 

IVB*  (HY-80) 

127  ksi 

78  ksi 

IVB*  (HY-130/150) 

194  ksi 

104  ksi 

XIVA  (70B) 

160  ksi 

85  ksi 

*  Fatigue  strength  for  10%  increase  in  deflection. 


From  this  tabulation,  the  axial  fatigue  strength  (Type  I)  is 
slightly  lower  than  the  fatigue  strength  of  the  Type  IVB 
plates  in  bending  for  HY-80  steel.  The  plates  in  bending 
appear  to  have  fatigue  strengths  5  to  8%  higher  than  that  of 
flat  plates  subjected  to  axial  loadings.  However,  when  one 
compares  the  fatigue  resistance  of  the  flat  plates  of  HY-80 
and  HY-130/150  steels,  it  is  found  that  the  fatigue  resistance 
of  the  higher  strength  steel  is  33  to  53%  higher  for  the 
higher  strength  steel.  Thus,  the  strength  of  the  material  has 
a  very  significant  effect  on  the  fatigue  resistance  of  the 
flat  plates,  a  much  greater  effect  than  the  difference  found 
when  comparing  the  bending  and  axial  loading  of  a  given 
material . 

A  comparison  of  the  fatigue  behavior  of  the  unwelded  Type  IVB 
and  Type  XIVA  specimens  gives  an  indication  of  the  effect  of 
one  axis  bending  versus  biaxial  bending  in  a  flat  plate.  In 
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this  instance,  it  is  found  that  the  fatigue  resistance  of  3/4- 
inch  plate  of  the  XIVA  specimens  was  9  to  26%  greater  than 
that  of  the  type  IVB  specimens.  However,  the  materials 
differed  in  thickness,  the  XIVA  specimens  were  3/4-inch  thick 
and  the  IVB  specimens  were  1-1/2-inch  thick  and  they  were  also 
of  a  different  material.  The  IVB  specimens  were  fabricated  of 
HY-80  steel  with  a  tensile  strength  of  108  ksi  whereas  the 
Type  XIVA  members  were  fabricated  of  70B  steel  with  a  tensile 
strength  cf  102  ksi.  Nevertheless,  it  appears  that  flat 
plates  subjected  to  a  transverse  loading  had  a  somewhat  higher 
fatigue  strength  when  subjected  to  biaxial  stressing  than  when 
subjected  to  uniaxial  stressing. 

A  second  group  of  members  examined  consisted  of  flat  plates 
with  a  transverse  butt  weld  which  were  subjected  to  transverse 
uniform  pressure.  These  members,  Type  V,  had  a  lower  fatigue 
strength  than  that  of  the  plain  flat  plates  of  Type  IV, 
previously  considered.  The  S-N  curves  were  also  steeper  thus 
indicating  a  more  severe  notch  effect  in  the  butt  welded 
members  which  produced  the  lower  fatigue  strength. 

In  the  case  of  the  Type  V  members,  both  A  and  B  specimens  were 
tested  and  thereby  give  some  indication  of  the  effect  of  a 
biaxial  loading.  As  indicated  in  Figures  4-7  and  4-8,  the 
biaxial  flexural  loading  again  provided  an  increase  in  fatigue 
resistance.  In  this  instance  the  fatigue  strength  of  the  Type 
VA  members  was  approximately  30%  greater  than  that  of  the  Type 
VB  members. 

The  third  group  of  flat  plate  members  subjected  to  lateral 
loadings  were  plates  with  transverse  welded  attachments 
(Type  VI) .  The  following  tabulation  presents  the  results  of 
the  Types  IV,  V  and  VI  test  specimens. 


Fatigue  Life  of  Different  Details 


Fatigue  Strength  at  Life  of 
Type  104  cycles  10s  cycles 


IV  B 

127  ksi 

78  ksi 

V  B 

106  ksi 

39.5  ksi 

VI  B 

60  ksi 

30  ksi 

V  A 

137  ksi 

51  ksi 

VI  A 

81  ksi 

39  ksi 
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It  is  readily  evident  that  the  effect  of  butt  welding  and  the 
addition  of  welded  attachments  has  a  significant  effect  on  the 
fatigue  resistance  of  the  flat  plates  subjected  to  transverse 
pressures,  the  greater  the  stress  concentration  of  the  weld 
detail,  the  greater  the  reduction  in  fatigue  resistance. 

As  in  the  case  of  detail  Type  V,  it  can  be  seen  that  the 
biaxial  stressing  in  the  Type  VI  members  (Type  VIA  vs.  Type 
VIB)  again  produced  an  increase  in  fatigue  resistance  of 
approximately  30%.  This  is  a  significant  increase  in  fatigue 
resistance  under  the  biaxial  stressing  condition. 

Another  evaluation  can  be  made  by  comparing  the  Type  III,  Type 
VI,  and  Type  VII  specimens.  These  include  plates  subjected  to 
lateral  loading,  beams  subjected  to  bending,  and  axial 
specimens  all  of  which  have  stiffeners  attached  to  them. 

Stress  ratios  R=0  and  R=-l  were  included  for  these  three  types 
of  members.  A  summary  of  the  results  is  presented  in  Figure 
4-19.  Here  it  may  be  seen  that  at  approximately  40,000  cycles 
of  loadings,  the  fatigue  resistance  of  the  three  types  of 
members  are  approximately  equal.  However,  the  Type  III 
specimens  have  a  somewhat  steeper  slope  to  their  S-N  curve. 
This  would  suggest  that  the  stress  concentration  for  the 
detail  was  somewhat  more  severe  for  this  type  of  specimen  than 
for  the  Type  VI  and  VII  specimens.  It  can  also  be  seen  that 
the  fatigue  resistance  is  approximately  50%  greater  under  a 
completely  reversed  stress  cycle  than  it  is  for  a  0  to  tension 
type  of  stress  cycle. 

The  fourth  type  of  member  to  be  considered  is  the  welded  box 
of  HY-100  steel.  These  boxes,  Type  VIII,  were  subjected  to 
internal  pressure  and  consequently  stressed  in  a  manner 
similar  to  the  butt  welded  plates  or  plates  with  attachments, 
Types  VA  and  VIA.  The  stresses  in  the  longitudinal  and 
transverse  direc*  on  of  the  boxes  are  in  a  ratio  of  1  to  2 . 

The  results  for  v..e  various  series  of  tests  on  boxes  are  shown 
in  Figure  4-14.  Their  behavior  appears  to  have  been  quite 
similar  to  that  of  the  stiffened  plates.  Type  VIA.  The  slopes 
of  the  S-N  curves  appear  to  be  quite  similar  and  their  fatigue 
strengths  comparable.  The  plain  butt  welded  plates,  Type  VA, 
appear  to  have  a  somewhat  higher  fatigue  strength.  In  one 
instance,  a  mild  steel  box  (identified  by  M)  was  tested  under 
internal  pressure  and  appeared  to  have  a  somewhat  lower 
fatigue  resistance,  as  shown  in  the  figure. 
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In  summary,  it  appears  that  biaxial  stressing  in  ratio  of 
approximately  1  to  2  provides  a  fatigue  strength  on  the  order 
of  25  to  30%  greater  than  that  of  members  stressed  axially  or 
in  bending  in  only  one  direction. 

4.3.2  Mean  Stress 


The  effect  of  mean  stress  on  the  fatigue  strength  of  the 
various  types  of  specimens  can  most  readily  be  studied  by 
examining  the  fatigue  behavior  of  the  members  at  various 
stress  ratios.  In  Table  4-5,  data  are  available  for  various 
stress  ratios  for  details  of  Type  I,  II,  III,  VII,  and  XI.  A 
summary  of  the  relative  fatigue  strengths  of  the  members  at 
various  stress  ratios  and  at  lives  of  104,  10s,  and  10*  cycles 
are  presented.  In  this  table  it  is  evident  that  a  significant 
difference  exists  in  the  fatigue  strengths  at  stress  ratios  of 
-1,  0,  and  +1/2. 

A  visual  indication  of  the  effect  of  stress  ratio  on  fatigue 
strength  can  be  obtained  from  the  fatigue  diagram  of  Figure  4- 
20.  Here  the  fatigue  strengths  for  the  various  details,  at 
various  lives,  have  been  plotted  and  clearly  show  the  stress 
ranges  decrease  significantly  as  the  mean  stress  is  increased. 
It  appears  that  the  fatigue  resistance  can  be  represented  by 
the  following  relationship. 


where 


Fr 


R 


Fr  =  Fro(l-0. 45R) 

fatigue  stress  range  at  stress  ratio,  R. 

fatigue  stress  range  for  a  stress  ratio  of 
0  (0-to-tension  stress  cycle) . 

stress  ratio,  ratio  of  minimum  to  maximum 
stress  in  the  stress  cycle. 


In  a  similar  study  [4-5],  a  similar  relationship  was  presented 
except  that  the  coefficient  for  R  was  somewhat  smaller. 
Nevertheless,  in  both  evaluations  it  is  seen  that  the  fatigue 
resistance  is  clearly  affected  by  the  mean  stress  to  which  the 
members  are  subjected. 


4.3.3  Stress  Gradient 


Relatively  few  data  are  available  concerning  the  effect  of 
stress  gradient  on  the  fatigue  resistance  of  welded  members. 
Some  data  are  available  on  small  rotating  beam  specimens  and 
are  reported  in  various  locations  in  the  literature.  For 
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TABLE  4-5 


EFFECT  OF 

MEAN 

STRESS  ON 

FATIGUE  STRENGTH  OF  VARIOUS 

MEMBERS 

SPEC. 

TYPE 

STRESS 

RATIO 

FATIGUE  STRENGTH 

AT  VARIOUS  LIVES 

104  105  106 

I 

R=0 

1.0 

1.0 

1.0 

R=-l 

1.30 

1.39 

1.49 

II 

R=+l/2 

0.73 

0.72 

0.71 

R=0 

1.0 

1.0 

1.0 

R=-l 

1.42 

1.42 

1.43 

III 

R=0 

1.0 

1.0 

- 

R=-l 

1.41 

1.43 

- 

VII 

R=0 

1.0 

1.0 

- 

R=-l 

1.45 

1.42 

- 

XI 

R=+l/2 

0.76 

0.74 

0.72 

R=0 

1.0 

1.0 

1.0 

R=-l 

1.33 

1.33 

1.27 

NOTE:  Values  are 

the  ratio 

of  the  fatigue 

strength  at 

the  stress 

ratio 

shown  to  the  fatigue  strength 

at  a  stress 

ratio  R=0. 
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MIN.  STRESS  ksi 

Modified  Goodman  Diagram 


example,  Reference  4-58  indicates  that  as  the  diameter  of  a 
member  increases  the  fatigue  resistance  of  the  material  can  be 
expected  to  decrease  under  a  flexural  loading.  There  are 
indications  that  up  to  a  diameter  or  thickness  of  approximate¬ 
ly  6  inches,  the  fatigue  resistance  continues  to  decrease. 
Beyond  that  thickness  or  depth  there  is  little  or  no  further 
decrease  in  fatigue  resistance.  This  is  in  agreement  with  the 
fatigue  resistance  of  12-in.  deep  beams.  The  fatigue  resis¬ 
tance  of  the  beams  is  approximately  the  same  as  that  of 
axially  loaded  plates  of  the  same  type  of  material,  indicating 
that  there  is  no  stress  gradient  effect  in  the  12 -in.  beam 
tests . 

If  it  is  assumed  that  axially  loaded  members  represent  the 
flexural  minimum  fatigue  resistance  of  deep  members  or  thick 
material  where  there  is  no  stress  gradient  effect,  then  some 
general  indication  of  stress  gradient  effects  can  be  obtained 
by  comparing  the  results  of  the  data  tabulated  in  Table  4-3. 

In  Figure  4-19,  data  are  compared  for  Type  III,  VI,  and  VII 
members.  This  comparison  indicates  that  for  a  life  of 
approximately  40,000  cycles  there  is  no  significant  difference 
in  fatigue  resistance  for  the  axially  loaded  members  (Type  III 
specimens)  and  the  plates  and  beams  1  1/2 -inches  thick  that 
are  loaded  in  flexure  (Type  VI  and  VII  specimens) .  However, 
the  S-N  curves  for  the  Type  III  specimens,  members  with 
attachments  on  both  sides  of  the  plate,  have  a  steeper  slope 
and  consequently  a  lower  fatigue  resistances  at  very  long 
lives. 

A  similar  comparison  is  shown  in  Figure  4-21  for  butt-welded 
plates  under  axial  loading  and  flexural  loading.  Again,  there 
is  a  difference  in  slope  of  the  S-N  curves.  In  this  instance 
however,  the  S-N  curve  for  the  butt-welded  plates  subjected  to 
flexural  loading  has  a  steeper  slope  than  those  for  the 
members  loaded  axially.  It  is  evident  that  much  more 
information  is  required  to  evaluate  the  question  of  stress 
gradient,  particularly  where  multi-axial  stresses  exist. 

4.3.4  Residual  Stresses 

Another  factor  considered  in  this  study  is  the  effect  of 
residual  stresses  on  the  fatigue  behavior  of  the  various  types 
of  test  members.  The  data  presented  in  Figures  4-10  and  4-11 
clearly  show  that  residual  weld  stresses  had  little  or  no 
effect  on  fatigue  life  of  Type  VIB  members;  members  that  were 
stress  relieved  had  approximately  the  same  fatigue  resistance 
as  members  tested  in  the  as-welded  condition.  This  should  be 
expected  for  the  short  lives  studied  in  the  investigations. 
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Comparison  of  Axial  Fatigue  of  Type  II  Specimens 
and  Bending  Fatigue  of  Type  VB  Specimens  of  HY-81 


However,  members  subject  to  lower  stress  ranges  (long  lives) 
stress  relieving  can  be  expected  to  have  an  effect  on  the 
fatigue  resistance  of  welded  test  members.  If  compressive 
residual  stresses  are  introduced  they  can  have  a  beneficial 
fatigue  effect  on  a  member,  even  at  the  high  stress  ranges 
employed  in  the  tests  reported  herein.  This  can  clearly  be 
seen  in  the  Figures  4-22  and  4-23  studies  wherein  test  members 
were  mechanically  peened.  This  mechanical  peening  not  only 
introduced  residual  compressive  stresses  in  the  peened  region 
but  also  changed  the  geometry  at  the  toe  of  the  weld  where 
this  mechanical  peening  was  applied. 

Detailed  examination  of  the  data  in  Reference  4-42,  wherein 
members  were  mechanically  peened  and  then  stress-relieved, 
indicates  that  the  stress-relieving  treatment  reduced  the 
fatigue  resistance  because  it  removed  compressive  residual 
stresses.  In  the  cases  where  a  stress-relieving  treatment  was 
applied  before  mechanical  peening  a  higher  fatigue  resistance 
was  obtained,  indicating  that  the  residual  compressive 
stresses  remained  in  the  member  after  mechanical  peening  in 
this  latter  case  and  the  higher  fatigue  resistance  resulted. 

From  the  limited  data  available,  it  is  evident  that  the  effect 
of  residual  stresses  on  fatigue,  whether  in  uniaxial  or 
multiaxial  test,  depends  upon  the  sense  of  the  residual 
stresses,  their  magnitude,  and  the  magnitude  of  the  applied 
stresses. 

4.3.5  Corrosion 

Although  corrosion  could  be  a  very  important  factor  for  ship 
structures,  relatively  little  information  is  available 
concerning  the  effect  of  salt  water  on  the  fatigue  behavior  of 
weldments  in  ship  steel.  In  Reference  4-22,  three  corrosion 
tests  were  conducted  on  Type  VIB  specimens.  In  these  tests 
salt  water  produced  about  an  18%  reduction  in  fatigue 
resistance  and  approximately  58%  reduction  in  fatigue  life  of 
the  test  members.  However,  since  these  tests  were  conducted 
at  a  frequency  of  20  cycles  per  minute,  the  test  would  only 
have  taken  approximately  50  hours.  This  raises  many  question 
concerning  the  effects  of  corrosion  under  long  life  conditions 
and  also  many  of  the  other  factors  that  are  involved  in 
corrosion  fatigue.  In  Reference 

4-59,  reference  is  made  to  many  pertinent  factors  that  affect 
corrosion  fatigue.  These  include  such  factors  as  temperature, 
oxygen  content,  electrolytic  potential,  load  ratios,  cyclic 
frequency,  etc.  It  is  readily  evident  that  much  more  needs  to 
be  done  to  evaluate  the  effects  of  corrosion  in  fatigue. 
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Cycles  to  Failure 


4.3.6  geometry.  Fabrication  Treatment.  Defects 


It  is  well  known  that  the  geometry  of  welded  structural 
members  has  a  significant  effect  on  the  fatigue  behavior  of 
such  members.  This  is  readily  evident  in  the  fatigue 
strengths  tabulated  in  Table  4-3.  Here  it  can  be  seen  that  at 
a  life  of  105  cycles  fatigue  strengths  for  a  stress  ratio  R=0 
are  found  to  range  from  32  to  54  ksi  and  for  a  life  of  10* 
cycles  the  fatigue  strength  range  from  17.5  ksi  to  57  ksi. 
Among  the  highest  values  are  the  results  from  tests  on  plain 
plate  material  and  cast  members  for  which  the  geometry  is  very 
simple  and  with  little  or  no  stress  concentration. 

As  noted  previously  and  shown  in  Figures  4-22  and  4-23, 
mechanical  peening  the  toes  of  welds  can  improve  the  fatigue 
resistance  of  a  welded  member.  However,  care  must  be 
exercised  in  evaluating  such  improvements  since  much  of  the 
improvement  obtained  in  the  weldments  shown  in  these  figures 
was  a  result  of  the  residual  compressive  stresses  produced  by 
the  mechanical  peening. 

Grinding  of  the  weld  toes  can  also  provide  some  benefit  to  the 
fatigue  resistance.  However,  in  this  instance  the  improvement 
is  a  result  of  the  improved  geometry. 

Another  geometrical  factor  that  may  affect  the  fatigue 
strength  of  structure  weldments  is  their  defect  or 
discontinuity  quality.  In  the  case  of  a  1/16-in.  undercut  in 
Type  VB  and  VIB  specimens  only  a  small  reduction  in  fatigue 
strength  was  observed.  However,  a  lack  of  fusion  was  found  to 
produce  a  greater  reduction  (see  Figures  4-9  through  4-11) . 

See  Reference  4-23. 

4.3.7  size  -  Thickness 

In  the  investigations  included  this  study,  materials  varying 
in  thickness  from  1/2-inch  to  approximately  1  5/8-inches  were 
included.  Nevertheless,  there  appears  to  be  relatively  small 
effect  of  member  thickness.  This  effect  appears  to  be 
overshadowed  by  such  other  factors  as  tensile  strength  and 
magnitude  of  stress  applied  in  the  tests. 

An  examination  of  thickness  in  fact  was  included  for  specimens 
of  Type  I,  VI,  VII,  and  VIII.  In  the  case  of  the  Type  I 
specimens,  plain  flat  plates  subjected  to  axial  loadings, 
thicknesses  ranging  from  1/2-inch  to  1  1/2-inches  were 
included.  In  the  case  of  3/4  and  1  1/2-inch  specimens  of  HY80 
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steel,  greater  lives  were  obtained  for  the  1  1/2-inch  material 
than  for  the  3/4-inch  thick  material  when  tested  on  a  stress 
cycle  0  to  50  ksi.  This  is  contrary  to  the  general 
observation  often  made  that  the  fatigue  resistance  decreases 
with  an  increase  in  thickness.  Upon  examination  of  the 
strength  of  the  two  material,  it  was  found  that  the  tensile 
strength  of  the  thicker  material  was  approximately  8%  greater 
than  the  tensile  strength  of  the 

3/4-inch  thick  material.  This  no  doubt  accounts  for  the 
increased  fatigue  resistance  of  the  thicker  material. 

In  the  case  the  Type  VI  and  Type  VII  specimens,  relatively 
little  information  was  available  concerning  the  effect  of 
thickness  and  no  firm  conclusions  could  be  drawn  concerning 
this  factor. 

For  the  Type  VIII  box  type  specimens,  material  1  inch  and  2 
inches  in  thickness  were  used  in  the  tests.  Neglecting  the 
tensile  strength  of  the  materials  it  is  found  that  both 
thickness  provided  approximately  the  same  strength  for  the 
test  members.  This  may  be  seen  in  Figure  4-14. 

4.3.8  Material  strength 

Material  strengths  ranging  from  mild  steel  to  HY-130/150  steel 
have  been  included  in  this  study.  These  various  materials 
have  been  examined  for  a  variety  of  specimen  type  and  in 
general  it  is  found  that  the  fatigue  strength  tends  to 
increase  with  the  tensile  strength.  However,  the  relationship 
of  the  fatigue  strength  to  tensile  strength  is  affected  by  a 
variety  of  other  factors  such  as  geometry,  stress  level,  etc. 

Combining  the  data  for  flat  plate  specimens  of  HY-80,  HY-100, 
HY-130/150  steels  and  in  thicknesses  of  1/2-inch  to  1  1/2- 
inches,  Reference  4-19,  one  finds  that  the  ratio  of  fatigue 
strength  to  tensile  strength  varies  from  values  of  0.80  to 
0.40,  depending  upon  the  life  of  the  member.  As  shown  in 
Figure  4-15,  the  ratio  of  0.80  is  obtained  for  a  life  of 
20,000  cycles  and  the  value  of  0.40  at  a  life  of  500,000 
cycles.  The  shortest  lives,  of  course,  are  obtained  at  the 
highest  test  stress  levels  and  the  longest  lives  at  the  lower 
test  stress  levels.  Thus,  it  would  appear  that  the  fatigue  to 
tensile  strength  ratio  is  primarily  a  function  of  the 
magnitude  of  stress  to  which  its  members  are  subjected,  this 
being  a  more  important  factor  than  the  strength  of  the 
material  or  the  thickness  of  the  material. 
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The  geometry  of  the  specimens  also  becomes  an  extremely 
important  factor  and  may  overshadow  the  effect  of  the  tensile 
strength  of  the  material.  It  is  evident  that  material  does 
have  some  effect  on  the  behavior  of  the  members  but  other 
factors  such  as  geometry  and  stress  magnitudes  play  a  much 
more  important  role. 

4.3.9  Random  Loading 

While  random  loading  has  been  proven  to  have  a  large  effect  on 
the  fatigue  life  of  structures,  to  date  little  or  no  data  are 
available  concerning  random  loadings  on  members  subjected  to 
multiaxial  fatigue  of  the  type  considered  in  this  study.  The 
type  of  random  loading  distribution  can  have  a  large  influence 
on  the  effective  stress  used  in  the  fatigue  analysis.  We  also 
know  that  marine  structures  are  subjected  to  random  loading 
where  principal  stresses'  proportions  are  also  random.  Random 
loading  is  a  factor  that  must  be  strongly  considered  in  any 
future  multiaxial  fatigue  testing. 

4.3.10  Stress  Phasing 

In  laboratory  studies  of  multiaxial  stresses  the  repeated 
stresses  are  generally  applied  in  a  proportional  manner,  that 
is,  the  proportions  of  flexure  to  shear  is  constant  or  the 
proportion  of  tension  or  flexure  to  torsion  is  constant. 
However,  in  marine  structures  such  as  a  ship  or  an  offshore 
platform,  proportions  will  vary.  Consequently,  in  a  stiffened 
laboratory  beam,  for  example,  the  direction  of  principal 
stress  at  a  given  point  in  the  beam  will  remain  constant  even 
though  the  magnitude  of  the  loading  is  changed.  In  marine 
structures  that  are  subjected  to  varying  flexural  and 
torsional  loads,  the  magnitudes  and  directions  of  the 
principal  stresses  will  change  with  time. 

To  date  there  does  not  appear  to  be  any  information  available 
concerning  the  above  types  of  stress  phasing  in  welded 
structural  details.  However,  considerable  research  has  been 
conducted  on  notched  specimens  as  discussed  in  Section  3.2. 
Until  such  information  becomes  available  for  welded  details  it 
will  probably  be  necessary  to  consider  the  fatigue  of  such 
structural  details  in  terms  of  the  maximum  principal  tensile 
stresses  without  consideration  of  the  oscillation  or  rotation 
in  the  direction  in  which  such  stresses  might  be  applied. 
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5.0  EVALUATION  OF  MULTIAXIAL  FATIGUE  RELIABILITY 


The  evaluation  of  multiaxial  fatigue  reliability  was  performed 
to  quantify  the  effects  of  multiaxial  fatigue  and  with  the 
purpose  of  developing  a  baseline  to  be  used  in  judging  overall 
improvements  in  reliability  through  the  use  of  multiaxial 
fatigue  design  procedures.  The  effort  included  a  review  of 
existing  reliability  formats  and  data  to  develop  a  baseline  to 
judge  overall  improvements  in  reliability  by  using  multiaxial 
fatigue  design  procedures  in  marine  structures. 

5.1  RELIABILITY  FORMAT  FOR  EVALUATING  MULTIAXIAL  FATIGUE  IN 
MARINE  STRUCTURES 

The  purpose  of  engineering  design  is  to  ensure  the  safety  or 
performance  of  a  given  system  for  a  given  period  of  time 
and/or  under  a  specified  loading  condition.  The  absolute 
safety  of  the  system  cannot  be  guaranteed  due  to  the  number  of 
uncertainties  involved.  In  structural  design,  these 
uncertainties  can  be  due  to  randomness  of  loads,  simplifying 
assumptions  in  the  strength  analysis,  material  properties, 
etc.  However,  through  probabilistic  analysis  we  can  limit  the 
risk  of  unacceptable  consequences.  The  major  benefit  of  a 
reliability-based  design  approach  which  utilizes  probabilistic 
analysis  is  that  a  designer  will  be  able  to  generate  an 
engineering  system  which  is  both  efficient  and  reliable  to  the 
level  specified  or  desired. 

The  concepts  of  the  reliability  approaches  for  most  types  of 
structures  are  relatively  simple,  but  complexity  is  associated 
with  execution  of  the  procedure. 

In  general,  the  structural  reliability  problem  can  be 
considered  as  one  of  load  and  resistance.  Failure  occurs  when 
the  supply  (the  resistance  or  strength  of  the  system)  is  less 
than  the  demand  (the  loading  of  the  system) .  For  a  structural 
system,  this  can  be  stated  as: 

Probability  of  Failure  =  =  P  (strength  <  load) 

If  both  the  load  and  strength  are  treated  as  random  variables, 
then  the  reliability  problem  can  be  treated  using 
probabilistic  methods.  In  order  to  perform  a  reliability 
analysis,  a  mathematical  model  which  relates  the  load  and 
resistance  must  be  derived.  This  relationship  is  expressed  in 
the  form  of  a  limit-state  equation.  For  the  simple  case  used 
above,  it  would  appear  as 

g(x)  =  R  -  L 
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where  R  and  L  are  the  resistance  and  load  effect  random 
variables.  Failure  is  represented  by  the  region  where  g(x)  is 
less  than  zero  while  the  safe  region  is  where  g(x)  is  greater 
than  zero.  The  line  g(x)  =  0  represents  the  boundary  between 
the  regions  and  is  thus  defined  as  the  limit-state  equation. 

The  implementation  of  reliability-based  design  methods  does 
not  mean  that  all  engineers  and  designers  need  to  be  deeply 
versed  in  probability  theory.  Rather,  the  design  criteria 
they  used  should  be  developed  in  a  format  which  is  both 
familiar  to  the  users  and  which  should  produce  desired  levels 
of  uniformity  in  safety  among  groups  of  structures.  This 
should  be  accomplished  without  departing  drastically  from 
existing  general  practice. 

One  of  the  more  popular  formats  for  including  probabilistic 
information  in  structural  design  is  the  Load  and  Resistance 
Factor  Design  (LRFD)  approach  as  recommended  by  the  National 
Bureau  of  Standards  [5-1].  This  approach  uses  load 
amplification  factors  and  resistance  reduction  factors 
(partial  safety  factors)  and  can  be  expressed  as: 

n 

0R  >  Z  $,L, 
i=l 


where  R  is  the  resistance,  e.g.  in  flexural  shear,  fatigue, 
etc. ;  L,  is  the  ith  load  effect,  e.g.  due  to  dead,  live,  wind, 
earthquake  loads,  etc. ;  <p  is  the  resistance  reduction  factor; 
6,  is  the  ith  partial  load  effect  amplification  factor;  and  n 
is  the  total  number  of  load  effects  considered  in  the  limit- 
state  design  equation. 

The  implementation  of  an  LRFD  format  for  reliability-based 
fatigue  design  has  been  proposed  by  Albrecht  [5-2],  Ang  [5-3], 
and  Munse  [5-4]  for  cumulative  damage  fatigue  approaches  and 
by  Wirsching  [5-5]  for  the  fracture  mechanics  approach  to 
fatigue  propagation. 

For  fatigue  cumulative  damage  estimates  of  structures,  the 
resistance  is  usually  represented  in  terms  of  the  mean  and 
standard  deviation  of  the  number  of  cycles  to  failure  at  a 
given  stress  range.  This  information  is  typically  derived 
from  constant  amplitude  fatigue  tests  data  for  specimens  or 
components  being  investigated.  According  to  this  approach,  a 
number  of  these  tests  are  conducted  and  the  results  are 
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provided  in  the  form  of  stress  range  versus  life  (S-N)  curves. 
The  question  arises  whether  an  equivalent  stress  or  strain 
approach  based  on  small  component  tests  simulating  the  load 
effects  of  the  complex  joint  can  be  used  to  reduce  the  scatter 
about  the  mean  line,  thus  reduce  conservatism  and  improve 
reliability  of  the  structure. 

Before  this  question  can  be  answered,  it  is  beneficial  to 
review  the  reliability  formats  in  greater  detail  to  illustrate 
the  role  of  fatigue  data  bias  and  scatter  in  the  reliability 
formats. 

The  standard  deviation  of  the  fatigue  life  data  can  be  found. 
However,  the  scatter  of  the  data  about  the  mean  fatigue  line 
is  not  the  only  uncertainty  involved  in  the  S-N  analysis.  A 
measure  of  the  total  uncertainty  (coefficient  of  variation)  in 
fatigue  life,  V„,  is  usually  developed  to  include  the 
uncertainty  in  fatigue  data,  errors  in  the  fatigue  model,  and 
any  uncertainty  in  the  individual  stresses  and  stress  effects. 

Ang  [5-3]  and  Munse  [5-4]  suggested  that  the  total  COV  in 
terms  of  fatigue  life  could  be  given  by: 


V„2  =  VN2  +  VF2  +  Vc2  +  (mVs) 2 

where : 

V„  =  total  COV  of  resistance  in  terms  of  cycles  to 
failure 

VN  =  variation  in  fatigue  test  data  about  mean  S-N 
line 

VF  =  variation  due  to  errors  in  fatigue  model  and 
use  of  Miner's  rule 

Vc  =  variation  due  to  uncertainty  in  mean  intercept 
of  the  regression  line;  includes  effects  of 
fabrication,  workmanship,  and  uncertainty  in 
slope 

Vs  =  variation  due  to  uncertainty  in  equivalent 
stress  range;  includes  effects  of  error  in 
stress  analysis 

m  =  slope  of  mean  S-N  regression  line. 
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Values  of  m  and  VN  can  be  obtained  from  sets  of  S-N  curves  for 
the  type  of  detail  being  investigated;  a  number  of  which  are 
tabulated  by  Munse  [5-4]. 

Reasonable  values  for  the  remaining  uncertainties  are 
available  in  the  literature  [5-3  through  5-8].  Typically  ,  Vs 
is  taken  to  be  0.1,  Vc  is  assumed  to  be  0.4,  and  VF  is  taken  as 
0.15. 


As  can  be  seen  from  the  discussion,  there  are  many  more 
uncertainties  other  than  VN  for  scatter.  However,  this  is 
where  the  primary  improvement  in  reliability  will  be  derived 
when  using  multiaxial  fatigue  design  procedures.  Typically, 
COV  for  scatter  VN  is  .6  for  welded  structural  details. 

To  determine  the  relationship  between  load  and  resistance  for 
fatigue  analysis,  the  value  of  COV  for  scatter  in  fatigue  life 
must  be  converted  to  total  COV  of  resistance  in  terms  of 
stress  range.  Using  the  properties  of  log-normal  distribu¬ 
tions,  resistance  transformation  in  terms  of  stress  range  is 
given  by 
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Reiterating  the  reliability  equation  where 

0R  '  >  Y  L ' 

where  the  load  and  resistance  factors,  Y  and  <p,  respectively, 
are  expressed  in  terms  of  stress  range  and  can  be  given  by 

Y  =  exp  ( 5BVJ ) 

0  =  exp  (-$BV„ ' ) 

where  S  =  a  splitting  function  and  B  is  an  acceptable 
probability  of  failure  [5-9]. 

The  reliability  equation  can  also  be  used  to  establish  a  level 
of  reliability  for  fracture  mechanics  fatigue  crack  growth 
methodologies.  Again,  there  are  numerous  uncertainties 
including  initial  flaw  size,  type  and  orientation,  effects  of 
experimental  procedures,  and  geometry  of  structure  to  name  a 
few.  These,  of  course,  contribute  to  the  bias  and  scatter  of 
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the  crack  growth  rate  parameters  on  the  Paris  equation 
discussed  in  Section  2. 

5.2  RELIABILITY  07  MULT I AXIAL  FATIGUE  RESEARCH 

With  this  background  we  now  have  a  basis  for  comparing 
multiaxial  fatigue  data  to  existing  approaches  and  determine 
if  there  are  improvements  in  bias  and  scatter  resulting  in 
more  reliable  designs. 

The  British  Department  of  Energy  (DOE)  [5-10]  published 
fatigue  design  rules  for  offshore  structures.  Reliability 
data  is  presented  for  T-joint  and  S-N  curve  and  is  said  to  be 
applicable  to  other  complex  joints  as  well.  The  standard 
deviation  of  special  interest  is  the  use  of  principal  stress 
and  maximum  shear  stress  at  the  hot  spot  rather  than  stress 
normal  to  the  weld  toe  as  applicable  to  the  API  rules  [5-11] 
where  the  standard  deviation  is  on  the  order  of  2  or  more. 

The  resulting  improvement  in  standard  deviation  is  apparent  by 
using  a  combined  stress. 

Nokleby  [5-12]  reviewed  criteria  for  prediction  of  high  cycle 
fatigue  under  multiaxial  stress  conditions  including 
octahedral  stress  criterion,  shear  stress  criterion,  internal 
friction  criterion  and  shear  stress  intensity  criterion.  The 
tests  were  conducted  on  smooth  steel  specimens  with  the 
results  as  follows. 


Criterion 

Bias 

COV 

Octahedral  stress 

1.021 

9.2 

Shear  stress 

1.080 

9.5 

Internal  friction 

1.023 

6.8 

Shear  stress  intensity 

1.018 

6.4 

A  positive  bias  indicates  a  conservative  prediction  and 
measurement  value.  From  the  paper,  it  appears  that  the  bias 
of  all  four  methods  are  on  the  conservative  side  most  of  the 
time. 

Munse  [5-4]  investigated  the  fatigue  initiation  and  crack 
growth  in  stiffened  plate  beam  similar  to  CVK  or  web  frames  in 
ship  structures.  An  S-N  curve  is  presented  in  which  the 
maximum  bending  stress  on  the  extreme  fiber  at  the  location  of 
failure  is  used  to  correlate  data.  These  data  are  shown  in 
Figure  3-27.  It  is  apparent  from  this  comparison  that  the 
maximum  bending  stress  at  the  fracture  section  does  not 
provide  a  consistent  relationship  in  these  tests  either.  A 
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much  better  correlation  is  obtained  when  the  data  are  analyzed 
on  the  basis  of  the  maximum  principal  tensile  stress 
(including  the  effect  of  shear)  at  the  point  of  failure.  Data 
plotted  on  this  latter  basis  are  presented  in  Figure  3-28. 
Although  there  is  still  scatter  in  the  test  results,  there  is 
a  well-defined  scatter  band  for  the  entire  range  of  lives. 

It  appears  that  the  shear  and  flexural  stress  in  the  web  of 
the  welded  members  and  the  resulting  principal  tensile 
stresses,  along  with  stress  concentration  gradients  at 
stiffener  welds,  produce  the  most  critical  fatigue  condition 
in  the  weldments  with  various  stiffeners. 

Moyar  and  Garg  [5-13]  re-analyzed  Munse's  plate  beam  data  and 
correlated  the  data  using  maximum  principal  stress,  modified 
maximum  range  of  shear  stress,  and  the  maximum  octahedral 
shear  stress  criteria. 

A  comparison  of  the  corresponding  predictions  for  the  three 
theories  versus  the  observed  cycles  to  failure  for  the  type  B 
beam  specimens  is  provided  in  Table  3-4.  The  mean  and 
standard  deviation  of  the  difference  of  predicted  minus 
observed  cycles  to  failure  are  presented.  The  maximum  shear 
range  theory  appears  more  satisfactory  in  that  the  mean 
difference  in  prediction  and  observation  is  smaller.  Overall, 
the  maximum  principal  stress  theory  is  the  least  conservative. 

Yung  and  Lawrence  [5-14]  evaluated  multiaxial  fatigue  design 
procedures  for  tube  plate  welds  under  combined  bending  and 
torsion.  The  multiaxial  fatigue  theories  of  Lohr  and  Ellison 
and  Kandil,  Brown  and  Miller  gave  reasonable  non-conservative 
estimates  of  fatigue  life.  Table  3-5  presents  the  fatigue 
life  data  and  results  of  predictions.  As  indicated,  the  mean 
values  for  maximum  shear  strain  are  in  best  agreement  with 
experimental  results.  The  COV  was  not  calculated,  however 
Figures  3-31  through  3-33  illustrate  the  scatter  in  data. 
Although  there  are  differences  in  bias  for  each  approach,  the 
scatter  is  similar  for  the  limited  number  of  points  shown. 

Lawrence  et  al  [5-15]  presented  reliability  data  for  their 
Initiation-Propagation  approach  for  predicting  fatigue  life  of 
welded  details.  He  compared  his  approach  with  that  of  Munse's 
and  others.  The  comparison  is  very  interesting  in  that  it 
shows  the  potential  gains  in  predicting  the  stress 
distribution  at  the  weld  toe  versus  nominal  stress  in  the 
structural  member  as  is  done  in  Munse's  approach.  However, 
there  is  a  key  tradeoff  between  the  additional  work  required 
to  predict  that  combined  stress  at  the  weld  toe  and  having 
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stress  data  presented  as  nominal  stress  with  additional 
inherent  data  scatter.  Also,  Lawrence  had  the  advantage  of 
knowing  the  notch  stress  concentration  factors  of  his  welded 
specimens  —  an  item  not  known  to  structural  designers.  This 
issue  is  central  to  the  need  for  multiaxial  fatigue  research 
and  the  possible  gain  in  reliability  versus  ease  of  use  for 
designers. 

5.3  EVALUATION  OF  MULTIAXIAL  FATIGUE  RELIABILITY  FOR  MARINE 
STRUCTURES 

To  put  these  coefficients  of  variance  data  in  perspective,  it 
is  helpful  to  estimate  the  effect  on  total  reliability  or 
safety  factor  for  gains  in  reliability  from  using  multiaxial 
techniques.  The  LRFD  format  presented  by  Noland  and  Albrecht 
[5-7]  is  used  to  makr  this  example  relatively  simple  and  easy 
to  follow  rather  than  using  more  complex  variations.  The  net 
input  will  be  similar  to  other  techniques  by  using 

Vs  —  .1 

Vc  =  .4 

VF  =  .15 

and  varying  VN  (scatter  of  S-N  data)  . 

The  impact  of  total  reliability  becomes  clear.  As  shown  in 
Table  5-1,  large  improvements  in  coefficient  of  variance  in 
data  scatter  produce  relatively  small  changes  in  total 
reliability.  This  small  change  results  because  the  magnitude 
of  COV  for  the  other  factors  or  uncertainties  in  the  total 
problem  of  fatigue  in  marine  structures.  In  other  words, 
there  are  many  large  uncertainties  in  the  overall  problem.  A 
large  improvement  in  one  parameter  only  produces  a  small  net 
effect.  This  finding  strongly  suggests  improvements  should  be 
made  across  the  board  rather  than  singly.  This  conclusion  is 
not  completely  fair  because  if  a  multiaxial  fatigue  procedure 
were  used  it  is  likely  that  other  uncertainties  would  be 
improved.  For  example,  the  level  of  accuracy  in  stress 
prediction  is  likely  to  be  greater,  helping  to  improve  the 
overall  level  of  reliability. 

Nevertheless,  this  is  a  good  illustration  of  the  potential 
overall  benefit  if  multiaxial  fatigue  procedures  are  used. 

On  a  more  deterministic  note,  one  of  the  original  questions 
often  asked  by  structural  designers  is  how  to  characterize  the 
complex  stress  distribution  for  fatigue  analysis.  With 
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TABLE  5-1 

POTENTIAL  IMPROVEMENTS  IN  OVERALL  RELIABILITY 
USING  MULT I AXIAL  FATIGUE  PROCEDURES 

%  of  %  of 


vn 

v*2 

v„' 

0/Y 

API 

Munse/ DOE 

MAF 

Munse 

.25* 

.385 

.143 

.723 

2.07 

33 

10 

/DOE 

.6 

.703 

.184 

.653 

2.3 

20 

** 

API 

1.3 

2.033 

.27 

.545 

2.75 

*  Representative  value  based  on  data  presented  for  welded  plate  beams 
and  welded  tube  plate  connections. 
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implementation  of  the  research  program  described  next,  these 
answers  will  be  provided. 
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6.0  CONCLUSIONS  AMD  RECOMMEND AT IONS 


General  multiaxial  fatigue  research  has  been  investigated  for 
simple  welded  structures  by  several  researchers.  There  are 
indications  that  multiaxial  fatigue  approaches  are  capable  of 
improving  fatigue  design  in  marine  structures.  However, 
additional  effort  is  required  to  fully  characterize  multiaxial 
fatigue  response  in  marine  structures  to  a  level  of  accuracy 
beyond  existing  techniques  that  are  based  on  component  tests  and 
nominal  stress  fields.  The  vast  majority  of  multiaxial  fatigue 
design  procedures  were  developed  for  unwelded  structures  where 
structural  loading  is  not  as  complex  as  for  marine  structures 
where  principal  stresses  and  mean  stresses  are  random  —  a  factor 
unaccounted  for  in  most  multiaxial  fatigue  procedures. 

Implementation  of  this  technology  should  certainly  not  be  limited 
to  component  testing  arena.  There  will  be  many  "spin-offs"  from 
the  advances  made  or  anticipated  in  multiaxial  fatigue  research. 
For  example,  engineering  design  relies  heavily  on  analytical 
procedures  such  as  finite  element  analyses  to  aid  the  design  and 
testing  procedures.  The  results  of  these  analyses  are  often 
quite  detailed  and  may  be  the  input  to  fatigue  life  predictions 
made  prior  to  construction.  Redefining  the  critical  parameters 
necessary  for  fatigue  life  predictability  will  certainly  affect  a 
large  base  of  existing  software,  not  to  mention  the  actual 
analytical  procedures.  The  designs  that  are  complex  enough  to 
require  finite  element  modeling  may  be  the  very  applications  that 
require  a  fatigue  life  predictive  approach  using  other  than  the 
existing  component-based  criteria. 

Additionally,  the  approach  must  be  used  in  conjunction  with 
predictive  techniques  that  are  capable  of  predicting  stress 
characteristics  in  the  weld  regions  in  great  detail.  While  there 
is  a  trend  toward  detail  in  stress  predictions,  there  is 
considerable  additional  effort.  The  designer  must  choose  the 
approach  which  best  fits  the  design  application  and  stress 
conditions.  Similarly,  when  detailed  fitness  for  service 
assessments  are  conducted,  these  multiaxial  approaches  will  be 
quite  useful,  especially  fracture  mechanics  approaches  that 
account  for  mixed  mode  Stage  I  crack  growth. 

6 . 1  CONCLUSIONS 

Conclusions  on  the  multiaxial  fatigue  approaches  are  outlined 
below.  These  conclusions  are  presented  to  summarize  the  findings 
cited  above  and  to  form  a  basis  for  further  research 
recommendations . 

1.  There  is  currently  no  universal  parameter  for  correlating 
cyclic  multiaxial  stress/strain  with  fatigue  life  for 
marine  structures.  Very  few  have  been  investigated  for 
welded  joints  as  a  group  and  additional  validation 
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efforts  are  required  before  they  can  be  recommended  for 
general  application  to  marine  structures. 

2.  Of  the  correlation  approaches  reviewed,  the  effective  or 
equivalent  criteria  such  as  use  of  maximum  shear  stress 
for  crack  initiation  and  maximum  principal  stress  for 
crack  growth  are  found  to  be  potentially  useful  tools  for 
extrapolating  materials'  responses  from  one  stress  state 
to  another  and  for  life  correlation  in  high  cycle 
multiaxial  regimes  and  appear  to  be  suitable  for 
predicting  fatigue  initiation  and  propagation  in  simple 
weld  configurations  (e.g.,  butt  welds  in  laterally  loaded 
plates) .  However,  additional  research  is  required  before 
these  approaches  may  be  integrated  into  a  design 
approach . 

3.  There  are  a  number  of  factors  whose  influence  must  be 
quantified  and  validated  before  multiaxial  techniques  can 
be  applied  to  welded  details  with  acceptable  levels  of 
confidence.  These  factors  include:  mean  stress,  residual 
stress,  stress  gradients,  and  random  load  and  phasing 
effects. 

4.  There  is  currently  insufficient  data  to  determine 
conclusively  if  there  is  a  significant  improvement  in 
fatigue  reliability  (bias  and  scatter)  over  techniques 
such  as  Munse's  approach  presented  in  SSC-318  or  the  UK 
DOE  code  for  offshore  structures.  By  characterizing 
detailed  stress  distributions  on  a  local  level, 
additional  variables  are  added  that  are  unquantified  and 
increase  uncertainty  in  overall  levels  of  reliability 
even  though  the  state  of  stress  may  be  known  on  a  detail 
level.  From  a  deterministic  perspective  the  equivalent 
stress  correlations  (e.g.  principal  stress)  are  useful 
for  structural  analysts  interested  in  predicting  stresses 
in  welded  details;  however,  many  other  variables  must  be 
considered  in  any  multiaxial  fatigue  analysis. 

5.  The  stress  distributions  in  ship  structural  details 
examined  as  part  of  this  study  are  located  in  stress 
fields  where  axial  primary  stress  dominates  the  principal 
stress  field.  There  are  secondary  biaxial  stress  effects 
at  frames  and  floors.  Cutouts  produce  axial  stress  at 
the  free  surface  in  the  global  biaxial  stress  fields. 

The  principal  stress  phase  relationship  is  random.  On 
the  other  extreme,  offshore  structures  have  complex 
tubular  joint  configurations  subjected  to  multiaxial 
stresses  and  strains  with  varying  out  of  phase  principal 
stress  fields. 
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6.  Both  ships  and  offshore  platforms  experience  both  a 
constant  mean  stress  (e.g.  dead  load  from  structural 
weight)  and  varying  mean  stresses  that  have  not  been 
accounted  for  in  existing  fatigue  design  procedures,  let 
alone  for  multiaxial  approaches. 

7.  Eventually,  more  extensive  characterization  of  time 
dependent  multiaxial  loading  spectra  of  actual  marine 
structural  components  will  be  useful  in  directing  the 
application  oriented  research  efforts. 


6.2  RECOMMENDED  MULTIAXIAL  FATIGUE  RE8EARCB 

Designers  of  marine  structures  need  information  on  the  fatigue 
response  of  complex  welded  details  under  loading  conditions  with 
variable  time-dependent  principal  stress.  Based  on  the  review  of 
multiaxial  fatigue  it  is  concluded  that,  while  there  are 
potentially  useful  tools  available,  none  have  been  developed  or 
validated  to  an  extent  that  will  provide  a  reliable  design 
without  extensive  component  testing  or  analysis. 

Although  multiaxial  fatigue  technology  has  been  developed  for 
relatively  simple  structures  this  technology  will  not  transfer 
directly  to  complex  marine  structures.  Additional  research  is 
required  to  validate  multiaxial  fatigue  technology  for  marine 
structures.  The  first  step  in  the  research  program  is  to  define 
the  stress  parameters,  primarily  in  ship  structures;  to  conduct 
follow  on  tests  based  on  realistic  stress  conditions. 


6.2.1.  DEFINE  SPACIAL  AND  TEMPORAL  CHARACTERISTICS  OF  PRINCIPAL 
STRESSES  IN  SHIP  AND  OFFSHORE  WELDED  DETAIL8. 

Principal  stress  characteristics  and  their  phase  relationships 
have  been  shown  to  have  a  substantial  influence  on  multiaxial 
fatigue  response  and  the  selection  of  appropriate  prediction 
technique.  These  characteristics  must  be  fully  defined  before 
multiaxial  fatigue  technology  can  be  validated  for  use  in  marine 
structures. 

The  characteristics  of  principal  stresses  include  magnitudes, 
phase  relationships,  gradients  and  mean  stress  components.  The 
random  nature  of  each  shall  be  considered.  These 
characteristics  shall  be  identified  for  typical  structural 
details  including  a  bottom  plate  butt  weld,  a  CVK  bracket  and  a 
longitudinal  cut  out  for  ships  and  multiplaner  KT  joints  with 
overlapping  and  non-overlapping  members  for  offshore  structures. 
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Typical  loading  conditions  shall  be  identified  and  applied  to  the 
marine  structure  and  detailed  stress  calculations  shall  be 
performed  to  a  level  of  detail  required  to  support  nominal  stress 
and  crack  growth  studies. 

6.2.2  CONDUCT  MULTIAXIAL  FATIGUE  TESTS  ON  MARINE  STRUCTURAL 
DETAILS 

Multiaxial  fatigue  testing  should  be  conducted  to  validate 
approaches  using  representative  stresses  and  strains  identified 
in  the  previous  project. 

This  effort  should  be  conducted  to  validate  multiaxial  fatigue 
approaches  in  marine  structures  based  on  their  unique  loading 
characteristics.  The  best  approach  will  likely  depend  on  the 
specific  structure  considered.  Even  in  marine  structures  as  a 
group,  principal  stress  characteristics  are  different.  The 
associated  multiaxial  fatigue  approach  will  reflect  the  stress 
characteristics  identified.  Fatigue  initiation  and  crack  growth 
tests  shall  be  conducted.  Sufficient  tests  shall  be  conducted  to 
estimate  bias  and  scatter  for  integration  into  reliability 
studies.  Analytical  correlation  shall  be  conducted  to  validate 
the  various  multiaxial  fatigue  techniques  that  apply  to  marine 
structures  including  stress  and  strain  based  fatigue  life 
criteria  and  crack  growth  criteria  (e.g.  K.rt) .  Studies  should 
concentrate  on  the  high  cycle  fatigue  range. 


6.2.3  INTEGRATE  MULTIAXIAL  FATIGUE  RESEARCH  INTO  A  RELIABILITY 
BASED  FORMAT 

The  objective  of  this  research  is  to  develop  analytical 
approaches  for  predicting  multiaxial  fatigue  response  in  marine 
structures.  The  approaches  will  likely  be  different  depending  on 
the  marine  structural  detail.  A  decision  tree  or  expert  system 
may  be  appropriate.  The  tasks  for  this  project  should  be 
developed  as  part  of  the  project  described  above. 
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EXAMPLES  OF  MULTIAXIAL  FATIGUE  LIFE  PREDICTIONS  FOR  A 
LATERALLY  LOADED  PLATE  WITH  A  TRANSVERSE  BUTT  WELD 

In  this  section  several  examples  will  be  developed  using  both 
uniaxial  and  multiaxial  approaches  to  illustrate  and  compare  the 
methodologies.  The  example  case  is  based  on  tests  of  HY-80  test 
specimens,  each  28"  x  56"  x  1.50"  thick  welded  on  the  centerline 
of  the  long  dimension.  The  test  specimens  were  simply  supported 
(on  the  upper  surface)  and  subjected  to  a  uniform  pressure 
loading  (on  the  lower  surface)  from  zero  to  maximum  pressure 
which  yielded  a  maximum  calculated  stress  at  the  center  of  the 
panel  of  80  ksi.  The  tests  were  conducted  at  the  Naval  Applied 
Science  Laboratory,  Brooklyn,  New  York  in  1968.  Results  are 
reported  in  Reference  [B-l]..  A  sketch  of  the  panel  is  shown  in 
Figure  B-l.  A  summary  of  the  results  is  given  in  Table  B-l.  The 
methodologies  used  for  comparison  of  stress-based  methods  are: 

(a)  maximum  principal  stress 

(b)  maximum  shearing  stress 

(c)  maximum  octahedral  shearing  stress 

The  principal  stresses  at  the  center  of  an  unwelded  plate  of  the 
dimensions  shown  in  Figure  B-l  have  been  calculated  by  the 
methods  of  Timonshenko  and  Woinowsky-Krieger  [B-2].  The  results 
for  all  four  test  articles  are: 

a,  =  ax  =  80,000  psi 

o2  =  ay  =  36,520  psi 

Tm.x  =  40,000  psi 

The  objective  is  to  predict  fatigue  life  for  the  butt  welded 
plate  using  simple  uniaxial  test  data.  The  S-N  curve  for 
unwelded  HY-80  is  plotted  on  the  basis  of  data  given  for  Type  I 
specimens.  These  specimens  are  standard  flat  plate  uniaxial 
fatigue  test  specimens.  The  data  is  plotted  in  Figure  B-2. 

Each  methodology  will  now  be  outlined  briefly  and  the  results 
compared  with  the  lives  shown  in  Table  B-l.  The  means  of  the 
observed  lives  shown  there  ate: 

Fatigue  life  to  first  observation  of  cracking, 

N/,  =  13,950  cycles 

Fatigue  life  to  10%  increase  in  deflection, 

N/2  =  36,138  cycles. 
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(a)  Maximum  principal  stress  This  method  assumes  that 
fatigue  lives  will  correlate  when  the  maximum  principal  stress, 
corrected  for  the  presence  of  stress  concentrations,  in  the 
structure  is  the  same  as  that  in  the  uniaxial  tests  at  the  point 
of  failure. 

The  fatigue  notch  factor,  Kf  of  Reference  B-3,  was  used  to 
account  for  stress  concentration  due  to  the  longitudinal  weld. 

Kt-1 

Kf  -  1  +  - 

1+a/r 

where 

Kt  =  6  [l  +  a  (t/r)2] 

From  Reference  B-3, 

x  =  0.165  (tan0)  167  for  butt  welds  in  bending 

B  =  1.0 

a  =  0.50 

From  Reference  B-3 , 

<t>m  =  18.25°  and  (tan0)m 167  =  .8274 

Assume  the  notch  root  radius,  r  =  .025  in.  and  that  a,  the  crack 
length,  is  also  .025  in.  Then  Kt  =  2.057  and 

Kf  =  1.529 

The  stress  at  the  midpoint  of  the  panel  and  adjacent  to  the  weld 
will  be 

Of  ~  (1.529) (80,000)  =  122,320  psi 

The  pressure  loading  on  the  panel  varies  from  zero  to  maximum, 
thus  enter  the  S-N  curve  for  R  =  0  with  of  =  122,320  psi  and  find 
that  N/  =  9,400  cycles. 

(b)  Maximum  shearing  stress.  The  maximum  shearing  stress 
criteria  assumes  that  failure  in  the  structure  will  occur  at  the 
same  life  that  would  be  observed  in  a  uniaxial  fatigue  test  at 
the  same  level  of  maximum  shearing  stress. 
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Strain  amplitude  (loq  scale) 


For  the  plate  and  for  the  uniaxial  test  specimen, 


r 


max 


o , 
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Applying  this  relationship  to  the  R=0  S-N  curve  gives  the  curve 
labeled  "r,,,"  in  Figure  B-2.  Using  the  same  value  of  K /  found  in 
(a)  above  yields 

Tmx  =  (1.529) (40,000)  =  61,160  psi 

In  this  case  the  same  life  is  obtained  as  in  case  (a) 

N /  =  9,400  cycles 

(c)  Maximum  octahedral  shearing  stress.  This  criterion 
assumes  that  failure  will  occur  at  the  same  lives  when  the 
maximum  octahedral  stress  in  the  structure  is  the  same  as  that  in 
the  uniaxial  fatigue  test.  Octahedral  shearing  stress  is  given 
by 


roct  =  -  \i°-Oz)Z+  (02-03)2  +  (C.-C^j  1/2 

For  the  uniaxial  test  specimen  a2  =  a3  -  0  this  reduces  to 

n 

To ««  =  -  <*i  =  *471  a, 

3 

Applying  this  constant  to  the  S-N  curve  yields  the  curve  labeled 
"Toct"  in  Figure  B-2. 

For  the  plate  =  80,000  psi,  o2  -  36,520  psi,  o3  =  0  which 
yields 

Toet  =  32,710  psi 

Applying  the  value  of  found  in  (a)  to  roct,  the  fatigue  life  is 
=  18,000  cycles. 
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A  comparison  of  methods  is  presented  below. 


Equivalent  N/, 


Stress,  ksi 

Cycles 

Maximum  principal  stress 

122. 

9,400 

Maximum  shearing  stress 

61. 

9,400 

Maximum  octahedral  shearing  stress 

32. 

18,000 

No.  of  cycles  to  initiation  from  tests 

- 

13,950 

Thus  it  appears  that  the  shear  stress  and  principal  stress 
provides  the  best  albeit  slightly  conservative  predictions  when 
using  uniaxial  data  to  predict  biaxial  fatigue  in  a  laterally 
loaded  plate  with  butt  weld.  Octahedral  stress  is  closest  in 
magnitude  but  under  predicts  fatigue  life,  an  undesirable  outcome 
from  a  designers  viewpoint.  The  effects  of  stress  concentration 
factors  for  weld  details  has  a  large  influence  in  these 
predictions  and  must  be  considered  carefully  in  prediction  of 
multiaxial  fatigue  response  in  welded  details. 

Strain-based  methods 


The  original  intent  of  this  section  was  to  make  a  direct 
comparison  of  stress-based  methods  and  strain-based  methods  using 
the  results  of  large  scale  fatigue  tests  on  a  shipbuilding  steel, 
thus  the  choice  of  the  Naval  Applied  Science  Laboratory  (NASL) 
tests  of  a  56"  by  28"  HY-80  butt-welded  specimen.  The  methods  to 
be  considered  were: 


(a)  maximum  principal  strain 

(b)  maximum  shear  strain 

(c)  Lohr  and  Ellison  critical  plane  method 

(d)  Kandil,  Brown  and  Miller  critical  plane  method. 

Each  of  the  above  involves  a  variant  of  the  strain-life  equation: 

A  €  a,' 

—  =  -  (2N)b  +  6,  (2N)C 

2  E 


A  e 

where  —  =  total  strain  amplitude 

2 


af  '  =  fatigue  strength  coefficient 

b  =  fatigue  strength  exponent 
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cf'  =  fatigue  ductility  coefficient 

c  =  fatigue  ductility  exponent 

2N  =  number  of  stress  reversals. 


A  graphical  interpretation  of  the  cyclic  constraints  af\  ef',  b 
and  c  is  shown  in  Figure  B-3  taken  from  Reference  B-5.  The  same 
reference  contains  a  summary  of  the  monotonic  and  cyclic  strain 
properties  of  selected  engineering  alloys  —  a  list  that  is  about 
as  complete  as  any  available  from  other  public  domain  sources. 
While  the  properties  of  a  number  of  steels  are  included,  HY-80  is 
not  among  them,  and  therein  lies  one  of  the  problems  in  applying 
strain  methods  —  lack  of  a  comprehensive  and  complete  tabulation 
of  multiaxial  fatigue  properties.  In  the  absence  of  reliable 
information  on  these  properties  for  HY-80  it  was  not  possible  to 
proceed  further  in  the  comparison  of  stress-based  with  strain- 
based  biaxial  life  predictions. 
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